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Some Techniques for Measuring Small Mutual Inductances 
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A method of measuring small mutual inductances is presented. 


was 0.1 wH. 
detail. 


The smallest inductor measured 


The circuit is a transformer-ratio-arm bridge with multiple balances and is described in 
Uncertainties are of the order of one part in 107. 


Several ideas for the design and construction of suitable mutual inductance standards are presented. 


Key Words: Auxiliary generator, bridge, coaxial choke, mutual inductor, stray magnetic field, 


transformer. 


1. Introduction 


A recent theoretical study by Page [1]! suggests 
the desirability of making accurate measurements of 
mutual inductance at the 0.1 wH levei and of stepping 
up to larger values. An experimental bridge has been 
constructed for the purpose of investigating some tech- 
niques which seemed appropriate in this range. 
Some of those techniques which appear to be particu- 
larly useful are described in this paper. The bridge 
is a transformer-ratio-arm bridge, having a ratio of 
10 to 1, and has been designed to compare mutual 
inductors through five orders of magnitude, from 0.1 
HH to 1 mH at a frequency of 1592 Hz. Sensitivity 
of a part in 107 is achieved at the 0.1 wH to 1 wH level 
provided the inductors carry 3A in their primary 
windings. 


2. Bridge Circuit 


A simplified version of the bridge is shown in figure 1. 
Mutual impedances Z, and Z2, which are both mainly 
inductive, are defined by the equations 


Z,=joM, +R, 
Z2=jwM.+R» 


If both inductors have the same current, i, in their 
primary windings,? the comparison of Z,; to Z»2 is 
achieved by measuring the ratio of the induced voltages 
e; and e. Since Z,=2 and Z,=2, then #3. 

i i Zo e2 
The ratio of e;/e2 is compared with the 10 to 1 voltage 
ratio of an inductive voltage divider, which may be 
calibrated by a modification of the capacitance-ratio 
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Colo. 80302 
' Figures in brackets indicate the literature references at the end of this paper. 


address: High Frequency Impedance Standards Section. 
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2 Ground capacitance current between the primary windings is negligible; it is discussed 


in section 4 of this paper 


method described by Cutkosky and Shields [2]. The 
bridge balance is accomplished by the adjustment of 
e3, which is an adjustable voltage source, small relative 
to e; and eg. 

The detector consists of a tunable amplifier and a 
tuned high-impedance preamplifier, coupled to the 
bridge through an impedance matching transformer. 
The transformer is tuned with a variable capacitor 
to give maximum output voltage. 

Accurate measurement of small inductances re- 
quires that all sources of stray magnetic field be well 
controlled. To that end, inductors and transformers 
are toroidal, bridge components are enclosed in copper 
eddy-current shields, components are interconnected 
by coaxial cables, and shield current of the coaxial 
leads is restricted to be nearly equal and opposite 
to current in their inner conductors. Figure 2 illus- 
trates shielding and the use of the technique of thread- 
ing a coaxial lead through a high permeability core 
[3] forming a coaxial choke to equalize current in 
the inner conductor and outer shield of the coaxial 
cable. 

Following construction of the bridge, tests were 
made to establish the effectiveness of the design in 
keeping stray magnetic fields small. Such tests were 
made by using two coils of wire, one connected to a 

















FIGURE 1. Simplified inductance bridge circuit. 
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FIGURE 2. Inductance showing shielding, junction points (A, B, 


C, D, E, F) of the mutual inductor, and use of high-permeability 
core to equalize coaxial current. 














detector and used to find sources of stray field, and 
the other connected to a generator and used to find 
places in the circuit that were sensitive to stray field. 
For example, the detector coil was moved around 
until a source of stray field was found, and the amount 
of detector off-balance was noted. The coil connected 
to the generator was driven so that its field gave about 
the same off-balance ai a comparable distance from 
the detector coil. Then the driven coil was used as a 
probe to determine the effect of such a field on the 
main bridge balance waen brought near various 
bridge components. 

The largest sources of — field were the coaxial 
leads in the loop that includes the inductor primary 
windings since the current is largest in that part of 
the circuit. It was determined that the coaxial choke, 
mentioned above, was not sufficient to equalize cur- 
rent in the inner conductor and outer shield of the 
coaxial cable. A greatly improved coaxial choke 
was realized by making it an active device using two 
cores in a manner described by Gibbings [4]. As it 
is used in this application, it may be called an “active 
coaxial choke.” The principle is shown in figure 3. 
If ep=Zpi, eg=Gep, and er=ep+eg, it can be shown 
that Zr7=Zp(1+G). Thus, the active choke reduces 
net current much more than a coaxial choke com- 
posed of a single magnetic core. The gain of the 
amplifier and the number of turns of cable through 
the device was adjusted to reduce net current by an 
amount which made the error negligible. 

Small stray magnetic fields due to slight nonsym- 
metries in the coaxial cables remained. Using the 
probe technique discussed above, three modifications 
were made. First, the copper can serving as an eddy- 
current shield for the inductive divider was found to 
be unsatisfactory because the lid was electrically 
attached to the can at only one point. The shield was 
improved by bolting the lid to the’ can firmly with 
screws spaced about 2 in apart around its circum- 
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FIGURE 3. ‘Active coaxial choke. 


ference. Second, an uncertainty in the value of mutual 
inductance arises from the coupling between the leads 
connected to the inductor primary windings and the 
leads connected to the inductor secondary windings, 
if the leads move relative to each other. Subsequent 
inductors were designed to have the primary winding 
terminals on the opposite side of the shield can from 
the secondary terminals. Third, the smallest induc- 
tors (for which the ratio of mutual inductance of the 
leads to mutual inductance of the inductor is largest) 
were modified as follows: The eddy-current shield 
was extended over the secondary leads by passing 
these leads through a three-inch-long copper cylinder 
having ¥s-in wall thickness (about two skin depths at 
the frequency used). 

The standards to be compared, as illustrated in 
figure 2, are defined as follows: Let i; be the current 
into the inner terminal of connector A, i; be the cur- 
rent out of the outer terminal, and e, be the voltage 
from the inner to the outer terminal. Currents and 
voltages are similarly defined at connectors B, C, D, 
E,and F. Then the mutual impedances of interest are 
defined by Zan=F 

B 

i,=t,=0: Zn =2 

ly 

ex =0, and ij,,=i;=0. The above conditions 

are similar to those used in the theoretical analysis 
by Cutkosky of four-terminal-pair networks [5]. 


with the conditions ig=i;', ep 


=0, and with the conditions 


7 
ig =p, 


3. Auxiliary Balances 


The techniques used to assure the realization of 
the conditions defining the inductors involve the use 
of several auxiliary generators and detectors. Figure 
4 shows the complete system. Each of the four auxil- 
iary generators consists of a decade capacitance net- 
work in parallel with a decade conductance network, 
driven by a 40-turn winding on the magnetic core of 
the bridge supply transformer. Figure 4 shows two 
such 40-turn windings, the sign of the voltage of one 
being opposite to that of the other so that each capaci- 
tance decade or conductance decade can be switched 
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FIGURE 4. Complete system for measurement of small mutual inductances. 


independently to reverse the sign of the voltage driving 
it. Stray capacitance from the 40-turn winding to the 
shield raises the potential of the shield and would 
result in current in the shield. This difficulty is over- 
come if an adjustable capacitor, Cg, is attached to the 
20-turn winding, the voltage of which is opposite in 
sign to that of the 40-turn winding. A capacitance 
bridge is formed by placing a detector between the 
shield and the bottom of the 40-turn winding, and the 
bridge can be balanced by adjusting Cg. A similar 
bridge is formed with C/, as one component. Perma- 
nent adjustments of C; and C{, were made following 
the completion of this component of the system. 


The load on the secondary winding of each inductor 
is supplied by an auxiliary generator. The loading 
consists of capacitance in the cable connecting the 


.inductors to the inductive voltage divider, T, and 


capacitance in the inductive voltage divider itself. 
Each of the auxiliary generators is adjusted until un- 
plugging the inner conductor of the coaxial connection 
at the inductor secondary shows no change in the null 
at detector No. 1. Such a condition can be expressed 
as i,=0 and i,=0. The impedances of the auxiliary 
generators are high relative to the inductor secondary 
windings so that convergence to balance is more 
quickly achieved. Without the auxiliary generators, 
the loading in the circuit is such that the voltage e, 
at connector A of a 100 wH inductor was changed by 
approximately one part in 10° when the load was dis- 
connected from the circuit by unplugging the lead to 
the secondary winding. Thus the auxiliary generators 
need to be adjusted to one part in 10* to obtain accuracy 
of one part in 108. 
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At the 1 mH level, an error arises from a capacitance 
change at the secondary terminal when a coaxial cable 
is connected to the terminal. When sucha connection 
is made, the inner connector of the coaxial cable slips 
over the inner member of the connector on the inductor, 
thus changing the capacitance from the inner connector 
of the inductor to its shield about 0.3 pF. The im- 
pedance of the secondary winding of the 1 mH inductor 
used in this work is about 20. A calculation shows 
that one could expect an error of about 6 in 108 due to 
such a capacitance change. This error can be made 
reproducible by using the same coaxial cable each 
time the 1 mH inductor is measured. 

The inductive voltage divider T is a two-stage device. 
Two stage devices have been analyzed by Gibbings 
[4] and by Cuikosky [6]. The auxiliary generator 
(gen. 4 in fig. 4) incorporated in the transformer sup- 
plies the excitation current in the ratio winding. With 
the bridge balanced, a detector winding on the second 
core assures that no excitation current is supplied 
by the inductor secondary windings. Thus this induc- 
tive voltage divider behaves like a device having infinite 
input impedance in that it draws no current. An 
additional feature of the two-stage design is a highly 
accurate 10 to 1 ratio [4). 

The conditions eg =0 and eg =0 are met by adjust- 
ment of the auxiliary generator 1] in the part of the 
circuit containing the primary windings of the mutual 
inductors such that detector No. 3 indicates null. 
Shunt admittances in conjunction with the impedance 
of the coaxial cable from the detector junction to points 
B or E may cause anerror. A calculation based on an 
assumed equivalent circuit and approximate values 
for the impedance and admittance indicates that such 
an error is in the order of 1 in 108 or less for typical 
length of coaxial cable used at this point. An experi- 
mental verification of the above calculation was made 
by observing the difference in the main bridge balance 
caused by changing the length of coaxial cable between 
B and E. 

Tests of interdependence of the detector balances 
were made by bringing the bridge to balance, and then 
changing one auxiliary generator voltage while observ- 
ing changes in the null conditions at the detector junc- 
tions. Such tests and experience in bringing all the 
detectors to a simultaneous null indicate some depend- 
ence of one balance on another. However, the time 
involved in obtaining convergence of all balances is 
not unreasonably long. 
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FIGURE 5. Method of assuring zero net current at inductor secondary 


terminal. 


The condition i,=i'/ can be assured by an auxiliary 
device illustrated in figure 5. One can adjust e,4 (to 
null Ds) to compensate for net current which can result 
from capacitance from the primary to the secondary 
winding of an inductor. However, it is possible to 
construct inductors in such a way that the change in 
the main bridge balance (detector 1) when e, is switched 
in and out is small enough so that the corresponding 
error is less than 1 part in 108. The 1 mH mutual 
inductor was altered by inserting a permanent short 
between the low-voltage side of the primary winding 
to the low-voltage side of the secondary. Inductors 
smaller than 1 mH were not fixed with a short as de- 
scribed above since the error corresponding to the 
change in the main bridge balance when e, is switched 
in and out in their case was smaller than 1 in 10° with- 
out the short. 


4. Main Bridge Balance 


The circuit used to obtain the small adjustable 
voltage, es; is shown in figure 6. That portion of figure 
6 shown within dotted lines represents a component 
part of an “‘a-c direct-reading-ratio set’’ which has 
been described by Cutkosky [6]. This is an active 
circuit with high input impedance and low output 
impedance. Hence Mz is not strongly loaded, and 
reasonably large loads can be driven by the circuit. 
The in-phase and quadrature components of es; are 
read directly from dials of two inductive voltage divid- 
ers, IVD#1 and #2 in figure 6. The mutual inductor 
Mz has nearly the same current in its primary winding 
as the inductor M,;. The voltage eg is related, through 
current i, to e; and e2. The relationship is controlled 
by the choice of the value of M3, which provides a 
scale factor. The value of M; has been calculated 
such that one step of the last dial is one part in 108. 
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FIGURE 6. Basic bridge and, in dotted lines, circuit used to obtain 
balance. 
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The range of es is +5 in 104 and the least count is | in 
108. The current in M3, therefore, must be no more 
than 2 in 10° different from that in M, for accuracy 
of 1 in 108. A calculation shows that ground capaci- 
tance between M, and Mz and in My; must not exceed 
1000 pF for this accuracy. The inductor M3 is con- 
structed to have ground capacitance well below 1000 
pF. Components of the apparatus shown within 
dotted lines in figure 6 were initially adjusted to an 
accuracy of a few parts in 10°. Past experience with 
the type components used in the apparatus indicates 
that at the time of these measurements the accuracy 
of the dial readings was better than 1 part in 10° of 
the dial reading. Thus if the inductors being com- 
pared (i.e., M; and M2) are adjusted so that their values 
differ by less than 1 in 10°, the measurements can be 
relied upon to 1 in 108. 


5. Measurement Procedure 


A single detector is used for all balances; it is 
plugged alternately into positions D;, De, and Ds (fig. 
4). Balancing starts by bringing D, to null by the dials 
IVD#1 and IVD#2 (fig. 6). Next, leaving the detec- 
tor at D,, the secondaries of the mutual inductors M, 
and M, are unplugged in turn. With the secondary 
M, unplugged, a null at D; is obtained by turning the 
dials of gen. 3 (fig. 4); similarly, the secondary of M2 
is unplugged, and the balance reached with gen. 2. 
The detector is then plugged into the D2 position and 
a null is achieved by turning dials of gen. 4. Finally, 
the detector is plugged into the Ds position and brought 
to null by gen. 1. 

All balances described above are then repeated 
several times until convergence is reached. 


6. Inductors 


Some of the requirements of the inductors are 
that they be stable, astatic, have small phase angle, 
be close to nominal value, and be capable of carrying 
substantial current in their primary windings. 

To make the inductors astatic, both the primary and 
secondary have toroidal symmetry about the same 
axis; the primary winding is outside the secondary 
winding and, therefore, has a larger cross section. 
The leads, as shown in figure 7, are brought out to the 
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FIGURE 7. 











Lead geometry and shielding of mutual inductor. 


coaxial connectors in such a way as to create a mini- 
mum loop. The inductor is mounted in a copper 
shield of /s-in thickness (about two skin depths at 
1592 Hz). The coaxial connector of one of the pri- 
mary leads is connected to the copper shield at one 
point inside the shield. The primary leads and the 
secondary leads are brought out on opposite sides of 
the toroid as shown in figure 7 to keep magnetic coup- 
ling to a minimum outside the toroid. 

If two windings have toroidal symmetry about a 
common axis, and if one of these windings is totally 
within the other, the mutual inductor thus formed has 
been observed to have the following property: The 
phase angle of the inductor is practically independent 
of the thickness of the wire used for the outside wind- 
ing; the phase angle depends on the thickness of the 
wire used on the inside winding, a larger phase angle 
being produced by thicker wire. The outside winding 
is, therefore, better suited for the primary, since the 
primary should have low resistance for minimum heat 
dissipation. 

A few inductors have been constructed using toroidal 
forms of Bakelite; they were potted in epoxy resin 
and enclosed in ¥/s-in thick copper shield cans. These 
inductors are useful for testing some of the features 
of the inductance bridge, but have temperature coefh- 
cients that are too large to make them useful as stand- 
ards in this work. 


A 0.1 wH inductor was made by painting silver 
windings onto a toroidal form of fused silica. The 


silica form, with the silver secondary winding, was 
baked and then a primary of copper wire was wound 
around the form but mechanically isolated from it. 
A critical problem is that of bringing the leads from 
the secondary out through the primary (see fig. 8). 
A cylindrical fused-silica rod was fused to the toroidal 
form as in figure 9 so that the leads could be brought 
out on a solid mount. An attempt to use silver painted 
leads was made but abandoned in favor of using #36 
copper wire. These lead wires were tightly twisted 
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FIGURE 8. Winding geometry of mutual inductors having fused 
silica toroidal form. 
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FIGURE 9. Fused silica toroidal form and lead mount. 


and attached to the fused-silica rod with alkyd resin 
paint, and electrical contact to the secondary was 
made with silver paint. The gage of the wire was 
chosen to be as small as possible to keep the loop 
area small, but large enough such that lead resistance 
was below the resistance of the painted secondary. 

A possible source of uncertainty due to relative 
motion between primary and secondary is the magnetic 
coupling from the “single turn” produced by the pro- 
gression of the turns about the toroid. This is an 
effective loop around the toroid in which an emf can 
be induced by a component of magnetic field parallel 
to the axis of the toroid. To reduce this uncertainty 
all the skew of the turns was concentrated at the inner 
surface of the toroid as shown in (fig. 10). After a thin 
insulating coating of alkyd resin paint was applied, 
a return loop was painted on top of the skewed wires 
forming a loop in the opposite direction to the loop 
formed by the skew. 

A nonuniformity of the secondary winding painted 
on the toroidal form could create a net loop which 
would cause the inductors in effect to be lopsided. 
Thus relative motion between the primary and sec- 
ondary windings could change the inductance. A 
calculation for the 0.1 wH inductor indicates that such 
a net loop of 1 mm? would dictate that for stability 
of one part in 108, the relative motion of primary to 
secondary must be 10-5 mm. 

It appears that improved techniques in painting 
the secondary winding onto a fused-silica toroidal 
form would result in inductors of good stability. 


(Paper 70C4—233) 





FicurRE 10. Illustration of painted secondary winding showing 


skew of turns concentrated in inner toroidal surface. 


7. Conclusion 


At the smallest level of inductance measured here, 
bridge sensitivity is one part in 10’. Sensitivity im- 
proves at larger values of inductance such that 100 wH 
can be compared with 1 mH with sensitivity of a few 
parts in 108. It appears likely that the bridge tech- 
niques described here can be relied on to yield errors 
that are less than one part in 107. An analysis of all 
errors in such a step-up from 0.1 wH to 1 mH must 
await the development of a more reliable set of mutual 
inductance standards. 


The author thanks R. D. Cutkosky for suggesting 
the scheme of measurement used here. 
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Bassali’s general theory for the flexure of the thin circular elastic plate supported at an arbitrary 
number of points and subjected to transverse load over an eccentric circle is specialized to the case 
of a centrally loaded plate supported at points equally spaced on a circle concentric with the center. 
Simplified methods for approximating the results predicted by the more complicated theoretical ex- 
pressions for deflection are presented along with the experimental results from 138 tests. Both the 
experimental results and the simplified equations are compared with the theory and agreement is 
found to be good. 


Key Words: Bassali’s theory, concentric loading, circular plates, deflection, elasticity, experi- 
mental, flexure, point supports, simplified approximate solutions, thin plates. 


1. Introduction 


The determination of the deflection of a centrally loaded circular plate supported at points 
equally spaced on a circle concentric with the center has long been an important structural analysis 
problem. In the past the analysis of this problem was usually limited by the assumption that the 
point supports were numerous enough to constitute a simple continuous line support. Nadai [1]? 
presented a theory for the deformation of a circular plate supported at several points with central 
point load or uniform load which was an improvement in that it recognized the errors involved in 
the aforementioned assumption. Unfortunately, Nadai’s point supports were located along the 
circumference of the plate. To some extent this limited the utility of the theory, as this method 
of Support is unusually difficult to realize in practical structures. 

More recently, Bassali [2] has given the solution of the problem of flexure of a thin circular 
elastic plate supported at an arbitrary number of points which may be located anywhere within 
the plate periphery, and loaded over a circular area lying anywhere within the boundary of the 
plate. Implicit in the work of Bassali is the solution of the problem of the centrally loaded plate 
supported at points equally spaced on a circle concentric with a central load. It may be noted 
that the theory accounts for the constraining effect of an annular region of the plate which overhangs 
the support circle and is otherwise free from restraint. 

This paper deals with the specialization of that part of the Bassali theory necessary to solve 
the particular problem described above, and presents the rather tedious theoretical expressions 
for the deflection at the center of the plate and at a point midway between supports located along 
the support circle. Since these expressions require considerable effort to evaluate, simplified 
methods of approximating the center deflection may be desirable for design purposes. Therefore 
simplified expressions for center deflection, based on the results of the exact theory, are given. 
Experimental results are then compared with those obtained analytically, and good agreement 
between the two is found for the range of geometries tested. 


' Figures in brackets indicate the literature references ‘at the end of this paper. 





1.1. List of Symbols 








x,y rectangular Cartesian coordinates 

r = Vx?+ 7’, radial polar coordinate 

0 arctan (y/x), polar coordinate angle 

z = re'®=x+ iy, complex coordinate of point (r, 0) or (x. y) 

z =x—iy, conjugate complex coordinate 

c radius of the plate 

b radius of transversely loaded area of plate 

a span radius, of circle on which support points lie 

h thickness of the plate 

m number of support points (m = 3) 

a =27/m, polar angle between support points 

As =sa, polar angle subtended by the sth support point (s=1,2, . . ., m) 

ds = 0-9, 

2s = ae's 

z inverse of zs with respect to C : z;Z,=c? 

be =z—2s 

Ze =z/—z 

pi = por"~*, transverse load intensity over 0<r<b 
(n= 2, produces uniformly distributed load) 

p2 = 0, transverse load intensity over b<r<c 

Po = 27pob"/n, total load on the plate 

P, = P»/m, reaction at 2s 

b’ = bV n/(n+ 2) 

p =rlc 

t =alc 

q = b/c 

q’ =b'/c 

i =c/a 

Rs(r, 6) = |Z,|, distance from z to zx 

Rir, 0) = |Z;j), distance from z to z; 

c =z/c = pe’® 

a =2z,/c = tes 

ts = te~*s, complex conjugate of C, 

Li: = z./c = te'*e 

v Poisson’s ratio 

K =(3+ v)/(v—1) 

E modulus of elasticity 

D = Eh?/12(1 — v?), the flexural rigidity 

W, deflection due to load, eq (1.01) 

W, deflection due to point load P, 

w(z) deflection of the plate within |z| < 6 

w2(z) deflection of the plate within 6 < lz|<c 

We =w);(0), deflection at the center of the plate 

w(t, 0) deflection at any point on support circle (r= a) 

Wa deflection at (r= a) and g;=(1—2s)7/m, s=1, 2... . , m where b<a 

Woc = w(0) 

C the edge of the plate: the circle r=c 
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Co the boundary of the loaded area of the plate: the circle r=b 


Be dimensionless deflection function defined by eq (4.02) 

B;. dimensionless deflection function defined by eqs (4.16) and (4.18) 
Ba dimensionless deflection function defined by eq (4.04) 

Ww dimensionless parameter defined by eq (4.17). 


2. Survey of Bassali’s Procedure 


No extensive recapitulation of Bassali’s paper [2] need be given here. Indeed, the theory is 
of such complexity that the mere presentation of the expressions for deflection is itself a lengthy 
process. It suffices to summarize the method used by Bassali and to note that his results were 
obtained through the use of the well-known theory of Muskhelishvili [3]. 

The center of a circular plate of radius c is taken to be at the origin of a Cartesian coordinate 
system, and the central plane of the plate to lie in the xy plane. The plate is assumed to carry a 
transverse load over the region? 0<r<b.3 The intensity of this load is 


p=pi=por"~? (n > 2) (1.01) 


(n=2 produces uniform load). The transverse load intensity is p=p2=0 over the remainder of 
the plate. The total load on the plate is Po=2apob"/n. The plate is taken to be supported at 
points equally spaced along the circle r=a (b <a Sc), viz, at z;=ae** (s=1. 2... . .m). From 
conditions of static equilibrium and symmetry the reactive force is the same at each support point 
zs, and is P,;=Po/m. 

To obtain the deflection it is necessary to find the particular solutions for the plate equation, 

for eq (1.01), and for the point loads Po/m. The former is easily seen by inspection to be 

9 — nt2 

W == (20k) * 1.02 
=n +22D oZo ( ) 


It is well known from the theory of plates [4] that for the isolated point load P, at zs 


Ps R2 In Rg. (1.03) 


Ws= aD 


Bassali considers the plate to be comprised of two separate plates. One of these is defined 
by the loaded circular area (0 <r <b), and the other is the annular plate (6 < r <c) with the point 
loads P;. The solution of the equation for the plate deflection with the load p; is then the sum of 
eq (1.02) and the solution of the biharmonic equation; the solution for the annular plate is eq (1.03), 
summed over all of the point supports, plus the solution of the biharmonic equation. Boundary 
conditions for the free edge r=c, together with kinematic and dynamic conditions of continuity 
along r=b6 [5] then make it posible to obtain the deflection, w, for any point on the plate. 


3. Deflection of the Circular Plate 


The deflection of the plate within the loaded region, |z| < 6, obtained from Bassali’s work is 


1) (2-994 2— pn+2) 


87D e n n(n + 2)2b” 


= Po 2 2 _l z 2 Rs ( = ¢ 
Ww; = Wo + lr +b") Inq Ty Re In +{1 (2.01) 


s=1 


2Bassali’s solution is more general, in that, the loaded area and supports may be arbitrarily located anywhere within the plate periphery. 
>For the reader's convenience a list of symbols is given in section 1.1 
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while in the region 6 <|z|<c it is 











7 Po ; ie m J 
w2=Wot oD Jie+o ) In p - 2 R? In -s (2.02) 
where 
_ Po a, Pea m(1 — p?) (a? + Kb") 
Wo = canal b> R In r +(k l)c Re {Ln()} = e+] +yixt yoy t ys 
(2.03) 
and 
Rp, 0)=a? (1 ie cos es) (2.04) 
R2 
2 Pept COS gs + t? (2.05) 
ae =1]1—2pt cos ¢.+ p’t?. (2.06) 


The quantity L»(¢) which appears in eq (2.03) plays an important role in the deflection, and is 
defined by 


m 1 _ -“ m = = = _ 
LunfQ=S) | (1— Seu) In (1 — f.fu) a S) [- olSso) + Seb +1 — S50) In 1 —%50)] (2.07) 
s=1 40 s=1 
and where g(Ls0) is defined as the dilogarithm [6, 7, 8] 


= RE a d = 
at=— | In (1—7) a \f.f|<1. (2.08) 


This function possesses the absolutely and uniformly convergent expansion [8] 





et= > SP fed<. (2.09) 


n=1 


Since GC = pteivs, eq (2.09) can be written as 


—_ nif oo (pt)” * 
ees) = > 7m ° 8 (2.10) 
and if €=%,, k=1, 2, . . ., m, eq (2.10) becomes 


= l 
Dd Poss == gle). (2.11) 


s=1 


From eqs (2.07) and (2.10) the real part of L,,(@) can be written 


Ms 


Re {L»(f)} =— s — Ss cos ng. +> (1—pt cos ¢s) In (1—2pt cos ¢s+ p?t?) 
n=1 s=1 


$s 


1 


—— oe pt sin gs 
pt 2 sin gs arctan 7 aan (2.12) 
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When (= &, = te'**, eq (2.12) becomes 


Re {Ln(Zx)} =-+ elem) +5 Ss (1—¢#? cos sa) In (1— 2t? cos sa+ t*) 


s=1 


m t? sin sa 
—? i se ase” sme (2.13) 
t 2 sin sa arctan 7 cos sa 
For the purposes of this investigation it suffices to consider the case in which the support 
points lie outside the loaded area, i.e., g<t. Since the central plane of the undeformed plate 
is taken to lie in the xy plane, y:=y2=0 in eq (2.03). Further, if the deflection is taken to vanish 
at a support point z,, then 


UZ) = WAZK) = O. (2.14) 


This, along with eqs (2.02) and (2.03) serve to determine ys: 


= —_ 2 2 x e 2 Rs 
¥3= 3D k(a? + b’2) In 1+< > R2(a, Ox) In : (a, Ox) 


m P Re (a, Ox) ] : : {f= 
2, pod ame oi P 5 ne 
2 Rifa, 0) In a a) mt De Re thuld)} 





t?)(a2+ “|. 


ak (2.15) 


1 
m 
It is a notational convenience to define A(t) and By(t) as 


p | m . Rs ie Rs (a, et 
wAnlt) =, Rian On) | In (as Be) I Oy 


and 


mt?B,,(t) = Re {Lm(Cx)}. 


Using eqs (2.04), (2.05), (2.06), and (2.13) the definitions of Am(t) and Bn(t) may be written as follows: 


MA nt) = b 3 (1—cos sa)[« In 2#2(1 — cos sa) —In (1— 22? cos sat t*)] (2.16) 


s=1 


: oe m ; ; ; 
mt*Bm(t)=5 2 (1—#? cos sa) In (1—2#? cos sat t4) 








=! 2m)... 42 ~ ; _# sin sa __ sin Sa ° € 
= g(t?")—t > sin s@ arctan ae peer (2.17) 
Using these definitions, eq (2.15) can be rewritten as 4 
— 2 2 
¥3= Po jar{An(e) +0 — x?)B,,(t)} — x(a? + 6”) In alt (a + xb | (2.18) 
877K a | 


The deflection of the plate is given by eqs (2.01) and (2.02), where wo of eq (2.03) is obtained by 
inserting y:=y2=0 and the value of ys; given by eq (2.18). 

The deflections at the center of the plate and along the support circle are of particular interest. 
At the center r=p=0, so that Rs=a, R;/c=t, and R{/r, =1, hence eq (2.01) yields 


4 Bassali’s equation for ys is in error in that he has + «(a?— b”) In ¢ where eq (2.18) has — x(a? +6") Int. His expression for Am(¢) is also in error, since the first 
#? as well as the s in both logarithmic terms of eq (2.16) are absent in Bassali’s work. 
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ee G In g—a? In r—(1 ~-| b? reesid 


87D n) — n(nt+2/2 (2.19) 


where 


] — t?) (a2+ «Kb’’) 


9 9 9 12 ( 
+ a®2{Ani(t) + (1 — k?)Bn(t)} — K(a? + 6”) In t— oe 





Woe = 





Py [a?+Kb” 

(2.20) 
87KD K+] 
This latter results from noting that eq (2.07) gives L,,(0)=0 and substituting this, and ys; from eq 
(2.18) in eq (2.03). The replacement of woe in eq (2.19) by the right-hand member of eq (2.20) 
gives 


We = Po |at{ Ant) + — x?)B,,(t)} — «(2a>+ b'*) Int 
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Pa?+Kb") _,,,(n+3) 


+ xb” In q+ asi K (n+ 2) 


|: (2.21) 


On the support circle p=t (or r=a), so that eqs (2.04) and (2.06) yield 


Rt, 0)=2a(1—cos ¢s) (2.22) 
_.n l—2F cos 9, +t (2.23) 
* 


s 


from which it follows that on the support circle eqs (2.02) and (2.03) become 





w(t, 0)=wolt, ee {mx (1 +f) Int—k 2 (1—cos gs) In [2021 — cos en)| 


and 


_ fe ~ ; m(1 — t?) kq” (x2 — 1) nn 
Wot, aero |S a—cos ¢s) In (1 — 20? cos gs + t4) += > — (1+ 2 ioe ram Re {L,,( te’ } 


s=1 


The value of ys; in the latter is 





2 42 — 12 
se [— me (1427) in -— 20) 2 (1+%¢ 


13 Br«Dm K+] t ) + m{Anlt) +1 — e0B m0} 


which is a slightly different form of eq (2.18). These expressions combine to yield 


Poa? m (x2 — 1) 
Jo) — _ P = 2 ~ = 4 
w(t, 0) als > (1—cos ¢s) In (1— 22? cos gs, + t*) + 


s=1 = 








Re {L,,(te®)} 


m 


—-«K = (1—cos ¢s) In [20211 — cos ¢s)] +m[An(t)+(1 — K)Bwit)]}- (2.24) 
s=1 


It is of interest to note that q’ does not appear in eq (2.24), so that the deflection along the support 
circle is independent of the radius of the central loaded area as long as b < a. 
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3.1. The Concentrated Central Load 


In this particularly important case the relevant form of wis uw» of eq (2.02), and b=b'=q=q'=0 
must be inserted in this and in eq (2.03) to produce the expression for the deflection. The final 


result is 





_ yf, i242... & 
w=wots—p (- In p 2 p> R? |n ) (3.01) 
where 
— | Spey Rs 2 ])-2 m(1 5 | 2 1 
wo= sD 2 R? In r! +(k?—1)c? Re {Ln(Q)} + 7 +3 (3.02) 
and 
_ Poa — oa _vi—#) : : 
¥3 = BoD [Ante +( K?)B,(t)—k« Int a | (3.03) 


For comparison with experiment it is necessary to have the deflection at the center of the plate. 
In this instance p=9, so that eqs (2.04), (2.06), and (2.12) yield 


Ul 


R,=a, x= 1, and Re {L(0)}=0. 


8s 


Combining eqs (3.01), (3.02), and (3.03) for this case gives 


Poa? | 
87KD 





Fe i 

We = w(0) = An(t) + (1 — «2)Bm(t) + —~ — 2x In t }- (3.04) 
K+] 

The configuration for which the supports are at the edge of the centrally loaded plate is also 


of interest. The expression for center deflection is obtained by setting t=1 in eq (3.04). The 
result is 


Poa? ] 
lw 1a 1 m S Sear 3.05 
Wet 8aKD 4 (1) +( K?)B (1) al (3.05) 
where 
mA »(1)= 4(k — 1) p> sin? = In (2 sin =) (3.06) 
: MAn(1) 1 & 5 Tr 
B» (1) =——— -= — s ————— = 
mB (1) %x—-1) 2 > (7 — sa) sin sa a (3.07) 


The last term in eq (3.07) comes from the fact that g(l)=— [6]. 
The case represented by eq (3.05) was considered by Nadai. In order to compare results 
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with Nadai’s work Bassali calculated w,; for m=3 (using v= 1/4) and obtained 


P 2 
Wer = 0.754 (3.08) 





which coincides with the value obtained by Nadai.® 


4. Simplified Approximation of the Deflection Formulas 


The deflection equations presented in the preceding section are of such complexity that it 
is a lengthy process to write the expressions for deflection at any point on the plate. Therefore, 
a more simple and direct way of presenting the theoretical equations is devised. To illustrate 
this, attention is directed to the expressions for deflections at two points: the center of the plate 
and the point on the support circle midway between adjacent supports. 

Inspection of eqs (3.01), (3.02), and (3.03) shows that by the substitution of D= Eh3/12(1 — v?) 
the deflection for concentrated central loading can be written in the form w= Poa?/Eh® times a 


dimensionless function of p, 6, t, v, and m. In particular, for the center deflection eq (3.04) 
gives 


Poa? 
We = Be Eh (4.01) 
where the center deflection function is 
_301—v’) Es a | 
Be —— [Ane + K )B y(t) + ae 2k Int (4.02) 


Similarly, the deflection at the point on the support circle midway between consecutive 
support points may be written 


Poa? 
wa= Ba pos (4.03) 
where the deflection function is 
_ X1—»*) _ (x2—1) ” 
me "SZ |An(e) + K°)Bm(t) + —— Re {Lite in )} 


+4 S {1—cos = (1-25) In {1-20 cos = (1—25)+| 
m m m 


s=1 


--3 {1 —cos = (1 —2s)} In 202 {} ~cos = (1 —23)}| (4.04) 


s=1 


This is obtained by inserting 6=277/m in eq (2.24), i.e., the angular coordinate of the point on the 
support circle midway between 6=0 and 6=2z7/m. 


5]t should be noted that the erroneous term obtained by Bassali in ys; will have no effect on eq (3.05) since both that term and the correct one given in eq (2.18) 
vanish when t=1. Even though it is beyond the scope of interest of this paper, it is noted here that an error also exists in Bassali’s expression for y3 for the case 
q2t. Asinthe case g <1 above, the erroneous term vanishes whent=1. Further, there is evidence that the error noted in Bassali’s expression for Ap(t) is merely 
typographical, since he gives the correct result in eqs (3.06) and (3.07). 
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THEORETICAL DEFLECTION FUNCTIONS 
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FicurE 1. Relationship of theoretical deflection func- Ve 
tions to number of supports. : 
FIGURE 2. Approximation of B. when m—> ~. 


In figure 1 B, and Ba given by eqs (4.02) and (4.04), respectively, are plotted against m, with t 
as a parameter. The value of v was taken to be 0.3. Examination of this figure shows that the 
deflection functions change rapidly with m for m<6. However, very little change in the deflec- 
tion functions for m > 6 was indicated by eqs (4.02) and (4.04). Thus the error involved in assum- 
ing (in accord with the practice mentioned earlier) that multiple point supports constitute contin- 
uous line support changes rapidly as m becomes smaller than 6, but changes very little for m > 6. 
As illustrated in figure 1 Ba approaches zero as m increases. Further, it is expected that when m 
becomes sufficiently large the support circle will tend to act as a continuous line support. 

Unfortunately the deflection function, B-, can not be evaluated directly from eq (4.02) by 
inserting m= to represent the case of the simple continuous line support. However, the solu- 
tion of eq (4.02) for this case can be approximated in the following manner. Theoretical values 
of B- were computed from eq (4.02) for m=3, 4, 5, and 6 with ¢t and v as parameters.. Figure 2 
presents typical results in the form of a family of curves for v=0.3 with t as a parameter. For 
each t these curves are assumed to be of the form 


Be = ao + aim-* + agm-§ (4.05) 


where do, ai, and a2 are constants. The method of least squares was used to evaluate the con- 
stants. To illustrate, when t=1 and v=0.3 eq (4.05) becomes 


Be =0.551 + 4.172m-3 + 17.493m-6. (4.06) 


Equation (4.06) is a mathematical expression of the t=1 curve shown in figure 2. Thus at m=, 
B-=0.551, and from eq (4.01) the center deflection becomes 
. ca! Poa? 
We = 0.551 Eh? (4.07) 
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Comparing these results with those of the well-known expression for deflection [4] 


_ (3+ v)Pa? 
Wmax ~ T6741 + v)D (4.08) 


and supplying D = Eh3/12(1 — v?) where v=0.3, the center deflection becomes 


Poa? 


Wmax = 0.551 Fp (4.09) 


The agreement of eqs (4.07) and (4.09) lends credence to the assumptions used to develop the 
approximate solution illustrated by eq (4.05). 

Since it is evident that the evaluation of the theoretical deflection functions presented earlier 
requires considerable effort, it follows that the equations may have limited practical usefulness. 
It was for this reason that a simplified method of approximating the more important deflection 
function, B,, (center deflection) was devised. 

The development of the simplified approximate solution was accomplished through a trial 
process of plotting and comparing the theoretical values of 8, against functions of the variables 
t, m, and v in eq (4.02) and applying the least squares method of fitting curves to the theoretical 
results. Various forms of equations were tried, including the exponential and power forms. 
However, the most satisfactory was the consecutive applications of linear least squares in the 
following sequence: for each value of m and v an equation was written in the form 


Be = ao + bot® (4.10) 


where symbols a and 6 with subscripts do not refer to the previously defined radii, a, and 6; then 
equations were written for ao and bo for each value of m in the forms 


ago ay + 6b,(1 —p") 


(4.11) 
bo = az + b2(1 — v?) 
finally, equations were written for the symbols in eq (4.11) in the forms 
a; = a3 + b3m-3 
b; =ay+ bsm-3 
(4.12) 
a2=a;+ bsm-3 
b2 = ag+ bgm~3. 
Combining this group of linear equations and substituting them into eq (4.10) gave 
Bc = a3 + bgm-* + (ay + bgm-’) (1 — v?) + [as + b3m-3 + (ag + bem-*) (1 — v?) ]t3 (4.13) 


or 
B-=— 0.0642 — 2.1900m-* + (0.5687 + 3.2542m—*) (1 — v?) 
+ [— 0.3793 + 11.0513m -* + (0.5223 — 7.8535m —)(1 — v?) |e. (4.14) 
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Since the simplified method of computing B, given in eq (4.14) was developed through the 
curve fitting process described above, this equation reflects the arbitrarily selected limitations of 
the parameters. These limitations were as follows: 


0.25 < vp <0.33 0.50 <t <= 1.00 m= 3, 4, 5, 6, and ©. 


The largest disparity between the results from eqs (4.14) and (4.02) ® computed within the above 
stated limitations was 1.26 percent. 

It is interesting to compare the expressions for center deflection when the load is applied 
uniformly (n= 2) to the plate through a centrally located circular area described by radius b and 
for center deflection when 6 is permitted to shrink to zero. Making use of the definition of a de- 
flection function, B, eq (2.21) can be written as 











ee Poa? - 
We = B. Eh (4.15) 
where 
3(1 — v?) Kq’ - Kq? Kq 5 Kq’ | 
=———— | A(t) + (1 — «7B, (t) +—> In + + —— —— 2x In t~—~ Int 
= ere on eee ae ae; ae 


(4.16) 


w- and B- for the concentrated load at the center of the plate are given by eqs (4.01) and (4.02), 
respectively. By taking the difference between eqs (4.16) and (4.02) all of the terms in eq (4.16) 
containing g can be collected into a separate expression and written as 


_31-Y)¢/ 2 4-3] = 
v= ra (ate! n (4.17) 





~ 


when b<a. This maneuver makes it possible to express B/. as 
B.=Bet+ (4.18) 


but what is more important, it presents an opportunity to extend the utility of the simplified ap- 
proximate solution for B, given in eq (4.14) to approximate the value of 8; within the stated limi- 
tations of eq (4.14). 


5. Test Specimens 


A total of eight thin plate specimens were made for use in the experimental portion of this 
investigation. All of these flat circular specimens were made from bare sheets of 7075-T6 alumi- 
num alloy. The modulus of elasticity for each sheet of aluminum was determined from tests of 
four tensile coupons. In all cases the direction of roll of the sheets was noted and two tensile 
specimens were oriented parallel and two perpendicular to the direction of roll. There were only 
small random differences between the moduli of elasticity with respect to the direction of roll. 
Therefore the average of the moduli parallel and perpendicular to the direction of roll of the sheet 
was selected to represent the moduli of elasticity of the plate specimens fabricated from the 
sheet. The physical properties of the plates are given in table 1. For computational purposes 
the Poisson’s ratio for all of the specimens was assumed to be 0.3. 


® Computations of deflection for m= were based on the approximation illustrated by eq (4.05). 
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TABLE 1. Physical properties of plate specimens 








Modulus 
Specimen | Thickness Radius of 
h c 


elasticity 
E 














in in psi 
A... 0.125 5.125 10.0 x 10° 
Bess 125 4.625 10.0 x 10° 
ee 125 4.375 10.0 x 10° 
= 161 5.125 10.8 x 106 
acces -161 4.625 10.8 x 10° 
WP issied .161 4.375 10.8 x 108 
= seer 250 5.125 9.2 x 10° 
H... | 132 5.125 10.0 x 10° 





The direction of roll was also noted on the plate specimens, and made to correspond to the 
0° to 180° orientation of the plates. The dimensions of the specimens are given in table 1. Fur- 
ther examination of table 1 will reveal that Specimens A, B, and C were fabricated from the same 
0.125 in. thick sheet, and Specimens D, E, and F were fabricated from a sheet 0.161] in. thick. 


Specimens G and H were fabricated from separate sheets of 0.250 in. and 0.132 in. thicknesses, 
respectively. 


6. Test Methods 


The specimens were supported by 1/4 in. diam threaded studs having a machine finished 
spherical shape on the end that bore against the specimens. These supports were equipped with 
lock nuts and inserted into tapped holes in the steel base block shown in figure 3. This base was 
18°/s in. square and 5%/s in. thick with a 3°/4 in. diam center hole. Positioning the support studs 
in this base made it possible to orient the supports to close tolerances for the cases of 3, 4, and 5 
equally spaced supports at the span radii listed in table 2. 

Deflection measurements were made with dial gages having a least division of 0.0001 in. 
These gages were mounted on the heavy steel base block described above using it as a datum plane. 

Two. different methods of applying load to the specimens were used. The diagrammetric 
sketches shown in figure 3 depict these methods. The dead weight loading technique was found 
to be the most convenient to use. It was necessary to drill and tap a 4-20 hole in the center of 
seven of the eight specimens to facilitate load application. A small eye bolt was screwed into 
this hole and locked into position to receive the loading pan. It was felt that this small disconti- 
nuity would not seriously effect the load-deflection characteristics of the specimens, and the 


7 TESTING MACHINE FRAME 


K 





























DIAL GAGES — — SPECIMEN 1 ; 

STEEL ™ = ’ 

BASE oe =] — SPECIMEN 
BLOCK SUPPORT = S 
FOR BASE . = 

BLOOK \ = =| —STEEL BASE 

LOADING = = 

PAN | BLOCK 

20 Ib WEIGHTS 
\— DIAL GAGES 
(a) 


FIGURE 3. Methods used to apply load to test specimens. 
3a. Load applied by dead weight. 
3b. Load applied by testing machine. 
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advantage of conveniently applying known load increments would outweigh the disadvantage of 
loading through the screw threads. To verify this assumption Specimen H was made as a com- 
panion to Specimens A, D, and G except that it did not have the small tapped center hole and 
was loaded in a testing machine frame through a 600-pound capacity elastic load measuring device 
as shown in figure 3. 

The test procedure in both methods of load application was the same. The load was applied 
in 20-pound increments with the center deflection and deflection midway between supports at 
r=a being recorded for each increment of load. 


TABLE 2. Span radii of equally spaced point supports 





Number of supports 
Specimen 





3 4 5 





in in in 


AGD.......... 5.00 5,00 


4.25 4.25 
3.50 3.50 


3.00 3.00 


B&eE.......... 4.50 4.50 
4.25 4.25 
3.63 3.63 


2.75 2.75 


Rien anim cs ans 5.00 5.00 5.00 
4.50 4.50 
4.25 4.25 
3.50 3.50 3.50 














7. Experimental Results 


The experimental deflection data of the plate tests are given in table 3 as the deflection func- 
tions, B- and Ba, defined by eqs (4.02) and (4.04), respectively. These experimental deflection 
functions were computed from eqs (4.01) and (4.03) using the deflection sensitivity, w/Po, as deter- 
mined from the test data. Typical load-deflection data are shown in figure 4. The reciprocals 
of the slopes of these curves are the deflection sensitivities in inches per pound. As mentioned 
previously the plates were tested on 3, 4, and 5 equally spaced point supports with the range of 
span radii as shown in table 2. Table 3 is arranged to show the support conditions as well as the 
experimental results for all of the plate specimens. 

Comparisons of the experimental and theoretical deflection functions are given in figure 5. 
This figure presents three coordinate systems representing the cases m=3, 4, and 5. The deflec- 
tion functions are shown as the ordinate and ¢, is the abscissa. The theoretical relationships 
between £6 and t from eqs (4.02) and (4.04) for a Poisson’s ratio of »=0.3 are represented by the 
solid lines. As indicated in the figure there was good agreement between the experimental and 
theoretical results. 
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TABLE 3. 


Experimental deflection functions, B..and By 





Number of supports 





















































Ratio of - 
Speci- | Span |span to speci- ] 
men radii men radii 3 
a -—— —— - 
B- Ba B. Ba * B. 
ia 
A 5.00 0.9756 0.6763 semecente 0.5447 0.1104 
5.00 9756 6771 0.3007 5486 1126 
5.00 9756 -6834 3382 5423 1087 
5.00 | eee ee OES 5496 .1079 
4.50 8780 .6468 ee 
4.50 8780 6546 2371 ieceeepieianasuiis oe 
4.25 8292 EE Bvesscscsucece 4935 .0649 
4.25 8292 5891 .1967 4908 0658 
4.25 eg SE REP Ba an wet -4986 .0693 
3.50 6829 5308 -1328 75 .0427 
3.50 6829 5308 -1366 4739 0419 
3.25 6341 5325 of 
3.25 6341 5103 A209 1 .... dene ane 
3.00 5854 5013 -1085 4516 .0430 
3.00 5854 4947 -1085 4531 0439 
2.31 4512 4765 -1095 
2.31 4512 4761 .1095 
B 9730 7163 .2940 
ES) AREA TE Rree OR Senn een 
9189 6691 7 eee & 
. _ £ BSS Fees 
7 3 ees 
-7838 -1722 
-7568 -1631 |. 
7568 .1634 |. 
5946 1247 
5946 -1139 
5000 -1095 
-5000 -1062 
( 4.25 9714 6627 2735 
3.50 8000 5547 1788 
3.50 .8000 5531 1796 |. 
3.00 6857 5197 -1440 |. 
3.00 6857 5186 -1451 
2.63 6000 5113 1267 
2.63 6000 5101 -1278 
2.31 5286 4958 1267 
2.31 5286 4867 7367 = 
D 5.00 0.9756 | 0.7649"| 0.3245 0.1262 
5.00 8756 . i= .1287 
5.00 975 - : w 1287 
4.50 8780 .7079 .2506 sa 
4.25 8292 6472 .2202 0774 
4.25 .8292 6455 2195 0779 
3.50 6829 5888 1586 0534 
3.50 -6829 5873 1579 0552 
3.25 6341 5654 1532 
3.25 -6341 5842 -1553 " finden 
3.00 5854 5574 1322 5193 0496 
3.00 5854 5534 1347 258 0551 
2.31 4512 5293 1114 
2.31 4512 5304 1161 
E 4.50 9730 .7304 -3339 | 0.5957 
4.50 9730 7419 3388 } .6037 
4.25 9189 -7021 : _ f Se Sees 5783 
4.25 9189 -7009 .2741 5756 
3.63 7838 6294 2018 | 5546 
| 3.63 -7838 6328 -2004 5443 
| 3.50 -7568 6157 .1840 5307 
3.50 7568 6090 (| as eee 5399 
2.75 5946 5650 -1401 |... 5120 
2.75 5946 5668 1377 |. 5191 
2.31 -5000 5582 -1308 
2.31 -5000 5491 -1313 
F 4.25 9714 7156 3252 
4.25 .9714 714 3252 
3.50 8000 -6083 -2031 
3.50 8000 6204 -2046 |.. 
3.00 -6857 3774 1573 
3.00 6857 O71 -1573 
2.63 6000 5 -1372 |. 
2.63 6000 5579 -1431 
2.31 5286 -5280 -1446 
2.31 5286 5268 1490 }.. 
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0.0530 
0543 
0408 
0411 
0294 


0.0657 
.0661 
0507 
.0494 
0343 
.0343 
0357 
.0346 
.0387 
.0393 


TABLE 3. 


Experimental deflection functions, B,. and By—Continued 





Speci- 
men 


LOAD, P, , Ibf 

































































Number of supports 
P >tio of iScccnes — ee es 
Span _ |span to speci- - 
radii men radii 3 4 5 
a t 
B. | Ba Be Ba Be Ba 
in | 
9756 -6688 y 
9756 6677 
8780 6113 
8780 5993 
8292 5949 
8292 5830 
6829 5497 
6829 5509 
6341 5521 
6341 5426 
9756 .6675 
.9756 -6644 
8292 5800 
8292 5764 
.6829 5478 
6829 5343 
3.00 5854 5201 
3.00 5854 5199 
ij T 
10OF- o 
r DEFLECTION MIDWAY BETWEEN 
SUPPORTS AT r= 3.50in 
80- “7 
6O0F CENTER DEFLECTION — 
40 = 
SPECIMEN A 
a2 3.50 in 
20 t= 0.6829 = 
ms: 4 
re) M 1 1 1 | 1 1 1 n l 1 1 1 1 
(@) 0.01 0.02 0.03 
DEFLECTION, W, in 
FIGURE 4. Typical load versus deflection test data. 
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FIGURE 5. Comparison of theoretical and experimental deflection functions. 


Since the maximum deflection, represented by B, is usually of more interest, a direct com- 
parison of the experimental and theoretical 8, values is made in figure 6 with an indication of the 
order of magnitude of error involved. 

All of the experimental data represented by circular symbols in figure 6 is presented in table 3, 
and represents the full range of this experimental investigation for m=3. The solid circular 
symbols labeled A and H were identified to show the effect of a tapped center hole in Specimen A 
on the load-deflection characteristics as compared with Specimen H which did not have a center 
hole and was approximately the same thickness. It can be seen from the data that the effect was 
so small that the disparity of experimental points from all of the other specimens was greater 


than that for companion Specimens A and H. The triangular symbols in figure 6 represent results 
reported by Nadai [1] on three very thin glass plates. 
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THEORETICAL B, 


8. Discussion 


Bassali obtained the solution for the problem of flexure of a thin circular elastic plate supported 
at several arbitrarily located points (the boundary of the plate was free), and transversely loaded 
over an eccentrically located circular area (the load was symmetrically distributed with respect 
to the center of the loaded area). The general solution provided by Bassali was specialized for 
this investigation of the problem of the concentrically loaded. plate supported at points equally 
spaced on a circle concentric with the central load, and located in the unloaded area, i.e., q < t. 

As noted previously, two expressions for deflection are required; w; being the deflection of the 
plate within the loaded region (r < 6), and we the deflection in the unloaded region (r= 6b). The 
deflection of the plate obtained from Bassali’s work is given in eqs (2.01) and (2.02), where wo is 
obtained from eq (2.03) by inserting yi = y2=0 and the value of y3 given by eq (2.18). It was noted 
that certain errors were apparent in Bassali’s work (see footnotes 4 and 5), but these were of a 
minor nature and stand corrected in eqs (2.16) and (2.18) of this paper for the case q St. 

Bassali’s work was specialized further to examine the case of the concentrated central load 
for the purpose of comparing his theory with the experimental results from 138 tests presented 
herein. It was found convenient to compare the measured deflections at the center and at a point 
midway between supports on the support circle with the theory in terms of the dimensionless 
deflection functions, B- and Ba. given in eqs (4.02) and (4.04), respectively. This comparison was 
made in figures 5 and 6, and good agreement between theoretical and experimental results was 
apparent. 

It is of some importance to note from figure 5 that both the theoretical and experimental 
results indicate that the deflection functions decrease as the annular portion of the plate over- 
hanging the support circle increases. It is of further interest to note from figures 1 and 2 that 
this constraining effect decreases as the number of supports, m, increases, but does not vanish 
for B- when m becomes sufficiently large to produce a support condition equivalent to that of a 
simple continuous line support. From this it would appear desirable to consider the effect of the 
overhanging portion of the plate in the design and analysis of practical structures of this type for 
both conditions of support discussed above. 
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Bassali’s expressions for deflection are of such complexity that they are probably unsuitable 
for design purposes, and there is no direct solution available from his work for the case of simple 
continuous line support. However, the family of curves shown in figure 2 along with the apparent 
y— intercepts suggested the simple method of approximating B, for the simple support condition 
as illustrated by eq (4.05). The result of this approximation given in eq (4.07), where t= 1, agreed 
favorably with that from the classical elementary theory for the case of simple support at the 
boundary of the plate. Thus, the assumptions used to develop the approximate solution given in 
eqs (4.05) and (4.06) were substantiated. 

Further interest in simplified design equations produced the expression for approximating 
the theoretical center deflection function, B-: The development of this expression, which approxi- 
mates the results of the theoretical expression in eq (4.02), is described by eqs (4.10—4.13), and 
presented in eq (4.14). As stated previously the largest difference between the results of the ap- 
proximate solution and the theoretical solution computed within the limitations inherent in eq 
(4.14) was 1.26 percent. Since such good agreement exists between the approximate and theoreti- 
cal results, eq (4.14) should prove useful for design purposes. It should be noted here that setting 
m~*=0O in eq (4.14) produces the approximate expression for 8, when there is a condition of 
simple continuous line support at t. This case was discussed earlier and illustrated in eq (4.05). 

A comparison of expressions for the theoretical center deflection functions, B, and B’-, given 
in eqs (4.02) and (4.16), respectively, indicated that the approximate solution for B, given in eq 
(4.14) could be modified to approximate B’, of eq (4.16) as indicated in eqs (4.17) and (4.18). Since 
Ww is obtained ‘directly from the theory, no further disparity is expected between the theoretical and 
the approximate £’. values. 


9. Conclusions 


1. In view of the good agreement between the theoretical and experimental results from the 
138 tests reported herein, Bassali’s theory adequately predicts the deflection of the plates over the 
range of geometries investigated. 


2. Both the theoretical and experimental results show that the deflection of the plate within 
the support circle is constrained by the annular portion of the plate which overhangs the support 
circle. This constraint is found to decrease as the overhang is decreased, and to increase as the 
number of supports are decreased. Since this constraining effect is not negligible for a plate 
overhanging a support circle containing many supports, it is an important consideration in the 


design and analysis of thin plate structures, and especially so when the number of supports are 
fewer than six. 


3. The simplified equations for deflection functions presented herein may be used for design 
purposes within the stated limitations of the equations. 
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NBS has needed a set of very reproducible germanium resistors that would be capable of main- 


taining temperature scales. 


This paper describes our procedure for selecting such a set. 


A group of germanium resistors from three commercial sources have been thermally cycled be- 
tween 4.2 °K and room temperature. The resistance-temperature calibrations at 4.2 °K were made 
with reference to a liquid helium bath so that reproducible temperatures could be determined from 


liquid helium vapor pressure measurements. 


Seven resistors out of 25 demonstrated reproducibilities 
(of the 4.2 °K calibration) of about 0.001 °K after undergoing the multiple cyclings. 


Guided by these 


results, we procured twelve, new, similar resistors which were cycled in a comparable procedure. 
After 85 cycles, in which 14 resistance-temperature calibrations were performed at 4.2 °K for each 
resistor, 10 of the 12 resistors demonstrated reproducibilities of approximately 0.001 °K. 

The resulting set of secondary thermometers have undergone calibrations in the temperature 
ranges, 2 to 5 °K and 2 to 20 °K; reports of this work will be published in the near future. 


Key Words:, Germanium resistors, germanium resistor reproducibility, low temperature thermom- 
etry, reproducibility of germanium resistors, thermal cycling at low temperature, 


thermometry. 


1. Introduction 


A major problem in low temperature work has been 
the lack of a highly reproducible, precisely calibrated 
secondary thermometer in the region of 1 to 20 °K. 
Carbon composition resistors found wide popularity 
after Clement et al. [1]2 had reported their repro- 
ducibilities and calibrations in the temperature regions 
of liquid helium and liquid hydrogen. However, 
carbon composition resistors show limited reproduci- 
bilities after cycling from room temperatures to 
4.2 °K [2,3], and exhibit an aging effect when main- 
tained in a highly stabilized liquid helium bath for 
several weeks [4]. 

Encouraging results for resistance thermometry in 
the liquid helium and liquid hydrogen temperature 
ranges were reported by Estermann [5] for doped 
germanium and silicon; by Friedberg [6] for poly- 
crystalline p-type germanium; and by the Bell Tele- 
phone Laboratories group [7, 8] using arsenic-doped 
germanium that was encapsulated in platinum cans. 
A very limited number of the B.T.L. specimens were 
distributed by the Calorimetry Conference to several 
laboratories for evaluation. Results from the several 
laboratories [9] indicated, in general, reproducibilities 
of 0.001 °K. Since that time, the properties of im- 
purity doped germanium resistors have been reported 


'M. H. Edlow is presently associated with the United States Patent Office 


2 Figures in brackets indicate the literature references at the end of this paper. 


by several commercial sources [10,11,12]. This 
paper is concerned with the determination of calibra- 
tion reproducibilities of commercial specimens? at or 
near the normal boiling point of liquid helium. The 
objective of the work was to obtain a group of 10 or 12 
similar resistors that exhibited reproducibilities of 
1 mdeg or better. This group of resistors will serve as 
secondary thermometer standards for maintaining 
temperature scales below 20 °K. Subsequent papers 
will cover their precise resistance-temperature 
relationships in the region 2 to 5 °K and 2 to 20 °K. 


2. Experimental Apparatus 


Our determinations of resistance are based upon 
d-c potentiometric measurements and the attendant 
techniques—four leads to a resistor are required for 
the measurement. A dual-six-dial potentiometer is 
used in conjunction with a special reversing switch 
that affords a reversal of current in all circuits —in- 
ternal potentiometer and external resistors (fig. 1)—to 
compensate for undesired thermal e.m.f.’s. The null 
detector, a millimicrovoltmeter, has a useful range of 
source resistance from 1 to 100,000 © with the capa- 
bility of detecting changes as small as one part in 
sixty thousand for measured resistances of 3,000 1. 
The germanium resistors are connected in series with 
1,000 and 10,000 2 standard resistors [13] to the exter- 


3 At the time that the resistors were procured these were, to the best of our knowledge, 
the only products of American manufacture available. 
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FIGURE 1. D-c electrical circuit for measuring germanium resistors. 
R, are standard resistors; R,, R,, and R,, the resistors to be calibrated. 
nal reversing switch; currents through this circuit are The potentiometer, external reversing switch, 
determined from potentials measured across the stand- batteries, null detector and standard resistors are 


ard resistors. Measuring currents of the order 1 to 
2 mA, are obtained by using two 50,000 0 voltage 
dividers connected in parallel across the battery. A 
double-throw, double-pole toggle switch connects the 
current supply to either the germanium resistor circuit 
or a 1 MO “bleeder” resistor; the bleeder resistor 
allows the battery to deliver a small, continuous cur- 
rent when not in use with the germanium resistors. 
A two-volt, low discharge battery (series connected to 
a resistance decade box) supplies current to the middle 
internal potentiometer circuit while six-volt, low dis- 
charge batteries supply current for both the external 
series circuit containing the germanium resistors and 
the upper internal potentiometer circuit. The bat- 
teries are surrounded by thermal insulation to minimize 
temperature changes of the electrolyte that would 
produce variations in the discharge voltage. Current 
standardizations need only be performed a few times 
per day, although more frequent standardization 
checks were routinely performed. 


supported on acrylic plates, which in turn rest on 
copper sheets that are interconnected by copper wires. 
A lead from one of the sheets is connected to a sur- 
rounding screen cage which is at earth potential. The 
double-wall, copper-screen cage completely encloses 
the d-c measuring apparatus and eliminates a-c joule 
heating that is caused by signals which emanate from 
nearby commercial transmitters [14]. 

The germanium resistors are immersed in a con- 
stant-temperature liquid helium bath contained in 
a 25-liter metal storage Dewar [15]. A vacuum 
tight Dewar cap has been designed to fit over the 
Dewar neck J (fig. 2) and permits routine thermal 
cycling of the resistors. A, B, C, and D are “quick 
couplings” and are hard soldered to the top E; E and 
a threaded brass cylinder G (internal component of 
quick coupling H) are soft soldered into the copper 
body F. The vapor from the helium bath passes 
through a side arm M that is connected to a 2.8 cm 
vapor-pressure pump-line. 
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FIGURE 2. Dewar cap. 


A stainless steel tube that passes through coupling C 
leads to a jacketed vapor pressure thermometer [16] 
that is positioned 0.3 cm from the bottom of the Dewar. 
A surface vapor pressure tube passes through coup- 
ling D and extends from about 10 cm above the coup- 
ling to about 8 cm below the neck of the Dewar. 
Both vapor pressure tubes consist of stainless steel, 
thin-wall (0.015 cm) tubing of 0.65 cm diam. A rubber 
line, when necessary, connects: an external helium 
cylinder to the surface tube and permits helium gas 
to be bled into the Dewar. The gas flows out of 
coupling A during resistor thermal cycling and pre- 
vents frozen air from accumulating on the inside 
surface of the Dewar neck. The coupling B con- 
tains a 0.25 cm quick coupling B’ that is closed by a 
metal rod: it will contain a He® vapor pressure ther- 
mometer in future work. 

The component of the cap that is displaced dur- 
ing thermal cycling of the resistors is shown in figure 3. 
A hollow copper cylinder R, 1.9 cm o.d., slides through 
the coupling A. A multi-tube Kovar-glass seal at 
S is set in a small recess machined at the end of the 
cylinder R. Six pairs of #38 A.W.G. Formvar-coated 
copper wire, T, are strung through the capillary 
tubes; and the region of the seal, tubes and leads 
is covered with wax to achieve a vacuum seal. The 
copper leads are permanently maintained in a com- 
pact bundle within the Dewar by applying a cement 
to the taut wires. The resistors (positioned at the 
end of the leads) touch the Dewar bottom. Spun 


















|| 


FIGURE 3. Dewar cap component that is moved during thermal 


cycling procedure. 


glass spaghetti insulation covers the connections 
between resistor leads and permanent connecting 
wires. Thread is tied around the spaghetti-covered 
leads to achieve a compact bundle which easily passes 
through the 1.9 cm opening of A during thermal 
cycling. Copper wire strung through holes X of the 
arms W and wound tightly around the Dewar handles 
secures the cylinder R tothe cap. The cap is secured 
to the Dewar handles in a similar manner. These pre- 
cautions are necessary in the work from 4.2 to 5 °K 
where the helium vapor pressure is greater than 
atmospheric, and exerts sufficient force to dislodge 
the Dewar cap. A complete picture of the storage 
Dewar and cap is shown in figure 4. 

Groups of three resistors are solder-connected to 
the lower extremities of the wires. (Experience with 
impurity-doped germanium resistors has shown that 
it is necessary to minimize the heat conducted to the 
specimen while soldering resistor leads to connecting 
wires. We have found that a reproducible calibra- 
tion at 4.2 °K has been changed by several milli- 
degrees when solder such as 50/50 lead-tin was used; 
consequently, Wood’s metal is used to solder the 
resistor leads to the connecting wires. A pencil iron 
is used to supply the minimum heat required to melt 
the solder.) Each resistor has four leads—two for 
current and two for potential—that extend from 
within the Dewar, out through the seal at S (fig. 3) 
to ambient conditions. The current leads are series 
connected external to the Dewar so that in event of a 
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FIGURE 4. Complete picture of storage Dewar and cap. 


resistor failure the faulty resistor can be quickly dis- 
connected from the current circuit and measure- 
ments continued on the other two resistors. The po- 
tential leads emerging from the Kovar-glass seal and 
the 1,000 © standard resistor potential leads are 
appropriately soldered to terminals of a rotary se- 
lector switch. The rotary switch is mounted on the 
lid of a glass jar so that the switch is completely 
surrounded by a glass wall to provide: a degree of 
thermal insulation from possible spurious tempera- 
ture fluctuations. The common terminals of the 
switch are connected to one pair of emf terminals of 
the potentiometer by #38 A.W.G. double-nylon 
covered, Formvar-coated copper wire. The poten- 
tial leads from the 10,000 © standard resistor are 
directly connected to the second pair of emf termi- 
nals on the potentiometer. 

The manometer for measuring the helium vapor 
pressure in this work is a digital readout 152.4 cm 
mercury manometer [4]. Either the vapor pressure 
thermometer-bulb or the bath surface vapor pres- 
sure can be connected to the manometer by means 
of insulated 0.6 cm copper tubing (A, A’, and G) and 
the vacuum bellows-type valves C and D (fig. 5). 
Valve F affords isolation of the manometer while the 
vacuum pumping line, through valve K, enables one 


to evacuate any of the components that lead to the 
manifold E. 

The gas filling system consists of a 1-liter storage 
can, L, and needle valve J through which helium gas 
can be metered into the system via rubber tubing H 
and vacuum valve I. Needle valve J’ permits refill- 
ing the can from a nearby helium cylinder. The mani- 
fold of figure 5, the mercury manometer, manifold 
evacuating apparatus, and helium bath vacuum 
“pumping” source are all external to the screened 
room. Tubing ingress to the room is afforded by 
portholes that have been built into a wall of the 
screened room. Our experience has shown that 
effective shielding by the screened room is destroyed 
when an undergrounded conductor is brought through 
the portholes; consequently, the Dewar pumping 
line has been appropriately grounded to the screened 
room. The space in porthole B (fig. 5), containing 
the two vapor pressure lines, is packed tightly with 
brass filament sponge and no evidence of R. F. sig- 
nals is detected. The shielding of the screened 
room is tested at the beginning of every week using 
a battery operated F.M. receiver. 

The well-known phenomenon of an oscillating he- 
lium gas column [18] was occasionally detected in 
the present setup. When this condition occurred in 
connection with the surface vapor-pressure tube, 
the liquid helium evaporation rate, as detected by 
a flow-meter, increased by a factor of 3 and oscilla- 
tions of the manometer mercury miniscus were 
readily observed. A sheet of fine copper screening 
rolled into a tight bundle and inserted in the open 
end of the surface vapor-pressure tube for a length 
of 7 to 10 cm eliminated the troublesome oscillations 
in this instance. 


3. Method of Cycling 


The following procedure has been found satisfac- 
tory (damage to the resistor leads and freezing of the 
leads or resistor specimens to the Dewar neck have 
been avoided) for repetitive thermal cyclings of 
resistors from room temperature to 4.2 °K. Rubber 
tubing connected a cylinder of helium gas to the 
surface vapor pressure tube and a slight overpres- 
sure of gas was accumulated in the Dewar. The 
nut N of coupling A (fig. 3) was loosened and the 
ensemble of cylinder, coupling nut N and portions of 
the connecting wires were withdrawn from the Dewar 
cap; the slight gas overpressure caused a stream of 
helium gas to continually pass through the opening 
of coupling A and thereby minimized any flow of air 
into the Dewar. After the resistors were raised to 
the vicinity of the Dewar cap, several layers of tissue 
paper (or raw cotton) were wrapped around the leads 
and positioned in the opening of the coupling A. 
The helium cylinder was then closed off. The re- 
sistors remained for three minutes within the Dewar 
cap, which is at approximately room temperature, 
and then were brought out into the room for approx- 
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imately 1 min.* (The resistors were not brought 
immediately out into the room in order to minimize 
the accumulation of condensed room moisture on 
the encapsulations. A sufficient quantity of mois- 
ture can lead to mechanical stressing of the resistor 
leads.) 

Before the resistors were cycled back into the 
liquid bath, the helium cylinder valve was again 
opened to produce a steady stream of gas-flow through 
coupling A. As soon as the resistors approached 
the surface of the bath, the helium cylinder valve was 
closed since the then-increased evaporation rate of 
the bath was sufficient to maintain a slight overpres- 
sure. The entire process of lowering the resistors 
into the helium bath consumed only a minute and 
inflicted a thermal shock on the resistors. 

The resistances of the specimens were measured 
soon after they were immersed in the liquid helium; 
at least 15 seconds were required to make measure- 
ments that would yield a value of resistance, which 
indicated an approximate condition of thermal equi- 
librium between the resistors and the surrounding 
liquid helium bath. However, the resistors were 
allowed to remain in the bath for about one minute 
prior to another cycle (see footnote 4). 


3.1. Determination of a Calibration 
Temperature 


When a well defined calibration point is desired, 
a meaningful thermal equilibrium must be achieved 


‘Our experience has always indicated that the germanium resistors, when warming 
to ambient conditions after a 4.2 °K exposure, require a lapse of only 2 or 3 min before 
attaining nominal ambient temperature. The reverse operation, cooling to a nominal 
1.2 °K from ambient temperatures, is accomplished 10 or 15 sec after immersion of the 
resistor in liquid helium. Private communications from Mr. Herder (Cryocal) and Mr. 
Halverson (Radiation Research) reveal, for comparable cooling and heating treatment, 
very similar experiences. 


K 
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Vapor pressure manifold. 


within the liquid helium bath and between the bath 
and resistors. The vapor pressure of the bath 
is placed under control of the pressure regulator 
[17] and values of vapor pressure (from the vapor 
pressure thermometer) and resistance are moni- 
tored until a steady state condition results. This 
requires about 45 min and on many occasions we have 
ascertained that the thermal equilibrium attained in 
this period of time will not change more than a few 
tenths of a millidegree over periods of days or weeks. 

In practice we do not experience difficulty in associ- 
ating a temperature (near 4.2 °K) with a resistance 
value for a reproducible germanium resistor. The 
surface vapor pressure of the liquid helium bath is 
easily controlled, using our pressure regulator [17], 
within 0.05 mm of mercury—this corresponds to 
0.00007 °K at 4.2 °K—so the reproduction of a bath 
surface vapor pressure does not constitute a prob- 
lem. We have previously reported [16] the existence 
of a pressure gradient in the vicinity of the liquid 
helium surface for the type of bath we employ: in 
the bulk of the liquid helium the temperature is 
relatively constant. Since the surface gradient is 
dependent upon the Dewar evaporation rate, and per- 
haps the Dewar geometry also, it is essential that 
a_ helium vapor pressure bulb and the resistors be 
located within the liquid helium bulk. (From the vapor 
pressure determination and the “1958 He? Scale of 
Temperature” [19], a temperature can be ascertained.) 
We have also found that the vapor pressure thermom- 
eter tube must be insulated so that it does not make 
thermal contact with the “colder” liquid helium 
surface through which the tube must pass. Vacuum- 
jacketing of the tube along the length that would make 
contact with the surface provides proper thermal 
insulation and avoids the effect of a “‘cold spot” 
[16, 20]. Radiant energy is prevented from reach- 
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ing the vapor pressure bulb by a radiation shield 
that is soldered within the tube and extends from 
7.6 cm above the top of the jacket upward for 10.2 
cm. When vapor pressure and temperature have 
been derived from a jacketed vapor pressure ther- 
mometer, and when an immersed germanium re- 
sistor has been stable, we have been able to associate 
resistance measurements with temperature to a 
constancy that is equivalent to 0.0003 °K. However, 
if temperatures are referenced to surface vapor 
pressure measurements, the above constancy can 
decrease to an equivalent of 0.001 °K. The de- 
crease, in part, can be attributed to variations in 
the liquid helium evaporation rate and the helium 
level within the Dewar. 

Within the bulk liquid helium a steady-state tem- 
perature can be routinely reproduced within 0.001 °K. 
For most of the resistors being reported, the well- 
defined calibrations have been conducted at tem- 
peratures that generally were reproduced within 
0.0003 °K (with reference to jacketed vapor pressure 
thermometer measurements). A few of the re- 
sistors were calibrated at reference bath-tempera- 
tures which differed by as much as 0.001 °K. We 
obtained dR/dT at approximately 4.2 °K for the in- 
dividual resistors by performing an additional cali- 
bration of resistance versus temperature (0.050 °K 
or 0.100 °K removed from the reference temperature 
for calibration after thermal cycling); this enabled us, 
through interpolation, to correct all of the nominal 
4.2 °K calibrations of a resistor to a particular tem- 
perature. 


4. Description of the Resistors 


The resistors were obtained from three manufac- 
turers who have periodically supplied us with speci- 
mens for experimentation [21]. Some of the prototypes 
have been described [10, 11, 12] but for purposes of 
adequate comparison we shall briefly describe each. 
Resistors from one source (Minneapolis-Honeywell) are 
encapsulated in a platinum can that contains hermeti- 
cally sealed helium gas. Four thin platinum wires 
from the resistor are brought through a glass header of 
the can and soldered to multistrand, plastic covered, 
copper leads. Mechanical support for the leads is 
achieved by “potting” an epoxy resin around the 
soldered platinum wire-copper lead connections. 
Resistors from a second source (Radiation Research) 
are also encapsulated in a platinum can which contains 
hermetically sealed helium gas. These resistors have 
four 0.025 cm platinum leads that are brought through a 
glass header for electrical measurements. Resistors 
from a third source (Texas Instruments) are in a soft, 
opaque glass envelope and have two 0.025 cm platinum 
leads for electrical measurements. (Envelopes or 
encapsulations may be perforated depending upon the 
users requirements.) Later specimens have been 
encapsulated in a metal can with a transistor header. 
Various manufacturers have used arsenic and gal- 
lium as the major doping impurities; both the re- 


sistivity and its temperature coefficient are effected 
by control of the impurity concentrations. 


4.1. Thermal Properties 


Joule heating, developed when current passes 
through a resistor, raises the temperature of the 
specimen to a value higher than that of its surround- 
ings. Since the usefulness of the germanium re- 
sistance thermometer depends upon its represen- 
tation of the correct temperature of the surrounding 
medium, this effect must be minimized. The results 
for two typical resistors follow. The resistance 
variation with current squared for a reproducible, 
hermetically sealed Minneapolis-Honeywell resistor 
in a stabilized liquid helium bath at 4.2 °K is shown 
in figure 6. The abscissa is the square of the cur- 
rent which is approximately proportional to the 
power being dissipated; at the left, the ordinate is 
the fractional change of resistance caused by joule 
heating, while at the right it is the corresponding in- 
crease in temperature, AT, obtained from the tem- 
perature coefficient of resistance. The increase in 
temperature, A7’, varies linearly with J? up to the 
value /? = 1.3 10-* amperes? which is in reason- 
able agreement with Blakemore’s result [11]. Tested 
with the above resistor was a Minneapolis-Honey- 
well resistor contained in a perforated encapsulation; 
its fractional change of resistance as well as the 
corresponding increase in temperature are also 
plotted as a function of [7 in figure 6. Comparison 
of the two plots shows relatively little or no joule 
heating for the perforated encapsulation, indicat- 
ing that the helium gas (in hermetically sealed speci- 
mens) is a relatively inefficient heat exchanger. A 
similar plot of AR/R, AT versus [? for a Texas In- 
struments encapsulated specimen is shown in figure 
7. The result indicates AT varying linearly through 
most of the /? range but exhibiting a noticeable 
deviation from linearity at /? = 1.6 x 10-°. 
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FIGURE 6. Joule heating for two Minneapolis-Honeywell resistors. 

(Ra2—R)/Ra2a, the fractional change in resistance, and its temperature equivalent, 
AT, are plotted as a function of /?. The nominal 4.2 °K resistance of the hermetically 
sealed and perforated encapsulations are 2800 and 4800 2), respectively. 
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FIGURE 7. Joule heating for a Texas Instrument resistor. 


(Ry2—R)/Ra2, the fractional change in resistance, and its temperature equivalent, 
AT. are plotted as a function of 7*. The nominal 4.2 °K resistance of the hermetically 
encapsulated resistor is 1800 2. 


The lowest points of the curve for the encapsu- 
lated, hermetically sealed resistors (figs. 6 and 7) 
show a negligible change of temperature due to 
joule heating; consequently, a power dissipation of 
about 0.02 ~W was taken as a reasonable maximum 
for joule heating and used as an upper limit for all 
subsequent measurements at 4.2 °K. From the 
linear portion of the curve in figure 6, we can de- 
termine the thermal resistance, Q=dT/d(/?R) (in 
degrees-per-microwatt), between the resistor ele- 
ment and its surroundings. For the Honeywell 
resistor Q=0.0045 °K/uW; Blakemore obtained, 
for a similar resistor, Q=0.0023 °K/uW [10]. From 
figure 7 we find for the Texas Instrument resistor 
Q=0.0017 °K/uW while Low obtained a value Q= 
0.001 °K/pW [12]. 


42. Reproducibility of Resistors 


If a germanium resistor is thermally cycled n times 
from 300 °K to 4.2 °K and systematic calibrations are 
made and referenced to a particular temperature, 
i.e., 4.214 °K, we define the calibration reproduci- 
bility of the resistor to mean: the difference, AR 
(ohms), between the maximum and minimum values 
of resistance that were measured during the cali- 
brations (properly referenced to 4.214 °K) made 
within the n cycles. From knowledge of the re- 
sistor’s dR/dT at 4.2 °K, the AR can be more perti- 
nently expressed in equivalent millidegrees—thus 
we refer to calibration reproducibilities at 4.2 °K 
in terms of millidegrees. 

The results of thermal cycling between 300 °K and 
4.2 °K for 25 resistors are listed in table 1. The 
resistors were obtained from three commercial 
sources and we were seeking information that might 
lead to the selection of a group of similar resistors 
that could be used as secondary standards in our lab- 
oratory. We established the requirement that to 


be useful for our purposes resistors should have a 
4.2 °K calibration reproducibility which was 0.001 °K 
or better. Resistors which failed mechanically 
during the thermal cycling havé not been included 
in table 1. Some of the resistors (numbers 19, 22, 
23, and 25) failed to meet our calibration reproduci- 
bility criterion rather early in the cycling proce- 
dure while others failed to perform satisfactorily 
after a greater number of cycles.’ In general, ac- 
curate 4.2 °K calibrations were performed for the 
first few cycles, and thereafter every 10 cycles. 
Seven of the twenty-five resistors were reproduc- 
ible to 0.0011 °K and have subsequently been used 
rather extensively in our laboratory. For example, 
resistors A and B have been incorporated with the 
acoustical thermometer for several years and will 
be directly referred to in forthcoming publications. 


TABLE 1. Reproducibility of germanium resistors from three com- 
mercial sources under thermal cycling from 4.2 °K 
temperature 


to room 

















Resistor | Number of Average Reproducibility 
serial # cycles R (ohms) 4.2 °K (milli- °K) 

1 61 2829.0 1.3 

2 123 2552.5 0.9 

3 275 2628.0 9 

4 1) 6260.0 - 

A 56 2868.3 1.0 

B 40 3086.0 1.1 

7 87 3945.4 1.3 

8 50 1074.0 1.3 

9° 51 1764.2 1.5 

10° 83 2743.0 2.3 

1] 36 367.75 1.1 

12 36 617.74 1.2 

13 36 292.75 1.5 

14 36 292.50 1.4 

15 64 236.84 6.6 

16 64 327 5.8 

17 40 327 17.4 

18 64 675.1 pC 

19 8 671.60 4.6 
20 30 1890 0.7 
21 30 720 9.0 

22 14 23100 10.0 (estimated) 
23 14 14500 10.0 (estimated) 
24 23 1825 2.1 

25 8 1940 3.8 





@ Resistor encapsulation was perforated. 


An exhaustive, accumulative cycling test was con- 
tinued for resistor number 3. It was thermally cycled 
275 times with equilibrium resistance measurements 
determined in 43 instances. The results are shown 
in figure 8. The first series of tests occurred in 
July 1963 and the resistor was then set aside. On 
January 7, 1964, the resistor was again connected 
and the cycling test continued during January— 
February 1964. The resistance measurements dem- 
onstrate a reproducibility of 0.9 mdeg at 4.2 °K. 

Figure 9 shows the equilibrium resistance-accumu- 
lated cycles plot for resistor number 1. We first 
measured this resistance in July 1963 and, subse- 


5 Resistor number 4 was needed for other experimentation so its thermal cycling was 
not continued. 
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FIGURE 9. Equilibrium resistance as a function of the number of 
accumulated cycles for #1. 
T= 4.2 °K. 


quently, made measurements after every 10 cycles 
until an accumulated total of 53 cycles was reached. 
The resistor leads were disconnected and subse- 
quently resoldered to the connecting wires five times 
over a period of six months. Despite this treatment, 
the resistance shows a reproducibility of 0.0013 
°K. Between the 56th and 57th cycles, the resister 
was inadvertently physically shocked. New con- 
necting wires were soldered to number | and addi- 
tional cycling demonstrated little change from the 
measurements made two and a half months earlier. 


5. Select Resistor Group 


After obtaining a group of resistors that offered 
promising calibration-reproducibilities, we adopted 
the following method as a criterion for determining 
resistance-temperature reproducibilities under thermal 
cycling. First, resistors were to be cycled three to 
five times from room temperature to 4.2 %K prior 
to a first measurement of resistance. This was 
“adopted when the accumulated data over a number 


proximately 1 mdeg (during the previous 10 cycles), 
resistors were to be cycled 75 more times with an 
equilibrium resistance measurement — performed 
after every 25th cycle. A cycle in which a resistance 
measurement is not performed consists of retaining the 
resistor for 3 min at the top of the Dewar neck, 1 mir 
in the room and returning it for one minute’s e» 
posure at the Dewar bottom. By this procedur 
accumulative effects of cycling, if any, were to | 

detected. Thus, a total of some 85 cycles wit 

14 resistance-temperature calibrations were to be 
performed in this “search.” The results, showr. 
in table 2, demonstrate that 10 out of 12 resistors were 
reproducible to approximately 0.001 °K; four t 
0.0007 °K or better. The equilibrium resistances for 
resistors L, G, and I as a function of accumulated 


thermal cycles are shown in figures 10 to. 12. Al 
three resistors were cycled simultaneously anc 
measured within several minutes of each othe 


during the test. Resistors L and I appear to bs 
good thermometers with reproducibilities of 0.000 
°K and 0.0008 °K while G exhibits a reproducibility 
of only 0.0019 °K. Another group of three resistors 
D, J, and N were tested simultaneously and their 
equilibrium resistances versus accumulated cycles are 
displayed in figures 13 to 15. Resistors D, J, and N 
exhibit reproducibilities of 0.45, 0.66, and 0.7 m°K, 
respectively. 
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TABLE 2. Reproducibility of a selected group of resistors under : aaa Saas i ? : : : 
thermal cycling from 4.2 °K to room temperature 
ee _ 27912F Ne 
Resistor | Number of | Average R Reproducibility 2 Ss 
cycles (ohms) 4.2 °K (milli- °K) 5 27908} 475 ¥ 
~ 4 a Ww ar 
3 j 
i 88 2800 1.1 = 27904 16 3 
D 86 2565 0.5 o ra 
E 85 2890 1.0 i 4 
F 88 2678 1. Baad ae 
G 86 2561 1.9 
H 87 2718 “13 27896 1 1 n L 1 n L fo) 
I 86 2790 0.8 0 10 20 30 40 #50 60 70 80 9 
J 86 2593 ye NUMBER OF ACCUMULATED CYCLES 
K 87 2620 1.1 q 
L 86 2672 - 0.7 Rete gih) pr kes se Py eee F 
FIGURE 12. Equilibrium resistance as a function of the number of 
M 87 2659 9 accumulated cycles for resistor I. 
N 86 2769 a T=4.2°K. 











* Resistor not cycled prior to a first measurement. 
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6. Measurement Errors 


The measurement errors for the potential and cur- 
rent readings are estimated to be approximately 
0.005 percent. At 4.2 °K the resistance-temperature 


coefficients (5 o) are approximately 0.06 percent 


per m°K. The error in a resistance determination, 
therefore, is equivalent to 0.17 m°K. The estimated 
error in pressure measurements is 0.0] percent; since 
1dP 
P dT 
0.1 percent per m°K, the estimated error in pressure 
determinations is equivalent to 0.1 m°K. Spurious 
temperature fluctuations within the bulk liquid helium 
bath are estimated as 0.1 or 0.2 m°K. Thus the total 


for helium vapor pressures at 4.2 °K is close to 


measurement error in determining resistance-tempera- | 


ture calibration reproducibilities at 4.2 °K could be- 
come as large as a temperature indeterminacy of 0.4 
or 0.5 m°K. Recognizing that the reproducibilities 
of column 4, table 2, are maximum calibration varia- 
tions (over an average of 14 calibration determinations 
for each resistor), the experimental reproducibilities 
(0.7 m°K or less) for the best resistors are in satisfac- 
tory agreement with the estimated error (0.4 or 


0.5 m°K). 


7. Conclusions 


The germanium resistors that have been reported 
in table 1 were, to our knowledge, not preferentially 
selected. Rather, they constitute representative 
resistors from three commercial sources. Our ex- 
periments have indicated that 14 of the 25 resistors 
reproduce a 4.2 °K resistance-temperature calibration 
within 2.0 m°K while only 7 of the group failed to 
reproduce within 5.0 m°K. This would seem to in- 
dicate a general availability of germanium resistors 
that are suitable for many low temperature purposes 
where the thermometry needs are not exacting. 

The selection of resistors reported in table 2 was 
preferential (guidance was furnished by the results of 
investigations on the first group—table 1) and the 
reproducibilities listed in column 4 clearly indicate 
the general superiority of this selection—7 of the 12 
resistors were reproducible within 1.0 m°K. 

This investigation has fulfilled our original purpose — 
to obtain a group of germanium resistance thermom- 
eters that offered the promise of being suitable for use 
as secondary thermometers. Since the investigation 
was completed nearly 2 years ago, there has been 
ample time for employing the resistors as secondary 
thermometers in subsequent endeavors. Resistors 
A and B (table 1) have been used extensively in 
acoustical thermometry; C and D (table 2) have been 
calibrated with reference to the acoustical thermom- 
eter and are the basic reference for calibrations be- 
tween 2 and 20 °K at the National Bureau of Standards. 
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Articles on analytical representations of germanium 
resistor calibration data are in preparation and will 
appear in the literature. 

Since some of the germanium resistors listed in, 
table 1 were obtained as long ago as 5 or 6 years, we 
have no opinions concerning the current availability 
of particular models. Without doubt some are now 
outdated due to advances in manufacturing technology. 
We have not had the opportunity to investigate groups 
of resistors which have been commercially available 
during the last 3 years. 
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An evaluation of two of the most widely accepted methods for calibrating vibrating-sample mag- 
netometers is given. The comparison method uses a material of known magnetization such as pure 
nickel. In the slope method, the magnetometer is calibrated from the low field linear slope of the 
magnetization curve of a sample of high permeability. 

The primary source of error in the comparison method arises from an uncertainty in the absolute 
magnetization of nickel and its dependence on environmental conditions. The study indicated that 
better accuracy can be expected from the slope method. The use of pure iron in this method was 


found preferable to high permeability ferrites. 
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1. Introduction 


The vibrating-sample magnetometer as developed 
by Foner [1]! is rapidly becoming one of the more 
commonly accepted techniques for determining the 
saturation magnetization of ferrimagnetic materials. 
The many improvements and refinements which have 
been made in these instruments in recent years have 
been thoroughly described by Feldmann and Hunt 
[2, 3, 4]. Although well designed magnetometers of 
this type prese.tly exist, their ultimate accuracy is 
still dependent on the calibration technique used in 
the measurement process. At the present time, two 
different calibration methods have received the 
greatest emphasis and have been included in a 
standard test method for saturation magnetization of 
nonmetallic magnetic materials by the American 
Society for Testing and Materials [5]. In the com- 
parison method the magnetometer is calibrated using 
a material of known magnetization which is usually a 
sphere of pure nickel. In the slope method the instru- 
ment is calibrated from the initial slope of the mag- 
netization curve of a spherical sample of high permea- 
bility. It is the purpose of this paper to describe the 
results of a study of the advantages and disadvantages 
of each of the above two calibration techniques. 

The more commonly used comparison method will 
first be discussed with a brief description of some of 
the precautions which we found necessary for im- 
proving the accuracy of vibrating-sample instruments 
in general. The experimental evaluation of the slope 
method will then be described. The advantages of 
using ultrapure iron rather than high permeability 
ferrites which were mentioned [5] for use in the 
slope method will also be discussed. 








* Radio Standards Physics Division, NBS Boulder Laboratories, Boulder, Colo. 
‘Figures in brackets indicate the literature references at the end of this paper. 


2. Comparison Method 


The saturation magnetization of a magnetic material 
using a Foner type vibrating-sample magnetometer is 
determined from the voltage generated in a set of 
pickup coils by a magnetized spherical sample vibrat- 
ing perpendicular to the applied magnetic field as 
shown in figure 1. The voltage generated in these 
coils is proportional to the magnetization of the sample 
[1]. We may thus write 


Eo= KV (47M), (1) 
where Ep is the voltage induced in the pickup coils by a 
sample of magnetization My and volume Vo, and K is 


the proportionality constant. The factor 47 has been 
inserted since the saturation induction, 47M, given 
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FIGURE 1. Sample and coil arrangement in electromagnet with 
x, y, and z coordinates. 
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in gauss, is normally quoted in the literature in describ- 
ing these materials when using CGS units. In the 
comparison method, K may be readily obtained from 
the measured voltage, E;, obtained using a standard 
sample of volume V, and known magnetization M,. 
Thus from eq (1), the magnetization Mo of an unknown 








sample using the comparison calibration method 
becomes, 
Eo D3 
4aMo=(=*) (58) 4M. 2 
a 2 all “ 


where the ratio of spherical volumes, V,; and Vo, have 
been replaced by the ratio of the corresponding cube 
of the diameters, D, and Do. 

It is apparent that the accuracy of any measurement 
depends directly upon how well one knows the magne- 
tization of the standard sample, M;.. However, before 
confronting this problem it is desirable to reduce the 
possible errors from the voltage and diameter terms 
to a minimum value. 


2.1. Sample Diameter Measurements 


It can be seen that the cube of the diameter of the 
spheres ‘appears in eq (2) which magnifies any errors 
in the diameter terms. It is thus necessary to grind 
highly spherical samples and then measure their 
average diameter to the highest possible accuracy. 

We have found that a sphere grinder such as de- 
scribed by Cross [6] will usually yield spheres that 
are out of round in the order of only 0.2 percent. 
Some spheres have been ground better than 0.05 per- 
cent out of round, while other spheres ground from a 
soft material such as nickel may be 0.5 percent out 
of round. The percentage of out of roundness is 


here defined as 100 —~** = <n Gam 
min diam 

The diameter of the spheres is usually determined 
from an average of 20 to 30 random measurements 
made on a given sphere using an electronic comparator 
calibrated with gage blocks. 

In some cases where the deviations in diameter are 
somewhat larger such as for nickel, it has been found 
more accurate to determine the volume and average 
diameter of the sphere from mass and density measure- 
ments. This procedure still requires obtaining at 
least one good spherical sample in order to accurately 
determine the density. This density is then used in 
determining the volumes of other less round spheres 
of the same material assuming the density is constant. 





2.2. Voltage Measurements 


Several factors will of course influence the voltage 
readings. Repeatability or precision of the voltage 
readings Ey and E, is, of course, of considerable im- 
portance. Since the above voltage readings are from 
two different samples, it is necessary to easily attach 
and remove each sample from the vibrating rod with 
the least disturbance possible. In one approach, the 


sample was glued to a plastic holder which in turn was 
fastened to the vibrating rod with a plastic screw. 
Repeatability varied from 0.2 to 1.5 percent. A better 
approach with a repeatability of less than 0.1 percent 
is to cement the sample to a small plastic holder which 
is in turn held on the end of the vibrating rod by means 
of avacuum. The samples and corresponding holders 
are attached or removed from the rod by merely 
controlling the vacuum. 


Another source of error associated with voltage 
readings is related to the position of the pickup coils 
with respect to the sample position, as has been de- 
scribed by previous authors [1, 3]. To study the 
effects of coil adjustment, the pickup coils in our instru- 
ment were mounted on a plastic beam which was 
fastened to a milling machine table to provide vernier 
adjustment in the three dimensions shown in figure 1. 
Experimental data indicate that the coils must be 
adjusted each time a sample is put in the instrument 
to within a few thousandths of an inch of the maximum 
voltage reading in the y and z direction, and a minimum 
reading in the x direction. For example, a change of 
0.010 inch in the z direction gave a change of 0.2 per- 
cent in output voltage; a 0.020 inch change gave an. 
0.8 percent change in output voltage. The distance be- 
tween the coils from center to center is approximately 
11/s inches. 


Still other precautions were found necessary for 
obtaining accurate voltage readings. (1) The mag- 
netometer was mechanically isolated from the elec- 
tromagnet. Mounting the vibrating rod system on a 
wooden table, independently supported by the floor, 
reduced the residual signal as read on a voltage ratio 
transformer from 0.0014 to 0.0001 at a field of 10000 Oe. 
(2) The sample pickup coils were mechanically 
isolated from the electromagnet. Mounting the pickup 
coils on a floor supported table, independent of the 
magnet, reduced the residual noise. (3) The empty 
holder voltage readings were subtracted from the 
corresponding sample readings. Failure to do this 
can result in a 0.1 percent error in a typical ratio trans- 
former reading of 0.1000 with an empty holder reading 
of 0.0001. (4) A high field was applied to the sample 
before taking measurements coming down the mag- 
netization curve to avoid any hysteresis ambiguities. 
This is of greater importance when measuring at low 
applied fields such as in the slope calibration method 
as will be discussed later in the paper. (5) The tem- 
perature of the magnetometer and sample were sta- 
bilized to a few tenths of a degree centigrade during 
a measurement. This was necessary for holding any 
variations in the calibration constant to 0.1 percent or 
less during a measurement. Any further instability 
effects due to temperature and other sources are 
reduced by calibrating the system both immediately 
before and after a measurement is made on an un- 
known sample and _ interpolating the calibration 
constant as a function of time. 

Another possible source of error in vibrating-sample 
magnetometers is related to magnetic image effects in 
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the pole caps of the electromagnet. Errors from this 
situation are apparently reduced by eddy current 
shielding in the pole caps as pointed out by Foner [1]. 
Using a 2-inch air gap in the magnet, he estimated the 
image effect to be less than 0.2 percent. In comparing 
data on iron samples using a 2-inch and 3-inch air gap 
at 9000 Oe applied field, we obtained agreement of 
better than 0.1 percent which likewise implies that 
image effects are small. 

In addition to the above efforts to avoid errors, a 
further check on the accuracy of the voltage readings 
was made by measuring E for several different sized 
samples of the same material at a given field and 
comparing the E/D® ratios. In order to avoid errors 
due to equipment drift, it was found advantageous to 
refer the E/D® readings of each of the samples to the 
corresponding E/D*® value of one of the samples ar- 
bitrarily chosen as a reference (R). The voltage from 
this reference sample was measured immediately 
before and after each of: the unknown samples and 
the reference voltage used for calculation was obtained 
from interpolation between these two values. Dis- 
agreement of slightly greater than 0.1 percent was 
obtained as shown in figure 2. These results in con- 
junction with the preceding discussion leads to the 
conclusion that it is possible to hold any errors in the 
diameter and voltage terms of eq (2) to less than a 
few tenths of a percent. 


2.3. Magnetization of Standard Sample 


The saturation magnetization of pure nickel is almost 
universally used as a standard when ferromagnetic or 
ferrimagnetic materials are considered. Few investi- 
gators have measured the absolute value of 47M, of 
nickel. The results of a literature search as shown in 
table 1 indicate several percent variation among those 
values which have been reported. Most of these 
values in the literature are given in terms of emu per 
gram. We have not converted these values to 47M, 
in gauss since the density of the reported material 
was not always known. Even if one chooses the cor- 
rect value, it is difficult to know if his own nickel sam- 
ple is of the same composition and is being measured 
under the same conditions that were present when the 
absolute value was determined. This problem results 
from the fact that the magnetization of any given 
sample of nickel depends upon purity, density, mag- 
netic field, temperature, strains, annealing, etc. 

The literature contains equations or graphs for 
correcting for density, temperature, and field effects 
when the experimental conditions differ from those 
related to the absolute values [12]. However, some 
of these corrections may not in themselves be entirely 
correct. A typical example might be the equation for 
the approach to saturation which could be used to 
predict the magnetization of nickel at some other 
field than where the absolute value is quoted. How- 
ever, assuming that accurate corrections could be 
made for temperature, density, and applied field, a 
problem still exists regarding purity and strain effects. 
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FIGURE 2. Comparison of E/D® ratios of pure iron spheres as a 
function of sample diameter. 


The data were normalized to the £/D® ratio of one of the spheres. arbitrarily chosen as 
pure iron reference sample. (R) 


TABLE 1. Values of saturation moment (emu) per gram 
for pure nickel 











Ref. Remarks 
emu/lg 
{7] | Single crystal. *Calculated emu. Reference gives 
—H=6500 gauss. Temperature and field not 
given. 
55.5 +1 [8] | Temp. 25 °C. Field not given. 
55.05 + .05 {9} | Temp. 15 °C. Extrapolated to infinite field. 
54.98 [10] | *Calculated from (B— H)x=6150 gauss for density 
trams 
8.901 em? and infinite held 
54.73° {11} | Temp. 20°C. Field 18000 (Oe.). 
54.55” {11} | Temp. 20°C. Field 8000 (Oe). 
54.39” {12} | Temp. 15 °C. Field not given. 
54.27” [5] | Room temp. Field not given. *Calculated from 
B—H=4rM,=6070 gauss and density=8.90 a 
<<" in ee ee ee ee et 
‘Data in reference given in B—H gauss. From this we calculated emu = B-H 


47 density 
»These values are based on Weiss values. 


In general, one hopes that his particular sample is of 
the same degree of purity as the material used in ob- 
taining the absolute values quoted in the literature. 
The problem of strains can become quite serious as 
can be seen in figure 3 which shows the magnetization 
curves obtained before and after annealing four 
spheres ground from an ultrapure rod of nickel with 
impurities of less than 0.005 percent. The variation 
in the data apparently results from small strains in- 
duced in the material in the grinding process. Such 
strains are probably removed by annealing since the 
curves after annealing are in much better agreement. 

The annealing process itself may be a problem since 
the literature is full of examples of rather large changes 
in the magnetization curves of nickel due to different 
annealing treatments. Fortunately such effects do 
not greatly influence the true saturation value at high 
fields, but they can become a problem if attempts are 
made to calibrate an instrument at lower applied fields. 
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FIGURE 3. Magnetization curves before and after annealing four 


pure nickel spherical samples. 


On the basis of the above arguments, it can be seen 
that the primary uncertainty in calibrating a vibrating- 
sample magnetometer using the comparison method 
may well reside in the value of the standard sample. 
As was shown, the measured diameters of the spheres 
and corresponding voltages in eq (2) can be deter- 
mined to a few tenths of a percent. However, unless 
extreme care is exercised in choosing the value of the 
standard sample, a further error exceeding 1 percent 
may be introduced into the problem. In the following 
section, we will show that the slope method was devel- 
oped primarily to avoid this problem of uncertainty 
in 47M, of the standard sample. 


3. Slope Method 


3.1. Evaluation of Error 


The slope method is based on the observation that 
the voltage in the pickup coils from a spherical sample 
is a linear function of the applied field over the lower 
region of the magnetization curve [13]. Typical curves 
for a high permeability ferrite and a pure iron spl.ere 
showing this situation may be seen in figure 4. As 
shown in the appendix, the slope of these curves in 
the linear region is relatively insensitive to the perme- 
ability of the material. It is also shown that this 
allows the equation for determining the unknown 
magnetization of a test sample at some field H using 
the slope method to be written as follows: 


EoD? 


1 a." 
P= +N| aA D8 





47M, = (3) 























1.00 T T T | T 
VOLTAGE VERSUS EXTERNAL ca 
0.90;- MAGNETIC FIELD FOR A Pad 
SPHERICAL SAMPLE OF ULTRA 
0.80} PURE IRON AND A HIGH 
PERMEABILITY FERRITE “ ; 
0.70} SPHERE DIAMETER. i 
= IRON = 0.07672" Z 
§ FERRITE = 0.07665" 
S 0.60 /- nae 
_ 
= ° 
0.50 7 
WwW 
Oo 
<a 
F 0.40 F | 
re) 
> 
0.30 fA HIGH PERMEABILITY 
? < FERRITE 
0.20 L.—+-—— 
0.10 F SLOPE = ae + 
c 
oe l 
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 
EXTERNAL MAGNETIC FIELD, Ho (Oe) 
. : : , : sample voltage 
FIGURE 4. Comparison of the linear region of the et 


: reference voltage 
ratio curves for an ultrapure iron and a high permeability fer- 
rite sample of equal size. 


where D,.=diameter of calibrating sphere, 


Do = diameter of unknown sample sphere, 


N=demagnetizing factor of calibrating 
sample, 

w=relative permeability of calibrating 
sample, 


AH,=change in de magnetic field applied to 
calibrating sample in linear portion of 
curve (see fig. 4), 
E,=coil voltage from unknown sample minus 
empty holder reading, 
AF,.=change in coil voltage from calibrating 
sample corresponding to AH,. 

For greater precision the two voltage readings used 
for determining AE,, were taken at H, and zero applied 
field. Although figure 4 implies that the voltage is 
zero at the origin, a considerably expanded plot of 
these results shows a small voltage corresponding to 
a slight remanence at zero applied field is present. 

We have already pointed out that an accuracy in the 
order of 0.1 percent to 0.2 percent is not unreasonable 
in the voltage and diameter measurements as was 
shown in figure 2. The three remaining values in eq 
(3) are the permeability, w, the demagnetizing factor N, 
and the field, H. used in determining the slope 
AE./AH,. The problem of comparing the two calibra- 
tion methods thus reduces to a comparison of the ac- 
curacy with which we can determine p, N, and H, to 
the accuracy with which we know 47M, of a standard 
sample as used in the previously described method. 

Any errors in the measurement of H,. may readily 
be reduced to a negligible amount by using a nuclear 
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resonance gaussmeter in conjunction with a frequency 
counter. AE,/AH, should be taken for several values 
of H, to be sure that the linear portion of the curve is 
being used. 


: oo : 1 
Consider now the remaining quantity Pex paul 
u— 


The permeability can, of course, have a wide range of 
values depending upon the type of calibrating sample 
used. However, as Frederick [13] pointed out the 
versatility of this method is considerably enhanced by 
choosing a material of high permeability. Under 





these conditions, the term is very small com- 


pared with N=1/3. For example, if we use ferrites 
whose permeability is at least 2000, then the difference 
in assuming ~=© compared to using the initial per- 
meability of 2000 in eq (3) would be only 0.15 percent. 
Thus for samples with permeabilities of 2000 or greater, 
the error contributed by uncertainties in the permea- 
bility term should be less than this. The fact that we 
do not need to know p accurately if it is sufficiently 
high in this slope method is the primary advantage of 
this technique over the comparison method. Further- 
more, materials having permeabilities greater than 
2000 are readily available and may easily be ground into 
spherical shapes. 

In the case of a perfect sphere, N= 1/3; however, 
any deviations from a true sphere would cause an 
error in N which would cause practically the same 
error in eq (3). The effect of out of roundness can be 
calculated, theoretically, if we assume that the out of 
round sample is still a prolate or oblate spheroid. 
Stoner [14] has shown that for the nearly spherical pro- 


late spheroid with a > b, 2 = -(?) <1 
— oe a 2 3 2 
Np=3-7€ a. . .) (4) 
we ee 
Nq=5 1—-Mp)=3 +75 € (14 ze. ; (5) 


and for the nearly spherical oblate spheroid with 


a<b,e=1-(2) <1 
1m 4, 
Np=3+ 5 (1+56 ee .) (6) 
1 ] 
Ny=5 (1-N))=3-75 @ +fe+. ‘ .) (7) 


where a=the polar semi-axis, 
b=the equatorial semi-axis, 
N,= demagnetizing factor along the polar axis, 
Nq= demagnetizing factor along the equatorial 
axis. 
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FIGURE 5. (a) Theoretical effect of small eccentricity, €, on values 
of N when magnetic field, H,, is applied first along the polar 
axis and then along the equatorial axis for oblate and prolate’ 
spheriods. 

(b) Total percent error in N expected between field along polar axis 
and along equatorial axis. 
(Data are obtained from the difference between corresponding curves in graph (a)). Ex 


perimental points represent percent error in voltage readings which are proportional 
to percent error in N. 


Equations (4), (5), (6), and (7) are shown graphically 
in figure 5a where 
€? Max diam — Min diam 


2 Min diam 





As an experiment, E, was measured along the polar 
axis and compared with E, measured along the equa- 
torial axis for two samples. The error (difference) in 
E. as shown in figure 5b, which is proportional to the 
corresponding error in N, was less than the theoretical 
curve predicted by taking differences in N in figure 5a. 
It was concluded that the calibrating spheres should 
not be out of round more than approximately 0.1 per- 
cent. The errors due to sample volume determination 
are actually greater than those due to misorientation. 
The diameter measurement has already been dis- 
cussed earlier in the paper. 

The above arguments thus imply that the error to 
be expected from the AH., uw, and N terms in eq (3) 
should be less than a few tenths percent. This result 
in conjunction with the 0.1 to 0.2 percent possible 
error in the voltage and diameter terms makes it fea- 
sible to expect errors of less than 0.5 percent in using 
eq (3). 
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samples with initial permeabilities ranging from 1000 to 3000. 


The data were normalized to the arbitrarily chosen pure iron reference sample (R) 
used in figure 2. 


3.2. Experimental Results Using High 
Permeability Ferrites 


In order to further confirm the above analysis, an 
experimental investigation of the validity of the cali- 
bration constant obtained using the slope method was 
carried out. The calibration constant may be readily 
obtained by writing eg (3) in the form 





1 Eo 
rly =—-— Th2? 
where K;, is given by 
(+4) AE, 
: a—1 
K, = D?AH. (9) 


On the basis of the previous discussion, it appears 
reasonable that we may set N=1/3 in this equation 
with negligible effect on the accuracy. With this 
assumption, the value of K; as a calibration constant 
was checked in the following ways: 

(1) For a given calibration field H,, K; was compared 
for several ferrite samples of different size from the 
same high permeability material (u assumed constant). 
For six samples ranging from 0.090 to 0.128-inch diam- 
eter, agreement was within 0.15 percent. No size 
effect trend could be detected. 

(2) For a given calibration field, K,; was compared 
for several samples of nearly the same size but for 
several different ferrite materials with initial permea- 
bilities, uo, that range from 1000 to 5000. As was 


previously discussed, we used an arbitrary reference 
sample before and after each measurement to improve 
the accuracy. Using known values of initial permea- 
bility [uo, in eq (9)], the variation in the K, ratios as 
shown in figure 6 for eight high permeability ferrite 
samples was 0.28 percent. Assuming 4 =, the varia- 
tion was 0.15 percent. This implies that the effective 
permeability is equal to or greater than the initial 
permeability as specified by the manufacturer as shown 
in figure 6. A ninth sample with the highest initial 
permeability (5000) gave 0.5 to 0.85 percent higher 
values for K,. 

(3) For a given sample, Ky was compared for several 
values of the field H,.. In the 350 to 600 Oe range, 
K, varied less than 0.1 percent and there was no defin- 
ite trend as to field dependence. The above procedure 
was repeated for several calibration samples of 
different materials of high permeability with similar 
results. 


3.3. Experimental Results Using Ultrapure Iron 


Because of the discrepancy caused by the highest 
initial permeability (5000) sample and because of the 
above 0.3 percent variation of the other ferrites, further 
investigation was carried out. As can be seen, an ideal 
calibrating material would be one with a high initial 
permeability of 10,000 or greater in which case the 
term. at in eq (9) could be dropped with negligible 


error. Data given in the literature, [15, 16], indicate 
that high purity iron can be produced with initial 
permeabilities in this range. 

A commercially available sample of high purity 
polycrystalline iron was obtained with a total impurity 
of 0.005 percent including less than 0.001 percent 
carbon impurity. Six samples were ground down to 
various sizes and measured with the results shown in 
figure 7. Again iron reference sample (R) was used. 
The less than 0.1 percent variation in K, for iron shown 
in the figure is less than the 0.15 percent variation 
observed between different diameter samples of a 
given ferrite as mentioned in the previous section as 
well as the variation observed between different 
ferrites as shown in figure 6. 

It can be seen from figure 7, that the iron reference 
sample (R) is fairly representative of the iron samples. 
The data in figure 6 were taken using the same iron 
reference sample (R). Therefore, figure 6 is a means 
of comparing how well several ferrites with different 
permeabilities compare with a representative iron 
sample using ~=po and w= for the ferrites and 
w= for iron in eq (9). It appears the correct uw 
value to use for the ferrites lies between po and ~. 
Some additional ultrapure iron samples were ground. 
annealed, and reground. Measurements comparing 
annealed and unannealed ultrapure iron samples 
agreed within 0.1 percent at applied fields of 600 Oc 
and 12.000 Oe. 

The close agreement in the above results indicates 
that high purity iron has several advantages when used 
in the slope method to calibrate a vibrating sample 
magnetometer. These are (1) high purity iron is 
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FIGURE 7. Comparison of Ky ratios for six ultravure iron spheres 
versus sample diameter for an external field of 1125 (Oe). 


The data were normalized to one of the spheres, arbitrarily clfosen as pure iron reference 
sample, (R) 


readily available commercially; (2) it is easier to accu- 
rately calculate K, since w« may be assumed to be 
infinity; (3) the properties of a pure material such as 
iron are more uniform than ferrites as to porosity, 
density, homogeneity, etc., which give a more constant 
K,; and (4) because the slope of the linear portion of 
the iron curve in figure 4 at an applied field of 5000 
Oe deviates only 0.5 percent or less from the slope of 
the curve at 600 Oe, a pure iron sphere could be used 
to calibrate a magnetometer, with some sacrifice in 
accuracy, at any field up to 5000 Oe using eq (9). 

It is apparent that the saturation magnetization of 
a given spherical sample of pure iron at high fields may 
readily be determined using the same spherical sample 
to calibrate the instrument at low fields using the slope 
method. This _ self-calibration procedure avoids 
diameter and sample position errors which may be 
present when the calibration sample and the sample 
under test are different. A measurement of the satura- 
tion magnetization of iron in this manner provides 
another means of checking the accuracy of the mag- 
netometer since this quantity may be compared with 
values given in the literature which are in good agree- 
ment. Using this approach, we found that the satura- 
tion magnetization data taken on five pure iron samples 
agreed to within 0.3 percent of measurements taken 
by Weiss and Forrer [17] and Danan [9]. We meas- 
ured a saturation magnetization value of 21465+20 
gauss for a pure iron sphere at 23.3 °C and an external 


field of 9000 Oe. 


4. Conclusions 


Experimental observations using the slope calibra- 
tion technique show that only a few tenths of a percent 
variation in K; may be expected. Ultrapure iron and 
most high permeability ferrites can be used for 
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calibration. However, ultrapure iron is preferred for 
several reasons as noted above. In an actual measure- 
ment of Mo. the additional error in Do and Ey would also 
be present. These errors as discussed earlier in 
the paper are of the order of one-tenth to a few tenths 
of apercent. These results thus show that an accuracy 
of better than 0.5 percent for magnetization measure- 
ments can easily be obtained with the slope method. 

In the comparison method previously discussed, we 
had the same one-tenth to a few tenths percent error 
in Ey and Dy. In addition, we had the error due to the 
uncertainty in 47M, of the standard nickel sample. 
It is possible that this might be known to a few tenths 
of a percent which would make this method com- 
parable to the slope technique. However, the sen- 
sitivity of nickel to its environment and the uncertainty 
in the true absolute value make this seem unlikely. 
It thus appears that the slope method is preferable to 
the comparison technique for calibrating vibrating- 
sample magnetometers. 


The authors express their thanks to Nolan Frederick 
who designed and built the vibrating-sample mag- 
netometer and to William McNaney for precision 
measurements of the specimen diameters. 


5. Appendix. Derivation of Equation (3) 


The relationship between the magnetization M and 
internal field, Hi, of a material in the demagnetized 
state is defined as [18] 


_w-1l,, 
M="_— i, 


where uw is the relative permeability. The internal 
field, Hi, in a finite specimen is related to the external 
field Hy by means of the demagnetizing factor N as 
follows: [19] 


H;=H )—N4aM. 


Substituting this value of H; in the first equation and 
solving for 477M, we get 


For a calibrating sample, this equation may be written 


on. 
eee Os 
c—) 


As in eq (1) in the text, we have 


E.= KV .An0M- 








or 


ne. —Ky tn, — BUSH 
—-—+N 
a-t 
and thus 
(+) AE. 
1 


K = 





V AH. 


For an unknown sample, Eyg=KVo47Mbp, as in eq (1) 
in the text, 








_ Fo 
47M, = KV, 
Substituting the value of K, 
47M, = EV 4 
l + v| AE. y 
p—-l | AH “ 


or expressing the volumes in terms of diameter, we 
obtain eq (3) in the text. 


ED! 


arr oe 
yt] aa D3 


47M, saa 








It should be noted that the above derivation is based 
on the assumption that the sample is in the demagne- 
tized state. It was previously mentioned that a very 
small voltage due to a slight remanent magnetization 
may be observed at zero applied field in the spherical 
samples. However, the excellent agreement in the 
data obtained for different materials implies that eq 
(3) is still essentially valid even though a slight hys- 
teresis exists in the samples. 
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The “law of propagation of error” is a tool that physical scientists have conveniently and frequently 


used in their work for many years, yet an adequate reference is difficult to find. 


In this paper an exposi- 


tory review of this topic is presented, particularly in the light of current practices and interpretations. 


Examples on the accuracy of the approximations are given. 


results is discussed. 


The reporting of the uncertainties of final 
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Introduction 


In the December 1939, issue of the American 
Physics Teacher, Raymond T. Birge wrote an ex- 
pository paper on “The Propagation of Errors.” 
In the introductory paragraph of his paper, Birge 
remarked: 

“The question of what constitutes the most reliable value to 
be assigned as the uncertainty of any given measured quantity is 
one that has been discussed for many decades and, presumably, 
will continue to be discussed. It is a question that involves many 
considerations and by its very nature has no unique answer. The 
subject of the propagation of errors, on the contrary, is a purely 
mathematical matter, with very definite and easily ascertained 
conclusions. Although the general subject of the present article 
is by no means new,!' many scientists still fail to avail themselves 
of the enlightening conclusions that may often thus be reached, 
while others frequently use the theory incorrectly and thus arrive 
at quite misleading conclusions.” 


Birge’s remark 27 years ago still sounds fitting today. 
For a number of years, the need for an expository 
paper on this topic has been felt by the staff of the 
Statistical Engineering Laboratory at the National 
Bureau of Standards. Frequent inquiries have to 
be answered, yet a diligent search in current litera- 
ture and textbooks failed to produce a suitable ref- 
erence that treats the subject matter adequately. 
The present manuscript was written to fill this need. 


In section 1, we consider the two distinct situations 
under which the propagation of error formulas can 
be used. The mathematical manipulations are the 
same, yet the interpretations of the results are en- 
tirely different. 


fined and the general formulas given. Frequently 


used special formulas are listed at the end of the ° 


section for convenient reference. In section 3 the 
accuracies of the approximations are discussed, 
together with suggestions on the use of the errors 
propagated. Section 4 contains suggestions on the 
reporting of final results. 


' See, for instance, M. Merriman, Method of Least Squares, pp. 75-79 (ed. 8, 1910). 


In section 2 the notations are de- . 


The “law of propagation of error” is a tool that 
physical scientists have conveniently and frequently 
used in their work for many years. No claim is 
made here that it is the only tool or even a suitable 
tool for all occasions. ‘‘Data analysis” is an ever- 
expanding field and other methods, existing or new, 
are probably available for the analysis and inter- 
pretation for each particular set of data. Never- 
theless, under certain assumptions given in detail 
in the following sections, the approximations resulting 
from the use of these formulas are useful in giving 
an estimate of the uncertainty of a reported value. 
The uncertainty computed from the use of these 
formulas, however, is probably somewhat less that 
the actual in the sense that no function form is known 
exactly and the number of variables considered usually 
does not represent fully the contributors of errors 
that affect the final result. 


1. Statistical Tolerancing Versus Imprecision 
of a Derived Quantity 


1.1. Propagation of error formulas are frequently 
used by engineers in the type of problem called “‘Sta- 
tistical tolerancing.” In such problems, we are 
concerned with the behavior of the characteristic 
W of a system as related to the behavior of a charac- 
teristic X of its component. For instance, an engi- 
neer may have designed a circuit. A property W 
of the circuit may be related to the value X of the 
resistance used. As the value of X is changed, 
W changes and the relationship can be expressed 
by a mathematical function 


W = F(X) 


within a certain range of the values of X. 

Suppose our engineer decides on VW = wu» to be the 
desired property of the circuit, and specifies X =x» 
for this purpose. He realizes, however, that there 
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will be variations among the large lot of resistors he 
ordered, no matter how tight his specifications are. 
Let x denote the value of any one of the resistors 
in the lot, then some of the time x will be below 
xo, While at other times x will be above x. In other 
words, x has a distribution of values somewhat clus- 
tered about x. As x varies with each resistor, so 
does w with each circuit manufactured. 

If our engineer knows the mean and _ standard 
deviation (or variance) of x, based on data from the 
history of their manufacture, then he can calculate 
the approximate mean and variance of w by the 
propagation of error formulas: 


mean (w) = F(mean x), and 


2 


variance (w) = | var (x), (4.1) 


Z| 


where the square brackets signify that the derivatives 
within the brackets are to be evaluated at the mean 
of x. The approximations computed refer to the mean 
and variance of an individual unit in the collection of 
circuits that will be manufactured from the lot of 
resistors. The distribution of values of w, however, 
is still far from being determined since it depends 
entirely on the functional form of the relation between 
W and X, as mathematical variables, and the distri- 
bution of x itself, as a random variable. This type of 
approach has been used frequently in preliminary 
examinations of the reliability of performance of a 
system, where X may be considered as a multidimen- 
sional variable. 

1.2 Let us consider now the second situation under 
which propagation of error formulas are used. This 
situation is the one considered in Birge’s paper, and 
is the one that will be discussed in the main part of 
this paper. 

A physicist may wish to determine the “true” value 
wo of interest, for example, the atomic weight of silver. 
He makes n independent measurements on some re- 
lated quantity x and calculates 


Al 
Xn=—— (x1 txt. Xe) 
n 


as an estimate of the true value xo and 


l 


n—1 





= 


n = 
> (%i— Xn)” 
1 


as an index of dispersion of his measured values. The 
physicist is mainly concerned in obtaining an estimate 
w of wo, and of the standard deviation of # as a measure 
of precision of his result. He therefore computes by 
the propagation of error formulas: 


w = F (Xp) 


i ee 
var (w)= aX n 
o2= V Var (wv) 


Often he assumes that @ is distributed at least approxi- 
mately in accordance with the normal law of error and 
gives probability limits to the statistical uncertainty of 
his estimate ® based on the standard deviation calcu- 
lated (,,) and this assumption. 

Cramér [1946] has shown that under very general 
conditions, functions of sample moments are asymp- 
totically normal, with mean and variance given by 
the respective propagation of error formulas.2 Since 
X» is the first sample moment, the estimate @ will be 
approximately normally distributed for large n. Hence 
our physicist is interested in the variance (or the 
standard deviation) of the normal distribution which 
the distribution of F(x,) approximates as n increases. 
(Note that both estimators # and Var (w) are functions 
of n.) For n large, the distribution of &® can be as- 
sumed to be approximately normal and_ probability 
statements can be made about w. 

1.3. Hence, we have the two cases: 

(1) The problem of determining the mean and vari- 
ance (or standard deviation) of the actual distribu- 
tion of a given function F(x) of a particular random 
variable x, and 

(2) The problem of estimating the mean and vari- 
ance (or standard deviation) of the normal distri- 
bution to which the distribution of F(x,) tends asymp- 
totically. 

As examples of problems studied under the first 
case, we can cite Fieller [1932] on the ratio of two 
normally distributed random variables, and Craig 
[1937] and Goodman [1962] on the product of two 
or more random variables. Tukey. in three Princeton 
University reports, extended the classical formulas 
through the fourth order terms for the mean and 
variance, and propagated the skewness and elon- 
gation of the distribution of F(x) as well. These 
reports present perhaps the most exhaustive treat- 
ment of statistical tolerancing to date. 

From now on we shall be concerned in this paper 
with the second case only, i.e., the problem of esti- 
mating the mean and variance, or standard deviation, 
of the normal distribution to which the distribution of 
F(X) tends as n increases indefinitely, and hence 
also the problem of using approximations to the 
mean and variance computed from a finite number of 
measurements. Since the mean and _ standard 
deviation are the parameters that specify a particular 
normal distribution, our problem is by its very nature 
less complicated than that of statistical tolerancing 
where the actual distribution of the function may 
have to be specified. We shall, however, utilize 
formulas given in Tukey’s reports to check on the 
adequacy of some of the approximations. 


w 


(1.2) 


? A brief summary is given in paragraph 2.2. 
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2. Propagation of Error Formulas 
2.1. Definitions and Notations 


(1) X, Y, Z in capitals stand for the mathematical 
variables to be measured; x, y, z in lower cases stand 
for the measured values of these variables; xj, yj, zx 
with subscripts stand for the particular values of the 
ith measurement on x, the jth on y, and the Ath on z, 
respectively. 

(2) W=f(X, Y, Z) is a continuous function of the 
variables X, Y, Z, with derivatives 

OW Pw 
aX’ aXaY’ 


etc. 


(3) All derivatives appearing in square brackets, for 
OW | [oW : , 
example EA! 54! stand for the values of these de- 
ax oY | 
rivatives evaluated at the means of x and y, if known, 
or at the sample averages of x and y, if the means are 
not known. 

(4) In order“to emphasize the fact that the mean M, 
variance o? and other population parameters are usu- 
ally not known, we list here symbols for both the 
estimators of population values and the population 
values. For a particular set of values of x, the values 


computed from these estimators are estimates, or 


computed values of these estimators. 











Estimators of parameters Corresponding population 
parameters 
12 , 
r=- VS x; M, (mean= first moment) 
ne 
] n 
s2 =—— YF (x; -— x)? o. (variance = second 
a-l#& sentral : 
i=1 central moment) 
L fog Nea? 
— {Xx?7-= 
n—1 | n 
] n 
Sry = Syr ak 5 (xi —X) (yi— Y) Ory = Cyr (covariance) 
= i=l 
bie (xi) (Lyi) 
= Li 
n—1 | n 
ions saab eee 
Fry = ys Pry (correlation coefficient) 
S0Su 
it o, (standard deviation of x 
about M,) 
l Or none 
St=— Sr o;= (standard deviation 
Vn Vn of the average X, 
or standard error) 
Ss Ox si: Si > ans 
Vy i — (coefficient of variation or 
x M, relative standard 
deviation) 





In addition, we use |Ax| to denote the bound for pos- 
sible systematic errors on the measurements of x. 
The bound of these errors, unknown in sign, is usually 
established or conjectured by the experimenter and 
its value is not based on the measurements in hand. 


2.2. General Theorem and Remarks 


As mentioned briefly in paragraph 1.2, the propaga- 
tion of error formulas are special applications of re- 
sults obtained in the study of properties of distributions 
of functions of sample moments. Doob [1935], Hsu 
[1949], and others have investigated the limiting 
distribution of functions of sample means relating to 
hypothesis testing. Curtiss [1943] derived the limit- 
ing means and variances of the several functions of 
variables in connection with transformations used in 
the analysis of variance. Cramér, in chapters 27 
and 28 of his classical treatise, proved two theorems 
and also discussed the asymptotic properties of distri- 
butions of functions of sample moments in detail. 
For convenient reference we shall phrase his theorems 
and remarks in terms of functions of sample averages, 
to serve as a basis of justification for the use of propa- 
gation of error formulas. 

THEOREM (Cramer, pp. 366, 352-356) 

If, in some neighborhood of the point X= M,, Y = M,, 
the function F(X, Y) is continuous and has continuous 
derivatives of the first and second order with respect to 
the arguments X and Y, the random variable & = F(x, y) 
is asymptotically normal, the mean and variance of 
the limiting normal distribution being given by: 


mean W=F(M,, My) (2.1) 


var = [5] a+ |e +2 |] Ey n 99) 





REMARK 1]. (Cramér, p. 367) 

It follows from this theorem that any functien of 
sample averages is, for large values of n, approxi- 
mately normally distributed about the value of the 
function determined by the mean values of the basic 
variables, with a variance of the form C/n, provided 
only that expressions (2.1) and (2.2) yield finite values 
for the mean and the variance of the limiting 
distribution. 

REMARK 2. (Cramér, pp. 367, 415, also Doob, Hsu) 

In general, the constant C in the expression of the 
variance’ will have a positive value. However, in 
exceptional cases C may be zero, which implies that 
the variance is of a smaller order than n-!. Then 
some expression of the form 

n’{i—F(M,, M,)}, p> :. 
may have a definite limiting distribution, but this is 
not necessarily normal. 
REMARK 3. (Cramér, pp. 366, 213-214) 

The function F(x. ¥) may be asymptotically normal 
even though the mean and variance of F(X. 7) do not 
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exist, or do not tend to the mean and variance of the 
limiting normal form. Generally, if the distribution 
of a random variable w depends on a parameter n, 
and if two quantities M and o can be found such that 


the distribution function of the variable ——— tends 
o 


to ®t) (normal distribution function with mean zero 
and standard deviation one) as n>, we shall say that 
w is asymptotically normal (M, oa). This does not 
imply that the mean and the standard deviation of 
w tend to M and oa, nor even that these moments 
exist, but is simply equivalent to saying that for any 
interval (a, 6) not depending on n, 


lim Prob. (M+ao<w< M+ba)=®(b)— Pa). 


n> co 


EXAMPLE: If x is from a continuous distribution with 
positive mean and a finite variance but with positive 
probability that some x can take negative values, then 
the function In x is not even defined for all values of x, 
and therefore the mean of the function In x does not 
exist; yet where the mean of x has a positive value, 
(2.1) and (2.2) give the mean and variance of the lim- 
iting normal distribution. 


2.3. Propagation of Error Formulas 


Fortified with the general theorem stated in the 
preceding paragraph, we shall proceed to derive the 
traditional propagation of error formulas in an ele- 
mentary manner, making some comments and as- 
sumptions that may be of interest. It will be helpful, 
however, to explain first what is meant here by the 
term “random error” in a measurement process. 


a. Random Errors 


In a measurement situation, we consider random 
errors typically to be the sum total of all the small 
negligible independent errors over which we have no 
control —interpolation in reading scales, slight fluctu- 
ation in environmental conditions, imperfection and 
nonconstancy of our senses, etc. Thus for a stable 
measurement process, we find that: 

(1) The measured values do follow a distribution, 
with small errors occurring more frequently than larger 
ones, and with positive and negative errors about 
balancing one another, and 

(2) there is no obvious trend or pattern in the se- 
quence of measurements. 

Let us denote the ith measurement of x to be 


xi=M,+ 6; 


where M, is the mean of all measurements for the 
measurement process, and e¢; the random error of 
measurement x;. Then for condition 1, we assume 
A;: The distribution of errors is symmetrical and 
bell-shaped, with mean zero and standard deviation 
Or, OF 


mean €¢;=0 
mean x;=M, 
var x =var €; 
= mean €?=o0%. 


And for condition 2, we assume 
A»: The errors in the measurements of x; (i= 1. 2. 
. n) are statistically independent: in particular 
these errors are not correlated or associated in any 
way, 1.e., 


mean (¢; - €)) = 0, L#j. 


ae e 
Thus for x =-—(x,;+x2+. . .+-x,), the mean of x is M,. 
n 


Furthermore, 


= _& tet. 


ep 
x—M,= . 





n 
By definition, the variance of x is 


€: tert. 


' 9 
eS sa | 
n 


=5 {mean (> ] T mean (> ce) 


core / 
i=1 (FJ 


2 


mean (x — M,)?= mean ( 





l be 
=— \" mean (€;)? + y 3 mean (€€;) ¢° 
n- i 
FJ 


Using assumpiions 4, and A», we obtain 


~ bs 
Var {x)——o-; 
il 


or the variance of the average of n independent meas- 
as : : oe 

urements is — of the variance of an individual meas- 
n 


urement.* 


Here the average x is a linear function of the in- 
dividual x’s, and the exact expressions of mean and 
variance of an average in terms of that of the individual 
values are well known. For functions that are not 
linear in the x’s, we expand the function about the 
mean of x by the Taylor series, and assume that the 
function in the neighborhood of the mean can be ap- 
proximated by the lower order terms. For example, 


let 
W=F(X, Y), 
x=M,+€e:, 
y=M,+€,, 


3 If, however, the measurements are not independent, then this formula is incorrect since 
the means of products (€€;) are not equal to zero. In that case let 


ao 

Z 
mean (€€) =pyo2 , and p= S. pi/nin— 1), then var (x)= — {1 +(n— lp}. 

n 


iej 
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where each of €; and €, satisfies assumptions A; and 
As, then we can write 


F(x, y= F(M,z, M,)+ Big at BA ey 


aX oY 
1 [[eF) 4.6 [ oF PF | 
+o (Lie) ¢+2 [sar] «+ [Fe] a] 


+ terms of higher orders in €; and €y. 


(2.3) 


Or, neglecting terms of higher order than €% and €3, 


F(x, y)— F(M,, M,) = Ei €xt+ Bi Ey 


l 4 or 5) oF o?F ge. 
+5 ‘Ka ers Ben Ex€yt 3 4 é| 


Since the means of €, and €, are 0, if we take the mean 
on both sides, 





mean {F(x, y)—F(Mz, M,)} 


ariibca a+2|25| ont |S ot (2.4) 


Thus the mean of a function of values always differs 
from the value of a function of means by a quantity rep- 
resented by (2.4), approximately. If the function of 
means F(M,, M,) is the value of interest, then to 
approximate F(M,, M,) by the mean of F(x, y) would 
introduce an error, or bias, the magnitude of which 
depends on the functional form, the variances of and 
the covariance between x and y. If, however, we 
use the function of averages, F(X, 7), then 


mean w= mean F(x, y) 


= F(M,, M,) 


1 PF oF - o?F Ory oF 
+3 (Sale +21 Fo n +(&* oa, (2.5) 


and the bias is only 1/n times that of the mean of the 
function of individual values. When n becomes large, 
this bias tends to zero, and (2.1) results. 

This bias can be calculated by (2.5) and compared to 
the standard deviation of #. In practice, if ao, and co, 
are small, the bias is often of a magnitude that is 
negligible. 

To propagate the variance, we note that if €, and €, 
are small in the sense that the second and higher order 
terms in (2.3) can be collectively neglected in compari- 
son to terms involving €,; and €, only, then 


Pe = a . [aF oF 
Fix. y) F(M,. M,,) = bite He 





and the variance of F(x, y) is, appreximately, 


4 F oF , 
mean {F (x, y)—F(M,, M,)}?= mean|| 5 | €-+ ai} 


OF oF oF || oF 
= | — 2 » 
Bal os + |S] 3+ 2 || lon 


And for @= F(x, y), the variance of @ is 


nya [EP oes [EPs EYE] 24 er 
varia) = [5] n (3) a {5 | vin 


the limiting form of which is (2.2). 
Finally, if 0%, 02, and ory are not known, we substi- 

















tute their estimators in formulas (2.5) and (2.7). 
resulting in: 
mean(@) = F(M,. M,)+ , Ea : 
2 (LaX?] n 
0?F | s? o-F |s 
is a+ 2) | 2 
Ea n Ee (2.8) 
and 
F}?s2 [oF]? 82, [aF ][aF]s 
ar(Q) = —— — — cal 3] a padi ry 9 
var | n Fa ia Bala a (2.9) 


If we assume further that the random errors in meas- 
urements of x and y are independent, 1eno,,=0. and 
the terms involving 0, in (2.5), (2.6), aud (2.7) vanishes. 
If this is the case, the terms involving s,, in (2.8) and 
(2.9) should also be dropped. This reduced version 
of the formula for independent x and y. 


dF |? 02  foF Po 
inate ie Ea aaa Fra ey n° 


is of the form given in Birge’s paper and in other text- 
books on statistical analysis of data [Mandel, 1964, 
pp. 72-76]. 

For W= F(X, Y, Z), there will be three variance and 
three covariance terms in (2.5) and (2.7). Extension 
to more than three variables presents no new problems. 


(2.10) 


b. Extension to More Than One Function of the Variables 
Let 
U=2(X, Y, 2), 
and 


V=h(X, Y, Z). 
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Then in addition to the above formulas, we have 


s oe x2 A ey al |o ; + | al a : 
Cur = oxX ox OF ay ay vy az aZ J; 
al Eveeva ea Eoaeal 
lax a) aY ax [PruFrPy 
+(e) +(e 
{ IY AZ aZ ak 


+ FS 2+ Fi Pall YO 
\LaZ aX inaelU 


Expression (2.11) may be convenient to use to get 
oiw) where IE =F(U. 1), and l 
tions of X. }. and Z. 





| | puso yo 


1) 


and J are known func- 


c. Some Frequently Used Formulas 


For convenience. a few special formulas for com- 
monly encountered functions are listed in table 1 with 
ve) assumed to be independent. These may 
derived from the above formulas. 


be 


2.4. Systematic Errors 

By a systematic error we mean a fixed deviation that 
is inherent in each and every measurement of x in a 
particular sequence of measurements. If the magni- 
tude and direction of the systematic error are known, a 
correction can be made such that W,= x». or the mean 
of the sequence of measurements is equal to the value 
sought after. If the sign of the systematic error is not 
known and the magnitude of the error can be only 
estimated to be within some reasonable bound | Av}. 
perhaps by experience or judgment. then V/, is within 
the limits x) — Av and x») + Av. 

For a function of two variables W = F(. Y. Y) then, a 
bound | Aw} for the systematic error in W is given by: 


_ || oF 1 2F |, 
Aw |= | Fea FE Jay}. 


assuming, as before, that Ax and Ay are small such that 
second and higher order terms in Ax and Ay are collec- 
tively negligible in the Taylor series expansion. Since 
ordinarily we do not know the signs of Av and Ay, we 
have no choice but to add the absolute values of the 
two terms together. even though the signs of the values 
of the partial derivatives evaluated are known. (If the 
signs of either Ax or Ay is known, this information, of 
course, should not be ignored.) If these derivatives 
are evaluated at the point x and y, then the random 
components of error of x and y are required to be small 
so that these derivatives take approximately the same 
values as when evaluated at x» and yo. 


(2.12) 


When there are a number of systematic errors to be 
propagated, one approach is to take | Aw| as the square 
root of the sum of squares of terms on the right-hand 
side of (2.12), instead of adding together the absolute 
values of all the terms. This procedure presupposes 
that some of the systematic errors may be positive and 
the others negative, and the two classes cancel each 
other to a certain extent. 

The treatment of inaccuracy due to systematic errors 
of assignable origin but of unknown magnitudes is 
discussed in detail in section 4.2 of Eisenhart [1963]. 
Since there is no generally accepted standard method 
for combining several systematic errors, Eisenhart 
advised and we quote 

“Therefore. methods for the “com- 


bination of errors’ should indicate explicitly which of these (or an 
alternative method) he has used.” 








anyone who uses one of these 


Information on the source and magnitude of each con- 
tributing elemental systematic error is, of course, also 
essential. 


3. Practical Accuracies at the Various Stages 
of Approximations 


From the preceding sections we observe that 
there are three stages of approximations: 

In the Taylor series expansion (2.3), terms higher 
than the first partial derivatives are considered to be 
negligible. 

(2) w is approximately normally distributed for large 
n. Is.the normal distribution still a good approxima- 
tion for small n? 

(3) If o% and o? are known,we obtain o# from (2.7), 
and we can use this value to construct a confidence 
interval * about « with the desired level of confidence 
(approximate) based on normal theory. If o% and o? 
are not known, and s®.and s? are cale ulated from a small 
number of measureme nts, what can we say about 
w using var(w) calculated from (2.9)? 

To get some numerical feeling for the closeness of 
these approximations, we shall simplify matters by 
making the following assumptions which do not seem 
to be too restrictive in measurement situations: 

B,: x and y are normally and independently distrib- 
uted, with the ratio M/o not less than 10.4 

B.: The functional forms used are the well-behaved 
ones that do not possess derivatives assuming 
unreasonably large values when evaluated at the 
averages of the individual variables. 

Thus for linear functions, such as 


W =AX+ BY, 
the second and higher derivatives vanish, 
exact. 


The adequacy of these approximations is studied in 
paragraphs 3.2 and 3.3 below. In paragraph 3.4 sug- 


and (2.6) is 


$4See Natrella [1963]. sec. | to 7. 

*For notational convenience, 
subsequent sections. 
substitution. 


also chs. 2 and 3. 
the symbols w. x. y. +, dy, ete., are used in this and the 
The corresponding symbols for the average could be used by straight 
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) TABLE 1. Propagation of error formulas for some simple functions 
re 


nd Approx. formula for var (1) Term to be added if x and 5 
te Function form of w* (x and y are assumed to be are correlated, and a 
statistically independent) reliable estimate of o>. 
s—. can be assumed 
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cycil (* + *) 
sin = cos|-—— 
| sin 2 Ss 3 





sin4 


x ) sin? - 


m (sj and s> in radians) 


100 52! coefficient w? _ (not directly derived | 
»d x of variation 2(n—1) from the formulas) ' 
in percent) 


++ 











It is assumed that the value of «@ is finite and real, e.g.. y # 0 for ratios with y as denominator. x > 0 for Vx and In x 
Weighted mean as a special case of Ax + By, with o, and o, considered known. 
* Distribution of w is highly skewed and normal approximation could be seriously in error for small n. 
* See, for example, Statistical Theory with Engineering Applications, p. 301, by A. Hald (John Wiley & Sons, New York. N.Y.. 1952) 
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gestions are made on the use of the standard deviation 
calculated for w when the standard deviations of x and 
y are not known. Readers may wish to go directly to 
paragraph 3.5 for a summary of the conclusions. 

3.2. For x. y 


trary F(x. y), the first correction terms to (2.6) are 


aF \[a?F 3, [OF | [oF . 
—— = ae o> a Bee ~vo | YyOuns 
ax |, ax2}”°" Lay JLay? )"" 
where y is a measure of skewness of the distribution.® 
Therefore these terms equal zero for x. y symmetri- 
cally distributed, a condition satisfied by assump- 
tion B,. 
The next order of correction terms involve a}, a} 
and oo? and are usually negligible compared to terms 
in (2.6). These terms are 


1 (AE Ip t+ [?F][*l rye] | 
Pari: x3 | vi iba | 
o7F aFl]2 
eral = ~ noi+ |= (T,- h| 


“4 
+{ Flare arav| *larsrilori}o%* 


For functions involving powers of x and y less than 
three, some of the partial derivatives also vanish. 
For example, if W=XY, the only nonzero term of 
this order is 0?-03, or 





(3.1) 





Pe 








Var (w) = M303. + M203, + o40%. 


The contribution of 30% is less than 1 in 200 if M/o 
is larger than ten. 

For functional forms such as quotients, roots, and 
logarithms, the accuracy is usually adequate since 
powers of the means of the variables appear in the 
denominators of the partial derivatives. 

For the exponential function W=e*, the variance 
of w as given by (2.6) is 


Var (w) =e*o?, 


whereas the exact formula® for the variance of w, 
when x is normally distributed, is 


Var (w) = e**e2™(e"’ — 1) 
o* . oa ) 


=e"%erM (o+F a+ ar 31 >. 


= e?Mg? fer (+3 += m1 + 31 a+. . .)} 


>For definition of y and | see (3.3) 


® See, for example, The Lognormal Distribution, p. 8, by J. Aitchison and J. A. C. Brown, 
Cambridge University Press, 1957. 


independently distributed and arbi- 


Here the variance of w as given by (2.6) underestimates 
the true variance by the factor given in the brackets, 
and the approximation could be seriously in error. 
(Note, however, the “‘exact”’ formula is correct only 
if x is exactly normally distributed. If x is only ap- 
proximately normally distributed, then both for- 
mulas are approximations.) 

For specific functions, formulas (3.1) and (3.2) given 
in Tukey’s report can be used to check on the ade- 
quacy of the approximation. We quote Tukey’s con- 
clusion in this respect: 

“The most important conclusion is that the classical propagation 
formula is much better than seems to be usually realized. Examples 
indicate that it is quite likely to suffice for most work.” 

3.3 Next we look into the adequacy of the normal 
approximation. For this purpose we will define the 
first four central moments of the distribution of w as 
follows: 


mean (w—M,,)=0 
mean (w—M,)?=oc 
mean (w— M,.)3 = yo? 
mean (w— M,)4=To". (3.3) 
If w is normally distributed, 
lowing Tukey, we shall define 


skewness = yo", and 
elongation = l'a* — 3a?: 


then beth skewness and elongation are equal to zero 
when w is normally distributed. 

If x and y are normally distributed as assumed under 
B,, then in general w=F (x, y) is not normally dis- 
tributed unless the function form -is linear. By a 
procedure similar to that used in the last section, the 
coefficients of skewness B; and excess B:2 of w can be 
calculated where: 


[skewness w]? 





a 


[var w]* 


elongation w 
[var w]? 





B2.= +o. 


If B: is close to zero and 2 is close to 3, the normal 
approximation may be considered as adequate. 

The terms up to order o* in the propagation of 
skewness for w=F (x, y), with x, y independent, are 


ay ae 2 3 
skewness w = sr yro} + EI Yyo% 
3 (OF PHF] a — nes 
2 S| Fl (.—l)od 
oF |? | 0?F ; 
lav] Fate ~ Des 


oF | [oF | | _0?F — 
bi Ea Ei Ea we 





(3.4) 


270 


y=0, and [=3.  Fol- 








A SRT TT 





For x, y normally distributed, only terms of order 
ot remain. If we take w=xy again as an example, 
then 


skewness w = 6M,M,07%.03, 


36M2Mio%+0%4 
[Mjo?. + M20? + 0307 ]* 





fi, = 


Neglecting o%.0? in the brackets in the denominator, 
and taking M/a=10, B; is computed to be 0.045. 
Hence, for &@=x y, where x and y are averages of four, 
the coefhicient of skewness is reduced by a factor of 
four or equals 0.011 approximately. 

Similarly, terms up to order o* for the elongation 
of w=f(x, y), with x, y independent, are 


oF} oF]! 
slongati ihe ee = Dot +1— ,—3)o4 
elongation u Ea (C,—3)o4 + Be (Ty—3)o} 
(3.5) 
which is zero for x, y normal. 


elongation w 
Hence 2. = ——— 








— >+3=3, and no correction 
(Variance w)” 


for elongation is necessary here. 
If we look up a table’ of flercentage points of dis- 
Mw.) . 
with given 
on 
B, and B., we note that the changes of values are rather 
sensitive to B,; and much less so to Bo. Thus the co- 
efhcient of elongation is usually not as much a source 
of worry in the normal approximation as is the coefh- 
cient of skewness. 


tribution of the standardized variate 


Formulas (3.4) and (3.5) and the table of percentage 
points allow us to check how good the normal approxi- 
mation is for a given number of measurements in the 
variables x and y. Table 2 gives some examples of 
results of such calculations. 

3.4 The third approximation concerns the use of the 
sample variance s* as an estimate of the population 
variance o?. If we know the precision of the proc- 
esses for the measurements of x and y, i.e., we know 
oa, and a,, O can be computed from (2.7) and a con- 
fidence interval about w can be constructed with the de- 
sired confidence coefficient | — a@ by using the table of 
the normal probability integral. If o, and oy, are not 
known, then even if 6% can be computed from (2.9), 
the constants to be used for constructing a confidence 
interval with confidence coefficient 1—a will be dif- 
ferent from those for known a. 

To offer some guideline to the solution of this prob- 
lem, we again assume measurements on x and y to 
be independently and normally distributed. If the 
number of measurements is large (a rule of thumb 
could be n> 30). then (2.7) can be used assuming 
o*.. 7, and oy, are known. 


See Table 42. Biometrika Tables for Statisticians. Vol. 1, edited by E. S. Pearson and 
H. O. Hartley, The University Press. 1958 — Also, pp. 79-84 


Of course one can always compute the half-widths 
of the respective 100 (1—a@) percent confidence inter- 
vals for M, and for M, by the use of the Student’s ¢ 
statistic, and use (2.12) to get the half-width of the 
interval for M,. i.e., set 


Sr Sy 
— and, Ay= ti, a a 
Vn ; op '-' VE 


Ax=t 


(eee BY oar 


and use (2.12) to get Aw. Then the interval w+ Aw 
is a confidence interval for M, for a confidence co- 
efficient of at least (1—a). This procedure, however, 
may be criticized on the ground of gross inefficiency 
in using the data. 

We may write (2.9) as 


var (%) = Ais? si i ost 


| [oF |? | [oF }? 
where A, =— |—>| and A»=~- |=>] are two constants. 
n | ax hay 
For given degrees of freedom for s;. n—1. and Sy. 
Ais? 
k-1. + 
AS. + Ags7z 
statistic have been tabulated’ for confidence coef- 
ficients of 0.99, 0.98, 0.95, and 0.90. The interval 


and given ratios of values of a “v” 





w+tvu Vs? + des? (3.6) 


is a confidence interval with confidence coefficient 
le 

These tables, however, do not contain values for 
“» for n and & less than 10, 10, 8, and 6 for the re- 
spective confidence coefficients, and hence cannot 
be used for smaller samples. In addition, they are 
useful only for two independent variables x and y. 

Alternatively Welch [1947] has proposed the use 
of “effective degrees of freedom” for the estimated 
variance of @ of the form 


var (%@) = LAjs?. 
The effective degree of freedom f is computed from 


(SAjis?)? 


X(A?s3/fi) 


F (3.7) 


where /; is the degrees of freedom for? s?. 

In general f will be fractional. The ¢ value with f 
degrees of freedom can be found or interpolated from 
the ¢ table and the confidence interval computed as 


M=Et . Or. 
Ua oe 


“See Table 11. Biometrika Tables for Statisticians, Vol. 1: also Further critical values 
for the two-means problem, W. H. Trickett. B. L. Welch. and G. 5. James. Biometrika 43, 
1956, pp. 204-5 


*1f s? is computed from n; measurements, the degrees of freedom is n;— | 
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TABLE 2. 


x. ¥ independently distributed, with y = 0. | 


— 
| 
u | Skewness from (3.4) 
| | 
T T 
| 
| | 
{x+ By 0 | 
ee | 
| Co 
Nan | 6M,M, — 
n° 
| 
n=] | 
j | 
n } | | 
n=10 | 
| 
ao 
e Me 
n* 
} 
| 
n=10 
| 
oe ig 
Y 162\P —; 
n= | 
n=4 
YW.3 (ot oo? 
x 6(—) (+4) 
: \M, Ms MEM 
n 10 
s of 
nx —— + — 
M* n? 
n 10 
- o* 
Ce 3e = 
nr 
Exact when x a ’ ormally distributed 


The approximate confidence intervals computed by 
the use of effective degrees of freedom were found to 
check the exact confidence intervals given by (3.6) 
very well over the range of the latter. 

3.5 In summary, the following may be concluded 
for practical purposes: 

(1) Terms of order higher than o? in the propagation 
of error formulas for variance, (2.6) and (2.7). can be 
neglected if (a) the standard deviations are small in 
comparison to their respective means, and (b) the 
second and higher order partial derivatives evaluated 
at the means do not give rise to abnormally large num- 
bers. This is usually true in the field of physical 
science, since errors of measurements are usually 
of the order of | part in 1000, or parts per million: 
furthermore. the functional forms used are usually 
the well-behaved ones. 

(2) The normal approximation will be adequate for 
large n, or if, in addition to (a) and (b) above, (c) the 
individual variables can be assumed to be normally 


Departures from normal approximations 
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3. and (Wio)= 10 


> ~~ Vi, 
Percentage point of —->— 
On 


Bi computed | 


| Lower 2.5% | l pper 2.5% 
| | 
| T 
| | 
0 1.96 + 1.967 
| | 
| } 
| | 
15 | 
100n 
0.045 | 1.84 + 2.06 
O11 | 1.91 | 2.01 
0045 1.93 | 1.99 
9 | 
100n | 
0.009 1.90 | + 2.00 
36 
100n 
0.09 1.80 | 2.09. 
| 
18 | 
100n 
| | 
018 1.89 03 
| 
3) 
100n 
0.009 1.90 | + 2.00 
9o~ Depends on o and n (both skewness and B, 
n underestimated for a/ Vn > 0.2). 
distributed. For particular functions, the approxi- 


mate values of the coefficients of skewness and elon- 
gation may be calculated and Pearson’s table can be 
used to check the adequacy of the approximation. 

(3) For the case where the standard deviations of 
the individual variables are unknown, and are esti- 
mated from the data, confidence intervals for the 
estimate f can be constructed either by the use of 
tabulated values of the “v” statistic or by the use of 
effective degrees of freedom. These confidence 
intervals can be considered as a form of “precision 
limits” in the sense that if one makes the same sets 
of measurements a large number of times under the 
same conditions, and constructs the confidence in- 
tervals each time by the same procedure, then a large 
proportion of the intervals so constructed, 100 (1 — a) 
percent, will bracket the mean of all these sets of 
measurements. When only one set of measurements 
will be made, the probability is 1 — a@ that this interval 
will bracket the mean. 
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4. Reporting of Results 


4.1. Suppose a set of measurement data is avail- 
able, and, by using the appropriate propagation of 
error formulas, the following are obtained for the 
quantity of interest, wo: 

(1) The estimate of wo, w, based on n values of XY, 
etc.3 

(2) the estimated standard error of @, G », and asso- 
ciated degress of freedom f; 

(3) limits to the systematic error in w, Aw. 

The estimated standard error of # gives a measure 
of precision of the experimental results, or a measure 
of scatter of the values of w from the average value of 
M, for repeated performance of the particular experi- 
ment. But this measure of precision does not indicate 
at all how close this average value is to the value wo 
intended to be measured. The estimation of limits to 
the systematic error is an essential part of an experi- 
ment and need not be discussed here [Youden, 1961}. 
One may remark generally that ‘systematic errors 
usually do not pose a serious problem when the 
“imprecision” is large, since these systematic errors 
are, so to speak, “swallowed up” by the random errors. 
The systematic errors, however, play an important 
role when the precision is excellent and is of about the 
same order of magnitude as the systematic error. In 
that case, it is essential that the systematic error, or 
errors, be reported separately from the imprecision 
part of the reported value, as measured by the standard 
error, or the confidence intervals, computed. 

In scientific literature, it is not uncommon to come 
across expressions of results in the form of M+e, 
where “M” is an average of some kind and ‘“‘e” repre- 
sents the uncertainty of “M” in some vague sense. 
This type of reporting proves to be most frustrating 
from the reader’s point of view. From the context 
alone the reader cannot possibly infer whether, ‘“e” 
represents probable error, 3-sigma limits, systematic 
error, or some combination of random and systematic 
errors. As a consequence, the quality of the results, 
and the validity of inference drawn from these results, 
are to a large extent left to the judgment and guesswork 
of the reader. Hence, the writer owes to himself, and 
to his reader, to specify clearly the meaning of ‘‘e” 
as he uses it. In particular, the number of measure- 
ments from which the measure of random error was 
computed and the manner in which the systematic 
error was estimated are both essential elements of the 
reported value and need to be included. 
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A footnote explaining the role of “e”’ is often very 
helpful. Several examples are given below: 


“In the expression of the form M+e, M is the average and e is the 
standard error of M based on n measurements (or based on v degrees 
of freedom).” 
“The indicated uncertainty limits for M are overall limits of error 
based on 95 percent confidence limits for the mean___and on allow- 
ances for effects of known sources of possible systematic error es 
“The uncertainty given represents 3-sigma limits based on the 
current accepted value of the standard deviation, known sources of 
systematic errors being negligible.” 


Chapter 23 of Natrella [1963] ‘“‘Expressions of the 
Uncertainties of Final Results” gives a thorough dis- 
cussion on this topic, and is an excellent reference 
for all physical scientists who have occasion to report 
numerical results of their experiments. 
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The Apparent Thermal Radiation Properties of an Isothermal 
V-Groove With Specularly Reflecting Walls 


Richard B. Zipin* 
(July 26, 1966) 


The reflection of a parallel beam of light incident in a V-groove is treated for the case when the 
walls are specularly reflecting. Methods of calculating apparent reflectances and absorptances 
are given. Experiments have been performed to test the theoretical calculations, and the agreement 
has been found to be satisfactory. 


Key Words: Reflectance, specular walls, thermal radiation, V-grooves. 


1. Introducticn 


Recently there has been a great deal of work re- Le 
ported in the literature on the problem of predicting 
the apparent thermal radiation properties of cavities 
and rough surfaces when the properties of ideally 
smooth surfaces of the same material are known. 
An interesting problem that has not been previously 
treated very completely is reflection by an isothermal 
V-groove, whose walls are specularly reflecting. 
Sparrow and Lin [3]! have reported some results 
for the effective absorptance of such a groove. Their 
method of analysis is followed here using a simpler 
procedure, extended somewhat, and results are 
given for effective directional reflectance. 

The main limitation on the results reported here is 
due to the assumption that the V-groove has specular 
walls. This assumption is valid when the wave- 
length of the incident radiation is much less than the 
characteristic dimensions of the groove, and at 
the same time much greater than the character- 
istic dimensions of the roughness of the groove 
walls. Experimental results reported by the author 
[4] have shown this to be true. 





zs 








2. Reflection Within a V-Groove With 
Specularly Reflecting Walls 


The problem to be considered is illustrated in 
figure 1. The V-groove has an opening angle of 26, / 
a depth A, and is infinite in the direction normal to 
the plane of incidence. The method of calculation F ’ 
is based on the fact that light reflected from a spec- IGURE 1. at ats ra specularly reflected within a 
ular surface appears to come from an image located 2 
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behind the surface at a distance equal to the dis- 
tance ef the source from the front of the specular 
surface [2]. The number of reflections, n, that an 
incident beam makes at the walls of the groove be- 
fore it leaves determines the apparent thermal radi- 
ation properties exhibited by the groove, and _ this 
may be determined by counting the number of mir- 
ror images of the first surface that are crossed by 
an undeviated beam. A _ single beam is _ incident 
in the groove at an acute angle @ measured with re- 
spect to the normal to the groove opening, and leaves 
the groove at an acute angle x measured with re- 
spect to the same normal. 


The number of reflections is given by 


where w& is defined in figure |, and the square brackets 
indicate that the quantity enclosed within them should 
be rounded off to the nearest higher integer. This is 
because there can never be a fraction of a reflection, 
and all beams incident in a groove must make at least 
one reflection. 

There are two separate cases to be distinguished 
here: 

(1) If 6 > 8, an incident beam can only strike the 
far wall on entering the groove, and, 

(2) if @ < 4 an incident beam may strike either wall 
on entering the groove. 

For beams incident on the far wall, the law of sines 


gives 


_h/cos o (h—yp/cos 0 
sin(@+q)  sin(7w—wW—O—@) 








(2) 


where y is the vertical distance from the groove open- 
ing at which a beam strikes the wall, and the subscript 
f indicates a quantity associated with a beam incident 
on the far wall. Now. 


. ) ‘ 
W= 7 —O—gd—sin (1%) sin(o+¢) (3) 
\\ oh 
and by geometry 
xXr=(— I)“ ar — hb — 2). (4) 


For beams incident 
nis used, and 


on the near wall, the subscript 


A cos @ 


sin(@—o& 


(h—- Yulcos 





sin(—7wt+w, +6+ 0) 





































so that “< 
. Vn . 
W,=7—O—gdt+sin (a) <ino—4)} (6) 
v/ 
and in this case 
Xn=(— 1)" “a+ b—2n,6). (7) 


A value of x obtained from either of eqs (4) and (7) 
which is equa! to or greater than 90° is obviously for- 
bidden unless 6=90°. 

By the use of eq (1), with eqs (3) and (4) or eqs (6) and 
(7), the reflections undergone by any beam incident in 
the groove may be determined. 

If the groove is fully illuminated and @ < 6 both 
walls will- have energy directly incident on them, but 
if @ > 6 only the far wall will receive energy by direct 
incidence. In the latter case the far wall will be di- 
rectly illuminated from y=0 to y= ymax. Where Vmax 
is the value of y at which the incident beam that just 
grazes the lip of the groove on entering strikes the far 





wall. To calculate y,,,,. the law of sines is used again, 
and 
h/cos 0 (h— Ymax)/cos 0 8 
a oe a eS need — (5) 
sin(7— 6— 9) sin(d— @) 
Therefore, 
2h cosd sind , 
Vmax [ee SASS Ta ( 9) 


~sin(0+@) 


and when @=8@, Vmax =/ as required. 

All beams incident on both the near and the far walls 
of the groove will undergo the same number of reflec- 
tions and leave the groove at the same angle as either 
the y=0 or the y= Ymax beams.* The subscripts 0 and 
max denote quantities related to the beams incident 
at y=Oand y= Ymax respectively, and y, will denote the 
critical value of y at which n changes from no to nyax 
and x changes from Xo to Xmax- A beam incident at 





y- will just graze the lip of the groove on leaving. For 
the far wall, the law of sines is again used, and 
h/cos 6 _(h—ye,pleos 6 , 
sin(@+)  sin(a7—2no, ,A—O— od) (10) 
so that 
sin(2ny, ,8+ 6+ d) 
Yes hi | (11) 


sin(@+ ) 





*For the near wall, y,, «x Must be A 
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For the near wall, 


Alcos 0 (h—Ye.nleos 0 








— — — +») 
sin(@—@) sin(7w7—2ny.,0—0+ 0d) (12) 
and so 
sin(2ny,,8+ 6—@) 
ee — ®)\. (13) 
| sin(@— d) j 


If ¢6=0, there is no distinction between the near and 
far walls, and eqs (11) and (13) give identical results. 

The angle yx will be said to be negative if it is meas- 
ured in the same direction from the normal as dis. and 
positive if it is measured in the opposite direction. 
When }> 86, x is positive if n is even and negative 
when nis odd. 


3. The Calculation of the Apparent Thermal 
Radiation Properties 
¢ 

If the average reflectance of the groove walls is p. 
and the walls are opaque, then the average absorptance 
of the walls is a=1—p. When a single narrow beam 
of energy per unit time /; is incident in the groove at 
an angle @ > 6, it will strike the far wall, be reflected 
n times within the groove, and then leave at an angle 
x. x will have one of the two values Xo or Xmax. accord- 
ing to the value of y at which the beam is incident. 

On the first reflection of the beam, (1 — p)/; is ab- 
sorbed and p/; is reflected. On the second reflection, 
(1—p) pl; is absorbed and p?/; is reflected. This con- 
tinues until the final, nth reflection, when (1 — p)p""'/; 
is absorbed and p"/; is reflected. The energy-rate of 
the beam leaving the groove after the n reflections is 
therefore /,.= p"/;. 

If the effective reflectance is defined as the ratio of 
the energy reflected at the angle x to the energy inci- 
dent at the angle @, it is given by 


! ? 
Pei. X)= 7 =p (14) 


where the value of n used must correspond to the 
value x considered. 

If 6< 9, a single narrow beam may strike either 
the near or the far wall of the groove. In this case 
x will have one of the four values, Xo.¢, Xmax.f+ Xo. 
OF Xmax.n Equation (14) may be used to find the 
effective reflectance. 

Now consider the groove to be fully illuminated. 
If @>86, energy is incident on the far wall from 
y=0 to y=Ymax- Let the incident beam carry energy 
per unit time /; per unit width in the plane of incidence 
and normal to its direction of propagation, or equiva- 
lently /; sin (6+) per unit width along the groove 
wall, so that the total energy incident in the groove 
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per unit time is 
I(d > ™=T; sin (0+ b)¥max/cos 0. (15) 

Recalling eq (9), and substituting for ymax, this becomes 
(hd > 0)=21];h cos @tan O. (16) 

The fraction of this energy that leaves after no reflec- 


; ie : —_— 
tions at an angle x» is —— p"o, while the fraction that 
Y max 


Ve 





leaves after 2 max reflections at Xmax is ( i= . p"max, 
max 


As before, the effective reflectance at an angle x 
will be defined as the ratio of the energy leaving the 
groove in the direction x to the total energy incident 


in the groove. In this case we have two such 

reflectances, 
— Ve n x 
Pet. Xo) = —— po (17) 

”’ max 
and 
_ Ye ~ 
Perl @.Xmax) ot Lae — he (18) 
) max/ 


Since the total energy reflected by the groove is 
the sum of these two reflectances multiplied by the 
incident energy. the total energy absorbed is 


{1 — ped b. Xo) — Perhb.Xmax) } (db > 4) (19) 
and so the effective absorptance is 


Qerfl o> A=1—-— Pert .Xo) — Pert ®.Xmax)- (20) 


This expression is only a function of @. 

If 6< @ and the groove is fully illuminated, energy 
is incident on both walls from y=0 to y=A. If the 
incident beam carries energy per unit time /; per unit 
width in the plane of incidence and normal to its 
direction of propagation, then it will have /; sin (0+ ) 
per unit width along the far wall and /; sin (6—@) per 
unit width along the near wall. The total energy 
incident on the far wall is /; sin (0+ 0) h/eos # and on 
the near wall is /; sin (@—0) h/cos 6, so that the total 
energy inci: i in the groove is 


l(b < #=21;h cos d tan 0=1(h > 0). (21) 


From the above it is easily seen that the fraction of 

the total energy incident in the groove that is directly 
sin (8— ) 

iS a cor eae ena ae 

2 sin A cos 

fraction of this energy that leaves after no,,, reflections 


incident on the near wall and the 


. Vevn ' . 
at Xo,n is = plo.n and that leaves after nmay.n reflec- 
1 

































. . Vewn = . ‘ » 
tions at Xmax.n IS (1 —%=2) p"max.n. The fraction of 
2 
the total energy incident in the groove that is directly 
sin (6+ ¢d) 


incident on the far wall is - The effective 


2 sin Ocos & 
reflectances are, therefore, given by 


sin(O@-@)  Ve.n 





n 
Me: Me. SS Ths 22 
Peo. X 2sin @cos @ h (<4) 
sin(@— d) Va Ge 
orl O. aw = a "pescpaticad 2: 
Pon Xm 2 2 sin 0cos h f (23) 
sin(O0+@) Yer no - 
Die) L "0.5 (24) 
2 sin 0cos d& h 
and 
sin(@+ ) Yes es 
orl OD. : jbwere aaa ae ymax. fi 5 
Pett Xmax. J 2 sin 6cos h f io) 


These results depend on the angles only, and are 
independent of h. 


The effective absorptance is 
Add < A= 1 — porld. Xo n)— Perk. Xmax. n) 


— pet, X04) — perd @. Xmax.s)- (26) 

If d= 6, eq (26) gives the same result as eq (20) for 
the @ >@case. This is expected because when d= 8@. 
no energy is actually directly incident on the near wall 
of the groove. If @6=0, there can be no distinction 
between the two walls of the groove, and they will 
both have the same quantity of energy incident on 


them. In this case 
Pett b.Xo.n) = Petdb,Xo.9) (27) 
and 
Peth®sXmax.n) = Pett! P.Xmax. s)- (28) 


4. Comparison of the Theoretical and Experi- 
mental Results 


Measurements have been made to test the theo- 
retical results presented above. The angular distri- 
bution of energy reflected in the plane of incidence 
by carefully prepared V-grooved surfaces, which were 
obtained from Bausch and Lomb, Inc., was studied 
as a function of angle of incidence. Details of the 
apparatus and the specimens may be found elsewhere 
[4]. 

Figure 2 shows the positions of the peaks in the 
reflected energy distribution as calculated with eqs 
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FIGURE 2. The positions of the peaks of energy reflected by the 


V-grooves with specular walls. 


(4) and (7), for V-grooves having included angles of 
90°, 120°, and 150°. The experimentally determined 
positions are also shown in the figure, and they agree 
very well with the theoretical results. 

In figures 3, 4, and 5, the theoretical and experi- 
mental values of the reflectance ratio, R. at the peaks 
of the reflected energy distribution are plotted as a 
function of angle of incidence, for the three values 
of groove angle. The reflectance ratio is defined as 


R= Leth Pix) (29) 
Pett bj,X;) 

In each figure, the point at which the experimental 
value of R was made unity by choosing the measured 
value of Pepldj,xj) as the reference is i :dicated. 
At these points the theoretical values of perddj,x;) 
are unity. 

Due to the nature of the apparatus, many of the 
reflected peaks could not be observed. The boxes 
in the figures are used to indicate which theoretical 
point the experimental points are to be compared 
with, wherever there is a possibility of confusion. 
The theoretical points at which the reflectance ratio 
vanishes at either grazing observation (x=90°) for 
nongrazing incidence, or at nongrazing observation 
for grazing incidence (@=90°) are the limiting cases 
obtained from eqs (4) and (7), but peaks cannot exist 
exactly at these angles, as is proven by Beckmann [1]. 
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FIGURE 3. Comparison of theory and experiment for V-grooves 
with specular walls, 20=90°. 


Except for a few cases of large deviations between 
the experimental and the theoretical values of R. 
figures 3, 4, and 5 show that reasonable agreement 
exists between them. The agreement is much better 
at 0.5m than at 1.5, since the assumption that the 
walls of the grooves reflect specularly becomes better 
as the wavelength is decreased. Theoretical curves 
for cases of more than one reflection were calculated 
using p=0.5. 


5. Discussion 


The assumption that the average reflectance of 
the groove wall surfaces is known is only rigorous 
when the wall reflectance is not a function of the angle 
of incidence on the wall. If the reflectance is a func- 
tion of the angle of incidence on the wall, the calcula- 
tions would have to be repeated, and térms such as 
p" would be changed to a product of n different values 
of p, each one the reflectance of the groove wall at 
the angle of incidence of one of the n reflections. 
Except for this the calculations are valid for mono- 
chromatic radiation of any wavelength, or for total 
radiation if the properties of the groove wall are 
independent of wavelength.* 


‘It should be noted that the only wavelength dependence in the analysis is due to the 
value of p used. 
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FIGURE 4. Comparison of theory and experiment for V-grooves 
with specular walls, 20= 120°. 
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FIGURE 5. Comparison of theory and experiment for V-grooves 
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If the groove walls have properties that are not 
independent of wavelength and calculations are 
desired for total radiation, another complication 
enters the situation. Let p(A) be the wall reflectance 
as a function of wavelength, and /(A) be the spectral 
distribution of the source, then any function G(p(A)) 
can be calculated for total radiation by the use of the 
following expression: 


. _ SE(p(A)) dr 
{G(p(A)) } otal = fIQdxr 





(30) 


where the integrations are performed over all wave- 
lengths emitted by the source. 

Any desired calculations can now be made for the 
apparent absorptance or reflectance when the walls 
are specular. It should be noted that sometimes an 
equation will yield absurd results which have no 
physical significance. Should this occur, the appro- 
priateness of the parameters under consideration 
should be reexamined to ascertain their validity for 
the case under consideration. 

A conclusion which may be drawn from the results 
is that the effective absorptance increases and the 
effective reflectances decrease as @ decreases. In 
the case @ > 80, there will always be a region at the 
vertex of the groove which will be in shadow, i.e., it 
will not be illuminated either directly or due to internal 
reflections. In the case @< 6 the whole internal 
surface of the groove will be illuminated. 


The author expresses his thanks to his major profes- 
sor, Dr. Y. S. Touloukian for his advice and encourage- 
ment during the course of this work, and to the Purdue 
Research Foundation for support in the form of a 
David Ross Grant. Thanks are also due Dr. D. P. 


DeWitt for suggestions and stimulating discussions. 


6. Nomenclature 


G = Any radiation function, as used in eq (30). 

h =The depth of a V-groove. 

I =Energy per unit time carried by a beam of radia- 
i tion. 
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n =The number of reflections within a specular 
V-groove. } 

R= Reflectance ratio. 

y =A distance measured vertically from the lip of a 
V-groove. 

a = Absorptance. 

0 =The half-angle of a V-groove. 

\ = Wavelength. 

p = Reflectance. 

do = The angle of incidence. 

x =An angle of reflection. 

Ww =An angle defined in figure 1. 

|.|=Brackets indicating a quantity to be rounded off 
to the nearest higher integer. 


7. Subscripts 


A critical value. 
An effective or apparent value of some quantity. 


I 


c 
eff 


f, n =A quantity associated with the far or near wall 


of a V-groove, respectively. 

i,j =Jndexes. 

i,r =A quantity associated with an incident or 
reflected beam, respectively. 

max =A quantity associated with the maximum value 
of » on the far wall of a V-groove. 

0 =A quantity associated with y=0. 
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Elastic constants of synthetic single crystal corundum, 
W. E. Tefft, J. Res. NBS 70A (Phys. and Chem.), No. 4, 277-280 
(July—Aug. 1966). 

The elastic constants of synthetic single crystal corundum (aluminum 
oxide) were determined from 0 to 900 °K by a resonance technique. 


Properties of aqueous mixtures of pure salts. Thermo- 
dynamics of the ternary system: water-calcivum: chloride- 
magnesium chloride at 25 °C, R. A. Robinson and V. E. Bower, 
J. Res. NBS 70A (Phys. and Chem.), No. 4,305—311 (July—Aug. 1966). 
Isopiestic vapor pressure measurements have been made on the 
system water-calcium chloride-magnesium chloride at 25 °C. 
Activity coefficients have been evaluated for each salt in the presence 
of the other in systems of constant total ionic strength. The excess 
free energy of mixing has been calculated and compared with the 
values for analogous systems. 


Properties of aqueous mixtures of pure salts. Thermo- 
dynamics of the ternary system: water-sodium chloride- 
calcium chloride at 25 °C, R. A. Robinson and V. E. Bower, 
J. Res. NBS 70A (Phys. and Chem.), No. 4, 313-318 (July—Aug. 1966). 
Isopiestic vapor pressure measurements have been made to deter- 
mine some thermodynamic properties of the system: water-sodium 
chloride-calcium chloride at 25 °C. Equations are derived for the 
activity coefficient of each salt in the presence of the other when the 
total ionic strength is kept constant. The limiting conditions in 
very dilute solutions are considered and ‘trace’ activity coefficients 
calculated. A comparison with previous work is made. The excess 
free energy of mixing is calculated. 


Use of the consistency check in the vector verification 
method, A. D. Mighell, B. T. Gorres, and R. A. Jacobson, J. Res. 
VBS 70A (Phys. and Chem.), No. 4, 319-321 (July—Aug. 1966). 
The consistency check is a type of frequency check in which full 
use is made of the space group symmetry. It is designed for applica- 
tion to a symmetry map, to a tentative atom map, to the results of 
superpositions, or to any pseudo electron-density map. Like the 
frequency check, its purpose is to aid in the analysis of one of the 
above maps by assigning a relative weight (relative weight equals 
the number of other peaks that the given peak is consistent with) to 
the peaks in these maps. To determine the relative weight, each 
peak in the map is compared with each of the others to see if it is 
consistent. Two peaks are considered consistent with each other, 
if and only if, the set of unique vectors between the two is present 
in the Patterson map. The method is equivalent to doing a full set 
of symmetry superpositions at each peak in a symmetry map or a 
tentative atom map; counting the number of points in the resulting 
map; and outputting this number at the initial peak position. No 
assumptions other than that for space group symmetry are required 
to apply the consistency check to a symmetry map. The resulting 
map, which has been calculated in a routine fashion directly from 
the observed intensities, may then be conveniently used as the 
starting point for the structure analysis of a unknown compound. 


A rotation cylinder method for measuring normal spectral 
emittance of ceramic oxide specimens from 1200 to 1600 
°K, H. E. Clark and D. G. Moore, J. Res. NBS 70A (Phys. and 
Chem.), No. 5, 393-415 (Sept.-Oct. 1966). 

Equipment was designed and constructed for measuring the spectral 
emittance of polycrystalline ceramic oxide specimens in the wave- 
length region 1 to 15 mw and at temperatures of 1200, 1400, and 
1600 °K. Specimens consisted of small hollow cylinders that were 
rotated at 100 RPM in a furnace cavity equipped with a water-cooled 
viewing port. Emittances were determined by comparing the radi- 


-ance of the specimen to that of a blackbody at the same temperature. 


A series of six measurements (two measurements each on three 
specimens) was made on commercially pure specimens of alumina, 
thoria, magnesia, and zirconia. All four materials showed similar 
behavior in that the emittances were low in the shorter wavelength 
regions and high at the longer wavelengths. The temperature 
coefficients of spectral emittance were positive for all four materials, 
but varied with the material; those for thoria were appreciably higher 
than the other three. In all cases coefficients were greater at short 
wavelengths than at long. 

Error sources were investigated, and the measurement uncertainties 
associated with each source were evaluated. It was concluded that 
at 1200 °K the true values of emittance at 2 yw for the particular 
specimens measured might be as much as 0.032 lower or 0.012 higher 
than the reported values. The data were in good agreement with 
reflectance measurements of the same material and with recently 
reported absorption coefficients. 


On the approximation of functions of several variables, B. 
Mond and O. Shisha, J. Res. NBS 70B (Math. and Math. Phys.), 
Vo. 3, 211-218 (July—Sept. 1966). 

The purpose of this note is to point out how a certain type of approxi- 
mation to functions of one real variable gives rise to similar approxi- 
mations to functions of several variables. Information in the rapidity 
of convergence in the one dimensional case, yields at once corre- 
sponding information for the multidimensional case. 


Finding a rank-maximizing matrix block, A. J. Goldman and 
M. Newman. J. Res. NBS 70B (Math. and Math. Phys.), No. 3, 219 
220 (July—Sept. 1966). 

An algorithm is developed for the following problem: Given three 
matrices of respective dimensions s Xs, s Xt, and t Xs, to find a 
t X¢ matrix such that the (s +t) X (s +t) matrix formed from the four 
blocks has maximum rank. 


NBS standard frequency and time services, NBS Misc. Publ. 
236, 1966 edition (1966). 15 cents. (Supersedes previous editions). 
Detailed descriptions are given of eight technical services provided 
by the National Bureau of Standards radio stations WWV, WWVH, 
WWVB, and WWVL. These services are: 1. Standard radio 
frequencies; 2. Standard audio frequencies; 3. Standard musical 
pitch; 4. Standard time intervals; 5. Time signals; 6. UT2 correc- 
tions; 7. Radio propagation forecasts; and 8. Geophysical alerts. 
In order to provide users with the best possible services, occasional 
changes in the broadcasting schedules are required. This publica- 
tion shows the schedules in effect on January 1, 1966. Annual 
revisions will be made. Advance notices of changes occurring 
between revisions will be sent to regular users of these services 
upon request. 


Cooperation, convertibility, and compatibility among infor- 
mation systems: A literature review, M. M. Henderson, J. S. 
Moats, M. E. Stevens, and S. M. Newman, NBS Misc. Publ. 276 
(June 15, 1966), $2.00. 

The purpose of the study of the literature on which this report is 
based was to examine those problems in the field of documentation 
and in the operation of information systems which could possibly 
be solved or alleviated by some greater measure of cooperation, 
convertibility, or compatibility among systems. particularly those 
systems for handling scientific and technical information supported 
in whole or in part by the United States Government. An account 
is given of early developments and general background information 
about organizations active in cooperative documentation efforts. 
Current cooperative activities are then discussed in terms of dis- 
semination and publication of secondary sources, acquisition and 
exchange of publications, analysis and identification, systematiza- 
tion and terminology contro!, storage and search, and standardiza- 
tion. General problem areas, special problems created by changes 
in the nature of the documentary materials to be handled and 


281 





special problems raised by the prospects for mechanization are 
then discussed. A final section raises questions with regard to the 
implications for future progress. 

Doors as barriers to fire and smoke, H. Shoub and D. Gross, 
Bldg. Sci. Series 3, (Mar. 25, 1966), 15 cents. 

A study was made of means for improving dwelling unit entrance 
doors as fire and smoke barriers. Existing combustible doors and 
frames could be modified to enhance their fire resistance, but it 
did not appear practical to raise them to the level of rated commercial 
fire door assemblies. Fire retardant paints, except those consisting 
of heavy, reinforced, intumescent-type coatings, provided little or 
no increase in fire resistance. 

Several modification of existing doors were not effective in preventing 
the transmission of smoke. However, controlling the pressure. 
levels on both sides of a door, as by suitable venting, appeared to 
offer a means of reducing smoke penetration into an area. 

It is recommended that current methods of fire tests of doors, and 
criteria relating to their fire and smoke transmission be improved. 


Weather resistance of porcelain enamels. Effect of ex- 
posure site and other variables after 7 years, M. A. Rushmer 
and M. D. Burdick, Bldg. Sci. Series 4, (May 2, 1966), 20 cents. 
An exposure test of procelain enamels at seven representative 
sites in the continental United States was initiated by the National 
Bureau of Standards and the Porcelain Enamel Institute in 1956. 
After seven years all specimens were returned to the Bureau and 
the changes in gloss and color determined. These changes were 
found to be different at all exposure sites except Pittsburgh and 
New Orleans. The most severe changes occurred for specimens 
exposed at Kure Beach, 80 feet from the ocean, while the least 
change occurred for specimens exposed at Los Angeles. The 
differences in behavior of the specimens correlated with both the 
relative humidity and the pH of the suspended particulate matter 
at the different sites. 

{ direct relation existed between the acid resistance of the enamels 
and weather resistance. However, enamels of different types, such 
as enamels on aluminum and steel, having the same acid resistance 
did not necessarily show the same weather resistance. 

Comparison with enamel specimens exposed for seven years in an 
earlier test showed that porcelain enamels produced after the end 
of World War II were equally resistant to changes in gloss at the 
Washington, D.C.. site as those produced before the war. 

As a group, the regular glossy acid-resistant enamels on steel showed 
the best weather resistance of the various types tested. No cor- 
rosion of the base metal was noted for any specimen on which the 
initial coverage was complete. 


Some properties of the calcium aluminoferrite hydrates, 
E. T. Carlson, Bldg. Sci. Series 6 (June 1, 1966), 20 cents. 

Calcium aluminoferrite hydrates in two series, 4CaO-(Al,O03,Fe20s). 
n HO (hexagonal plates) and 3CaO-(Al2O3,Fe203)-6 H2O (isometric), 
were prepared from the anhydrous aluminoferrites by hydration in 
the presence of Ca(OH)2. The hexagonal phase was stable below 
15 °C, the isometric above 35 °C, in contact with solution. The end 
member, 2 CaQ-Fe2O3, did not produce an isometric hydrate, 
however. X-ray diffraction patterns of the hexagonal series in the 
19 H2O stage of hydration were indistinguishable from patterns of 
4 CaO-Al,O;-19 H.O. After drying to the 13-hydrate stage, slight 
differences in the patterns were observed. The isometric series 
shows a definite shift in x-ray pattern with change in Fe2O3/AbOs 
ratio, the unit-cell edge increasing from 12.573 for 3 CaO-Al,O3°6 
H:O to 12.716 for the member having Fe2O3/AloO3 close to 5. DTA 
traces show endotherms near 100 and 200 °C corresponding to de- 
hydration stages of the hexagonal series and one near 300 °C rep- 
resenting the first dehydration stage of the isometric phase. The 
hexagonal hydrates are converted to the isometric during the course 
of the DTA run. The peak attributed to the isometric phase shifts 
from 323 to 290 °C going from the pure aluminate to the member 
richest in Fe2O All the aluminoferrite hydrates, both hexagonal 
and isometric, reacted with CaSO, solution to give products of the 
monosulfate and trisulfate (ettringite) types, the latter predominating 
at later stages. The rate of reaction decreased with increasing 
FeO; content. 


Transistorized building blocks for data instrumentation, 
P. G. Stein, VBS Tech. Note 268 (May 28, 1966), 60 cents 
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A fourth family of etched-circuit logic cards has been developed at 
the National Bureau of Standards. : It closely parallels the design 
criteria and philosophy embodied in previous work, but differs in 
many respects. Each circuit description is separated physically 
from the others to facilitate use as an application manual. Specifica- 
tions, typical uses, graphic symbols, and reliability are discussed 
in brief. 


Key Words: Digital circuits, logic cards, etched-circuit package, 
digital instrumentation, logic design symbols, com- 


puter logic. 


Selected values of chemical thermodynamic properties. 
Part 2. Tables for the elements twenty-three through thirty- 
two in the standard order of arrangement, D. D. Wagnmn. 
W.H. Evans, I. Halow. V. B. Parker, S. M. Bailey, and R. H. Schumm. 
VBS Tech. Note 270-2 (May 6, 1966), 40 cents. 

Technical Note 270-2 is the second of a series of publications 
containing material prepared as a revision of the tables of Series | 
of National Bureau of Standards Circular 500, Selected Values of 
Chemical Thermodynamic Properties. This Note contains re- 
vised Tables 24:to 32, covering the elements silicon, germanium, 


tin, lead, boron, aluminum, gallium, indium, and thallium, and 
compounds of these elements following the Standard Order of 
Arrangement. 


Tables for the first twenty-three elements and their compounds, 
and a discussion of the symbols,.conventions, units of energy. and 
conversion factors used in this series are given in National Bureau 
of Standards Technical Note 270-1. 


Nuclear and radiation standards of importance to the Na- 
tional Atomic Energy Program, H. W. Koch, J. H. Donnert, 
W. W. Seng dr. G. Rogosa. and L. Rosen, VBS Tech. 283 
(Mar. 31, 1966), 35 cents: 

A oe approach to physical measurements in nuclear physics 
requires the establishment, availability, and use of standards. 
The standards of measurement may consist of standard cross sec- 
tions, foils, counters, and evaluated nuclear data. Because of the 
fundamental importance of measurement standards, the Nuclear 
Cross Section Advisory Group has examined the definitions, re- 
quirements, characteristics, and availability of nuclear and radi- 
ation standards of importance to the national atomic energy program. 


Spark source mass spectrograph program July 
June 1965, P. J. Paulsen and P. E. 
(Apr. 15, 1966), 25 cents. 

This report describes the activities of the spark source mass spec- 
trograph project for the period from July 1964 through June 1965. 


1964 to 
Branch, VBS Tech. Note 286 


The subjects covered include instrument modifications which 
resulted in lower backgrounds, modifications planned for the 


near future, an outline of the procedure used in an analysis, a 
description of data computation techniques, and experimental 
results obtained for some NBS standard reference materials. 
Photographs of the instrument with its modifications as well as 
photographs of the improved mass spectrum are also included. 


Measurement philosophy of the pilot program for mass 
calibration, P. E. Pontius, VBS Tech. Note 288 (May 6, 1966), 
30 cents. 

The Pilot Program for mass measurement is the result of a consid- 
eration in which the values produced are thought of as the prod- 
ucts of a mass measurement process. The collective perform- 
ance of elements of the mass measurement process results in 
establishing the process precision which, under certain conditions, 
can be described quantitatively by pertinent performance param- 
eters. The uncertainty attached to the product of the process, 
the measured value, is computed from these parameters and reflects 
the total performance of the process rather than the immediate 
measurement which might have produced the value. Interpre- 
tations of uncertainty and surveillance tests are discussed. The 
Pilot Program in mass measurement, whereby suitable process 
performance parameters can be established for precise mass 
measurement processes in other facilities, is discussed. 
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Status report, National Standard Reference Data System 
April 1966, ed. E. L. Brady, NBS Tech. Note 289 (June 1, 1966), 
50 cents. 

This report summarizes the status of the activities of the NBS 
Office of Standard Reference Data as of March 1966. It pro- 
vides a detailed review of the data compilation activities within 
the seven broad categories within the technical scope of the pro- 
gram: (1) nuclear data, (2) atomic data and molecular data, (3) 
solid state data, (4) thermodynamic and transport data, (5) chemical 
kinetics, (6) colloid and surface properties, and (7) mechanical 


properties. Plans for an information services program are also 
outlined. Certain problem areas are identified. The appendix 


includes a listing of properties with which the program is con- 
cerned and a listing of the advisory panels of the Office of Standard 
Reference Data. 


Key Words: Nuclear data, atomic and molecular data, solid state 
data: thermodynamic and transport properties, chemical kinetics, 
colloid and surface properties, mechanical properties; informa- 
tion services, standard reference data. 


Calculation of the admittance of a parallel plate capacitor 
containing a toroid-shaped sample, E. G. Johnson, Jr., VBS 
Tech. Note 336 (Apr 21, 1966), 30 cents. 

The theoretical admittance of a toroid-shaped sample ina parallel 
plate caggcitor system is developed. This calculation represents 


a first an’ .econd order perturbation correction to the usual quasi- 
static formula. The FORTRAN listing of the machine program is 
included. The analysis is done under the condition that the sam- 


ple is not in contact with the capacitor plates and can be described 
by a complex permeability, wu, and complex permittivity, €y. 

plex | ys plex 1 : 
Key words: Fortran, lossy materials, parallei 
permeability, permittivity, theoretical admittance. 


plate capacitor, 


Advances in ionospheric mapping by numerical methods, 
W. B. Jones, R. P. Graham, and M. Leftin, VBS Tech. Note 337 
(May 12, 1966), 45 cents. (Formerly the Central Radio Propaga- 
tion Laboratory of the National Bureau of Standards, now the 
Environmental Science Services Administration, Boulder, Colo.) 
Paper describes recent progress made at ITSA (formerly CRPL) 
to improve numerical methods for mapping and predicting iono- 
spheric characteristics used in telecommunication. Two major 
problems are considered: (1) tendency .of maps to smooth out 
physical properties of the ionosphere, particularly at low latitudes 
and (2) ambiguous values at geographic poles and resulting dis- 
tortions in immediate surroundings. The second problem is 
overcome by means of a universal time analysis. Significant 
improvement is made in the first problem by the use of a modified 
magnetic dip coordinate. A general description is given of the 
new procedures for forming numerical maps, including a number 
of illustrations. Also included are a discussion of the use of nu- 
merical maps in the revised form and a Fortran Program for evalu- 
ating numerical maps of ionospheric characteristics and for cal- 
culating the MUF(ZERO)F2 and MUF(4000)F2. 


The vaterite-type ABO; rare-earth borates, W. F. Bradley, 
D. L. Graf, and R.S. Roth, Acta Cryst. 20, Pt. 2, 283-287 (Feb. 1966). 
Approximate models are proposed for the room temperature and 
the high temperature modifications of vaterite-type rare earth 
borates. A low temperature model based on three-membered 
rings of borate tetrahedra is found to be reconcilable with the 
X-ray powder diffraction data, the optical properties, and the 
infrared absorption character. A high-temperature model con- 
taining triangular borate ions is relatable to the CaCQs; vaterite 
modification. 


Time-resolved electrical measurements in high current 
discharges, E. C. Casidy, S. W. Zimmerman, and K. K. Neuman, 
Rev. Sci. Instr. 37, No. 2, 210-214 (Feb. 1966). ; 

A method for measurement of the resistive component of the 
instantaneous voltage across a sample installed in a high voltage. 
high current circuit is described. 
the current permitted time resolved determination of electrical 
energy dissipation, power, and resistance of the sample. The 


Simultaneous measurement of 


system was calibrated calorimetrically, and measurements were 
made with exploding wire samples. 


Accuracy of standards and measurements, W. A. Wildhack, 
Acta. Imeko 4, 63-81 (1964). 

This paper reviews the accuracy of measurements and of cali- 
bration at the National Bureau of Standards, based on “taccuracy 
charts” prepared for each of various areas of measurement. Ac- 
curacy charts, consisting of plots of accuracy vs. magnitude of the 
measured quantity for each of the more precise standards, are 
shown in connection with a discussion of the present status and 
future goals for measurement and calibration capabilities. Meas- 
urement areas included are length, mass, temperature, pressure, 
vacuum, force, resistance, capacitance, inductance and _ voltage. 
The concept and utility of the accuracy chart itself are stressed, 
as a technique for a basic appraisal of capabilities and goals. The 
need for a more general agreement on expressions for accuracy 
is discussed, and recent work in this direction by NBS is reviewed. 


Experimental verification of superachromatism, R. E. 
Stephens, J. Opt. Soc. Am. 56, No. 2, 213-214 (Feb. 1966) 

To verify superachromatism experimentally a five-element prism, 
to deviate a beam of light approximately 30 degrees without dis- 
persion or distortion, has been designed and constructed from three 
types of Schott glass. When tested on a spectrometer, viewing 
a narrow slit with a telescope magnification of approximately 
30 times, no dispersion could be seen. Pointings on 7 isolated 
spectrum lines separated by a monochromator show a maximum 
dispersion of 3 sec, 1 part in 36000. 


Interpretation of Brillouin spectra, R. D. Mountain, J. Chem. 
Phys. 44, No. 2, 832-833 (Jan. 15, 1966). 

When a fluid consists of molecules having internal degrees of 
freedom weakly coupled to the translational degrees of freedom 
of the fluid (phonons), it is not correct to refer to the phonons as 
isentropic fluctuations in the density. This is relevant to the 
interpretation of Brillouin scattered light in terms of the prop- 
erties of the scattering fluid. An approximate expression for 
the ratio of the intensity of the unshifted component to that of 
the Brillouin components is presented and applied to CS». 


Spectral distribution of scattered light in a simple fluid, 
R. D. Mountain, Rev. Mod. Phys. 38, No. 205-214 (Jan. 1966). 
The spectral distribution of light scattered by density fluctuations 
in a dense, monatomic, one-component fluid is calculated from 
the time dependence of the density fluctuations predicted by the 
linearized hydrodynamic equations of irreversible thermodynamics. 
The results of Landau and Placzek are verified and a procedure 
for deriving correction terms is discussed with the dispersion in 
the velocity of thermal sound waves obtained as an illustration. 
Particular attention is paid to the critical region. The proper- 
ties of carbon dioxide are used to estimate the spectral distribu- 
tion of critical opalescence. A comparison is made between 
light-scattering and  sound-propagation experiments. Space 
dispersion near the critical point in the pressure and the thermal 
conductivity is examined briefly. Finally, some of the experi- 
mental problems involved in measuring the spectral distribution 
of the scattered light are discussed. 


Search methods used with transistor patent applications, 
J..R. Cornog and H. L. Bryan, Jr., IEEE Spectrum, pp. 116—121 
(Feb. 1966). 

As the world has become more technically oriented, the number of 
patent applications has .been increasing also too rapidly for the 
Patent Office to assimilate them comfortably with current tech- 
niques. When an application for a patent is received, it must be 
evaluated as to novelty by a specialist who searches the prior art 
for similar patents. Previously, all such searches were done 
manually, which meant that the examiner had to rely entirely 
on his knowledge and experience. In an effort to remedy the 
situation the Office has instituted mechanized search methods. 
In order to ascertain the differences in patterns of thinking associ- 
ated with manual and mechanized searches, a study was carried 
out in which a patent application in the transistor art was searched 
both manually and by a mechanized method. The mechanized 
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search in this case permitted more patents to be analyzed more 
quickly but, being completely literal, it does not allow for hunches 
or browsing. 


Composition, work and interrelation of international and 
national organization engaged in the standardization of 
dental materials and therapeutic agents, G. C. Paffenbarger 
and M. Kumpula. /ntern. Dent. J. 15, No. 4, 571-587 (Dec. 1965). 
There are two international bodies, two regional, and about 12 
national groups and agencies that are concerned with specifications 
for dental materials. On the international scene there is some 
jurisdictional conflict between the Federation Dentaire Interna- 
tionale (FDI), composed of national dental associations, and the 
International Organization for Standardization (ISO), composed of 
the national standardization bodies. The American groups have 
have been the most active for a longer period of time and, hence, 
their influence has been paramount. 


Dental applications of polymers: a review, G. M. Brauer, 
J. Am. Dental Assoc. 72, 1151-1158 (May 1966). 

The properties of acrylic polymers make them especially useful 
for denture base materials and plastic teeth. Plastic fillings have 
not been too successful due to their high coefficient of thermal 
expansion and lack of bonding to tooth structure. A cavity liner 
impervious to the mouth fluids that bonds to tooth structure is 
needed. If such a liner were available, the need for a filling ma- 
terial that matches the thermal expansion of the tooth would be 
less critical. 

For taking impressions, silicones and polysulfides have been em- 
ployed successfully. The possible application of other plastics 
such as epoxy, aziridino, polycarbonates and semi-organic poly- 
mers is reviewed. 


Dislocations and stacking faults in rutile crystals grown 
by flame-fusion methods, D. J. Barber and E. N. Farabaugh, 
J. Appl. Phys. 36, No. 9, 2803-2806 (Sept. 1965). 

Single crystals of rutile, grown by the Verneuil method, have been 
chemically thinned and examined by electron transmission mi- 
croscopy. The density and distribution of dislocations agree with 
that expected from etch pit data, but the density is at least two 
orders of magnitude below that for rutile films formed by oxidizing 
titanium carbide. The etching of dislocations and stacking faults 
is discussed. 

Some preliminary results of high-temperature studies in the micro- 
scope are reported. These include the discovery of a hitherto 
unknown slip system, namely < 100 > [001]. 


Effect of reactions with the atmosphere during fatigue of 
metals, J. A. Bennett (Proc. Sagamore Conf. Fatigue, Sagamore, 
New York, Aug. 15, 1963), Book, Fatigue: An Interdisciplinary Ap- 
proach, pp. 209-227 (Syracuse University Press, Syracuse, N.Y., 
1964). 

Data in the literature show that most structural metals are sub- 
ject to corrosion fatigue in normal indoor atmospheres. Pro- 
tective films of polar organic compounds can mitigate the effects 
of the environment on many metals, particularly during crack 
growth. Although a high humidity atmosphere is detrimental to 
the fatigue strength of aluminum alloys, recent investigations have 
shown that there is an initial period during which the presence of 
water vapor has no effect on the number of cycles to fracture. 
After this initial period, moisture has a marked effect on the life, 
apparently due to a reaction with the metal that results in the 
release of a relatively large volume of hydrogen. This change in 
chemical behavior is believed to result from a change in the nature 
of the surface deformation. 


Electrometric method of pH determination, R. G. Bates, 
Book, Standard Methods of Chemical Analysis, 6th Edition, IIIA. 
Instrumental Analysis, ed. F. J. Welcher, ch. 26, pp. 521-532 (D. 
VanNostrand Co., Inc., Princeton, N.J., 1966). 

This chapter discusses the nature of the definition of pH and the 
operational pH scale, describes the most useful pH-responsive 
electrodes and the most common reference electrodes and salt 
bridges, and sets forth the characteristics of modern pH instru- 
mentation. Careful instructions for standardization of pH assem- 


blies and for the measurement of the pH of unknown or test sol 
tions are given, the limitations of experimental pH numbers a 
described, and the mode of interpretation of the measured pi 
is outlined. 


Melting, flow, and thermal expansion characteristics of 
some dental and commercial waxes, M. Ohashi and G. C. 
Paffenbarger, J. Am. Dental Assoc. 72, 1141-1150 (May 1966). 
The literature on pertinent properties of waxes that may be suit- 
able in dental impression materials is meager. To provide more 
data, the melting point and range (ASTM D 87-57), some flow 
characteristics (ADA Specification No. 4), and the linear therma’ 
expansion (ADA Specification No. 4) were determined on 66 dent: 

and commercial waxes. Instead of plotting the time-temperatur. 
curve the more sensitive temperature-rate cooling curve was em- 
ployed. These curves were useful in identifying waxes and vari 
ations among lots or batches. The arrests in the curves were a: 
high as 141 °C. and as low as 41 °C. Supercooling was observed 
in six instances. The flow of wax cylinders (6 mm. high) under 
a 2,000-gram load at 30, 37, 40, and 45 °C. had a variable relation 
ship with melting points. The coefficients of linear thermal e» 
pansion ranged from 1,000 x 10-®/°C. (25-30°C.) to as low a: 

110 x 10-8/°C. (25-30°C.). 


Rubber and rubber products, M. Tryon and E. Horowit: , 
Book, Standard Methods of Chemical Analysis, 6th Edition, IIIb 
Instrumental Analysis, ed. F. J. Welcher, ch. 59, pp. 1664-1763 
(D. VanNorstrand Co., Inc., Princeton, N.J., 1966). 

A compilation of selected methods for analysis of rubber and 
rubber products by instrumental techniques. 


A precision pulse-operated electronic phase shifter and 
frequency translator, J. Barnes and A. Wainwright, Proc. IEEE 
33, No. 12, 2143-2144 (Dec. 1965). 

A reference signal is processed through a voltage variable delay 
line and a narrow band filter. After a voltage pulse of proper 
amplitude and shape is applied to the voltage variable delay line, 
the output signal from the narrow band filter has exactly one cycle 
of phase less (or more depending upon pulse shape) than the ref- 
erence signal. By suitable processing, it is possible to generate 
precise submultiples of the one cycle (360°) of phase and thus 
realize precise and stable phase shifts. 

This system has great utility in the generation of small but pre- 
cise frequency offsets such as are used in the generation of the 
coordinated universal time scales. To this end a rather compact 
unit has been constructed to offset the frequency of a 100 KHz signal 
from 0 to —500 parts in 10’ in steps of 50 parts in 10". 


Can infrared improve visibility through fog? A. Ashley 
and C. A. Douglas, Illum. Engr. 61, No. 4, 243-250 (Apr. 1966). 
The fundamental system parameters-atmospheric attenuation 
detector sensitivity, radiant intensity of source, and effects o 
background are examined to determine whether or not the use o. 
infrared sources and detectors, instead of visible sources and the 
eye, for marking airport runways and taxiways offer sufficient 
potential to warrant further consideration of their use. as airport 
marker “lights.” The evaluation is based on theory and previous 
experimental work. At present there are no advantages in the 
use of infrared sources and detectors as a replacement for the 
current visible sources and further consideration is not warrantec 
at this time. 


Condensate collection as a measuring technique for study- 
ing the cooling capacity of air-conditioning units, P. R. 
Achenbach and J. C. Davis (Proc. XI Intern. Congress of Refrig- 
eration, Munich, Germany, 1963), Book, Progress in Refrigeration 
Science and Technology Wl, 1281-1286 (Pergamon Press Inc. 
New York, N.Y., 1963). 

The condensate produced by an air-to-air heat pump operated as 
an air-conditioning unit was weighed simultaneously with the 
measurement of total cooling capacity by the refrigerant flow- 
meter method and the sensible and latent cooling capacities by the 
psychrometric method. Comparative values of the latent cool- 
ing capacities and of the total cooling capacities were in good 
agreement. The investigation indicated that a careful collection 
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,and weighing of the condensate is useful in evaluating other labora- 
Jory techniques and could be used to simplify field measurements. 


Cryogenic coil for megajoule energy storage, V. Arp, Proc. 
Intern. Symp. Magnet Technology, ed. H. Brechna and H. S. Gordon, 
pp. 625-629 (U.S. Atomic Energy Commission, 1965). 

A liquid hydrogen cooled induction coil is being constructed for 
short time energy storage of 10° joules. The coil is made of very 
pure aluminum having an electrical resistivity of about 3 x 10-9 
ohm-cm at 20 °K. The coil parameters are dictated by its pro- 
posed use for the Army Missile Command, and are as follows: 
inductance of 50u henries, charging time of 10 seconds or less 
from zero to peak current of 200,000 amperes, and maximum 
terminal voltage during discharge of 12,000 volts. The coil will be 
housed in a non-metallic dewar, and powered from batteries. De- 
tails of the coil and cryostat design and electrical circuits will be 
discussed. The basic design is capable of being pushed to higher 
energy storage without major difficulty. 


Design considerations for a Kossel microdiffraction camera, 
D. L. Vieth and H. Yakowitz, Rev. Sci. Instr. 37, No. 2, 206-209 
Feb. 1966). 

A Kossel microdiffraction camera has been adapted to the National 
Bureau of Standards electron probe microanalyzer. Complete 
design concepts for any Kossel camera are discussed and evalu- 
ated. The details of their conversion to the National Bureau of 
‘Standards Kossel microdiffraction camera are given. The major 
concepts adopted were: (1) Transmission camera, (2) Inclusion of 
microgoniometric capabilities, (3) Film cassette in air rather than 
vacuum, (4) Variable sourte to film distance of 6 to 11 em. Trans- 
mission pseudo-Kossel patterns of an aluminum crystal employ- 
ing Cu radiation are shown to illustrate the capabilities of the 
camera. 


Heating performance of air-to-air heat pumps at two Air 
Force housing projects, P. R. Achenbach, C. W. Phillips, and 
W. T. Smith, ASHRAE J. 70, 390-407 (1964). 

The heating performance of the split-type air-to-air heat pumps 
in five sample dwellings at one base and of the unitary air-to-air 
heat pumps in five sample dwellings at another base of the U.S. 
Air Force was investigated during a one-month period of cold 
weather at each site. The test results showed that the compres- 
sion systems of the 2.5- and 3.6-hp split systems provided 51 to 
82 percent of the heating requirements of the sample houses at one 
air base, whereas the compression systems of the 5-hp unitary 
systems provided 99 to 116 percent of the heating requirements 
at the other air base at‘ design conditions. The supplementary 
resistance heaters used about 10% as much energy as the com- 
pressor system during the tests. The coefficient of performance 
averaged about 1.8 at an outdoor temperature of 20 °F for the 
split systems and about 2.0 at an outdoor temperature of 14 °F for 
the unitary systems. The electric energy used by the heat-pump 
;compressor and resistance heaters averaged 52 and 55 percent of 
the total energy used in the sample houses at the two air bases; the 
water heaters used 25 and 30 percent of the total; and the remainder 
was used by the cooking ranges, clothes dryers, and miscellaneous 
devices. 


Kinematic locator for crystal alignment, B. G. Simson and 
R: D. Deslattes, Rev. Sci. Instr. 37, No. 3, 300-301 (Mar. 1966). 
!A kinematic locator with one rotational degree of freedom is de- 
scribed. This device makes possible the transfer of aligned 
crystals within 2 seconds of arc. 


Laboratory study of the effect of solar radiation on cool- 
ing load of stationary refrigerated vehicles, P. R. Achen- 
bach and C. W. Phillips (Proc. XI Intern. Congress of Refriger- 
ation, Munich, Germany, 1963) Book, Progress in Refrigeration 
Science and Technology Ul, 1631-1635 (Pergamon Press Inc., 
New York, N.Y., 1963). 

Solar radiation is a significant factor in the overall cooling load 
of refrigerated trucks used for distributing frozen food within a 
city because such deliveries are made principally during the day- 
time and involve frequent stops. A laboratory apparatus for simu- 
lating the effect of solar radiation on the cooling load of refrig- 
erated trucks is described. The apparatus consists of stationary 


banks of electric heating elements with adjustable power input 
mounted in fixed parabolic reflectors facing each of the longer 
walls and the roof of the vehicle under test. The absorption of 
solar radiation by a dark-painted truck was simulated by dupli- 
cating the rise in surface temperature observed on a similar ve- 
hicle on a bright sunny day. Studies made on three trucks indi- 
cated that the maximum increase in cooling load, due to solar 
radiation on stationary dark-painted vehicles, ranged from 13 to 
20 percent of the steady state value when the percentage increase 
was based on the highest average cooling load for four consecutive 
hours. The test apparatus could be used to simulate solar radi- 
ation effects on vehicles with other surface characteristics by 
selecting an appropriate pattern of surface temperature rise. 


Measuring the cooling load of refrigerated vehicles by free 
evaporation of liquid nitrogen, P. R. Achenbach, C. W. Phil- 
lips, and R. W. Penney, Suppl. Bull. Intern. Inst. Refrigeration, 
Commission VII, 3, 115-124 (1965). 

The total steady-state cooling load of an insulated structure, sim- 
ulating a truck body, was determined by two alternate methods, 
one, using a chilled brine coil mounted inside the structure, and 
the other, by release and evaporation of liquid nitrogen in the re- 
frigerated space. There was a significant difference in the values 
obtained by the two methods. The measured air leakage of the 
structure was shown to account for most of the difference between 
the two values. 


Thermodynamic property computations for systems analy- 
sis, J. G. Hust and R. B. Stewart, ASHRAE J., pp. 64-69 (Feb. 1966). 

Methods for calculating thermodynamic properties from equations 
of state and zero pressure properties are reviewed. The calculation 
of entropy of gases, saturated fluids, and liquids is considered in 
detail. Particular attention is given to the determination of the two 
phase boundary and thermodynamic consistency of the calculated 

properties. 

A method of iterative solution for implicit variables and its application 
to thermodynamic calculations are discussed. In addition methods 
of least squares are considered. Normal equations resulting from 
the usual minimization procedure are compared with corresponding 
equations of the constrained least squares problem. Discussion is 
also included on a technique for fitting related thermodynamic data 
simultaneously. 


Critical configurational entropy at glass transformation, 
S. S. Chang, A. B. Bestul, and J. A. Horman, Proc. VII Intern. Conf. 
Glass, Brussels, Belgium, 1965, pp. 26.1-26.15 (Gordon & Breach, 
New York, N.Y., 1966). 

When an equilibrium supercooled liquid transforms to a glass, the 
amount of configurational entropy which it possesses at the glass 
transformation is trapped into the glass. By two different methods. 
and a distinguishable variation of one of them, it is shown experi- 
mentally that the value of this configurational entropy at glass 
transformation for most glasses for which data exist to evaluate it, 
is near 0.7 cal deg-'! per “bead” as defined by Wunderlich. These 
beads are equivalent to molecular chain links for polymers. The 
above stated result is obtained first by evaluating configurational 
entropy from experimental data on molecular configurational re- 
laxation times, using a reinterpretation of the Williams-Landel- 
Ferry equation, for the temperature dependence of such relaxation 
times, in terms of configurational entropy rather than free volume. 
Secondly, the configurational entropy at glass transformation is 
evaluated independently as the calorimetrically determined residual 
entropy of the glass at absolute zero temperature. In addition the 
above results imply a fixed value for the ratio of the glass transforma- 
tion temperature to the reference temperature at which configura- 
tional entropy would vanish for an equilibrium supercooled liquid. 
This implication would provide a third method of evaluating con- 
figurational entropy at the glass transformation if the necessary 
reference temperature could be determined experimentally by a 
non-calorimetric method. The implications of the configurational 
entropy reinterpretation of the Williams-Landel-Ferry equation 
are examined. 

Field emission ultramicrometer, R. D. Young, Rev. Sci. Instr. 
37, No. 3, 275-278 (Mar. 1966). 

A simple, contact free, ultrasensitive distance and displacement 
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measuring instrument has been investigated. Experiments demon- 
strate that the instrument is capable of operating at spacings as 
small as a few hundred angstroms. Calculations indicate that 
distances of 10-* to 10-* cm and less can be reproduced to within 
about one part in 105. With suitable calibration, distance measure- 
ments in the 10-° to 10-® cm range can be expected to have ac- 
curacies limited only by available calibration techniques. The 
instrument would be most useful as a null or differential distance 
measuring device. It has the unusual property that resolution 
improves over several orders of magnitude as the null or measure- 
ment point is approached. Since the instrument contains no optical 
or mechanical lever systems, or delicately balanced bridges, it has 
inherent long term stability. Proposed applications include: (1) 
measurement of ball and hole diameters (contact free), (2) differen- 
tial thermai expansion cell, (3) mechanical vibration sensor, (4) sur- 
face profile delineator (contact free), and (5) surface contour delinea- 
tor (contact free). 


Microscopic visualization of dose distributions, W. L. Mc- 
Laughlin, Intern. J. Appl. Radiation Isotopes 17, 86—96 (1966). 
Colorless cyanides of triphenylmethane dyes when suitably ac- 
tivated can be made into films that become deeply colored upon 
irradiation with short-wave ultraviolet and ionizing radiations, but 
are not sensitive to near ultraviolet radiation or visible light. The 
response range for x rays, y-rays, and electrons is approximately 
10° to 108 rads. Experiments show that, because of their stability, 
low energy dependence, and high spatial resolution, these solid 
systems have excellent potential for visual inspection of high-level 
dose distributions on a microscopic scale. Since molecular excita- 
tions due to low-energy secondaries are most important to the over-all 
radiation effects, the relatively large sensitivity of such dye systems 
expected in the intermediate and far ultraviolet may represent an 
advantage in dosimetry. 


Multiple scattering corrections for proton range measure- 
ments, M. J. Berger and S. M. Seltzer, Natl. Acad. Sci.-Natl. Res. 
Council Publ. 1133, #5. Studies in Penetration of Charged Particles 
in Matter Nuclear Science Series Report 39, pp. 69-98 (1964). 
Three approximate methods are described which may be used to 
calculate differential and integral distributions of projected range, 
median ranges and curves of ionization vs. depth (Bragg curves) for 
protons traversing thick absorbers. One of the methods is based 
on the application of a multiple scattering detour distribution ob- 
tained by Yang; the second is due to Bichsel and Uehling; the third 
makes use of random sampling. Two types of statistical fluctuations 
are taken into account: (a) energy loss straggling in collisions with 
atomic electrons, and (b) the “wiggliness” of the track (detours) 
due to multiple Coulomb scattering by atoms. All three methods 
were used to extract the values of the proton c.s.d.a. range, r, and 
of the mean excitation energy, Jaa;, of the medium from Bragg curves 
in lead and copper measured by Mather and Segré and by Zrelov and 
Stoletov. They give results for r and Jaq; that are in good agreement 
with each other, but do not account entirely for shape of the meas- 
ured Bragg curves. 


Nuclear giant quadrupole resonance, R. Ligensa, W. Greiner, 
and M. Danos, Phys. Rev. Letters 16, No. 9, 364-367 (Feb. 28, 1966). 
The description of the giant quadrupole resonance of deformed 
nuclei is extended to include the coupling between the different 
collective modes. The interaction with the odd particle is also 
included. The results are compared with recent experimental 
photon absorption cross sections. 


Nucleation in polymers, F. Gornick and J. D. Hoffman, Ind. Engr. 
Chem. 58, No. 2, 41-53 (Feb. 1966). 

This article is a feview of current theories and experimental data on 
homogeneous nécleation and crystal growth with chain folding in 
linear polymers. 


A quick, direct method for the determination of activation 
energy from thermogravimetric data, J. H. Flynn and L. A. 
Wall, Polymer Letters 4, Pt. B, No. 5, 323-328 (May 1966). 

A quick simple method for determining energy of activation from 
thermogravimetric weight-loss vs. temperature data at two or more 
constant heating rates is presented. The activation energy is ac- 
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curately obtained trom the slope of a logarithm heating rate vs. 
reciprocal absolute temperature plot at constant degree of conver- 
sion. The independence of activation energy from degree of con- 
version and temperature may be tested by determining it at various 
values of constant weight loss. 


A standard reference material for Méssbauer spectrometry 
of iron and its compounds, J. J. Spijkerman, F. C. Ruegg, and 
J. R. DeVoe, Technical Reports Series 50, Applications of the 
Mossbauer Effect in Chemistry and Solid-State Physics, pp. 254-259 
(International Atomic Energy Agency, Vienna, Austria, 1966). 
The rapid development of the Méssbauer effect has resulted in a new 
spectrometric method for chemical structure analysis. The various 
sources used in Modssbauer spectroscopy requires a reference 
material to provide Méssbauer data on a uniform basis. A series 
of standard reference materials for Méssbauer spectroscopy will be 
made available from the National Bureau of Standards, U.S.A. The 
first in this series is a single crystal of disodium pentacyanonitrosy]- 
ferrate dihydrate, Na2{Fe(CN);NO]-2H2O. The crystals are supplied 
by a commercial manufacturer, and calibrated by NBSon a Méss- 
bauer spectrometer, using an optical fringe counting technique. 
The Mossbauer spectrum of a single crystal absorber, cut along the 
be plane, is a well resolved, symmetric doublet. The center of this 
doublet is defined as the zero reference point for the differential 
chemical shift, with the containing absorber at 25.0 °C. The ab- 
sorber crystals supplied will be lem x lem, with 25.0 mg/cm? of 
natural iron. 


Crystal growth of bone mineral, W. E. Brown, Clin. Orthopaed. 
44, 205-219 (1966). 

Formation of hydroxyapatite in vitro occurs either (1) through direct 
precipitation, or (2) by the initial formation of octacalcium phosphate 
followed by hydrolysis in situ to hydroxyapatite. The two processes 
lead to crystals with distinctive morphologies. A third process for 
the formation of hydroxyapatite (alternate precipitation of octacal- 
cium phosphate and hydrolysis to hydroxyapatite on the unit-cell 
level) could lead to crystals with either morphology. Octacalcium 
phosphate plays a vital, albeit transitory, role in two of these proc- 
esses and allows the interpretation of properties of apatitic materials 
without assumptions regarding the existence of arbitrary and ill- 
defined crystalline species. The effects of fluoride ions, in par- 
ticular, on the properties of apatitic materials become understand- 
able in terms of octacalcium phosphate as an intermediary in the 
crystallization of hydroxyapatite. 

The role of octacalcium phosphate in the formation of tooth and bone 
is likely to be equally as great as that in vitro. Some of these 
possibilities are discussed for the first time. 


Infra-red spectra of hydroxyapatite, octacalcium phosphate 
and pyrolysed octacalcium phosphate, B. O. Fowler, E. C. 
Moreno, and W. E. Brown, Arch. Oral. Biol. 11, 477-492 (1966). 
Infrared spectra of hydroxyapatite, octacalcium phosphate, and 
pyrolyzed octacalcium phosphate were studied in the range 4,000 
to 400 cm~', and probable band assignments are given. Unreported 
bands were found for both octacalcium phosphate (13 bands) and 
hydroxyapatite (5 bands); the hydroxyapatite band at 631 cm~' was 
reassigned to the librational mode of the hydroxyl group. 

The pyrolysis reactions of octacalcium phosphate were studied by 
IR absorption, weight loss, and pyrophosphate analysis. The 
major products between 325 and 600 °C were hydroxyapatite and 
B—CaP2,0; and B—Ca;(POx)2. Simple dehydration was con- 
tinuous from 50 to about 400 °C. Formation of B — CazP20;, instead 
of y—CazP20;, was detected at a lower temperature (325 °C) than 
has been reported. The maximum amount of pyrophosphate formed, 
about 40% at 500 °C, was intermediate between the amounts pre- 
dicted by reactions postulated in other studies. From 700 to 
900°C, the pyrophosphate content approximated the expected values. 


Phase rule considerations and the solubility of tooth enamel, 
W. E. Brown and B. W. Wallace, Ann. N.Y. Acad. Sci. 131, Art. 
2, 690-693 (Sept. 30, 1965). 

The phase diagram is a particularly useful means of displaying 
solubility information because it makes it easier to take into account 
several variables at one time. The use of the phase diagram is 
illustrated in connection with the possible reactions that can take 
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place when NaF —H;PQ, solutions are used to treat enamel. A 
critical quantity proves to be the concentration of Na* ion in the 
treatment solution. 

Calcium and phosphorus ions may not diffuse through enamel at the 
same rates. A situation in which calcium diffuses faster than phos- 
phate would lead to increased concentrations of calcium, phosphate 
and hydrogen ions within the enamel. This in turn would lead to 


increased rates of diffusion and dissolution, thus helping to account. 


for subsurface dissolution of enamel. 


Soft tissue response to implants of gallium alloys and silver 
amalgam alloys, H. W. Lyon, R. M. Waterstrat, and G. C. Paffen- 
barger, J. Am. Dental Assoc. 72, No. 3, 659-664 (Mar. 1966). 
Twenty implants of gold-gallium alloys and 20 silver-amalgam con- 
trols were placed subcutaneously in ten NMRI-D strain rats for six 
months. Silver-amalgam control specimens elicited a mild, innocu- 
ous response. In contrast, gold-gallium alloy implants not only 
caused a severe foreign body response (fibrosis) from the host tissues 
but also fragmented and partially disintegrated. These results were 
exactly similar to those obtained in a previous experiment wherein 
palladium-gallium alloy implants were subjected to identical pro- 
cedures. Disintegration of the experimental alloys was probably 
due to corrosion of the AuGaz phase rather than an electrolytic 
reaction between this phase and the gold. 


A statistical model of flicker noise, J. A. Barnes and D. W. Allan, 
Proc. IEEE 54, No. 2, 176-178 (Feb. 1966). 

By the method of fractional order of integration it is shown that it is 
possible to derive the characteristics of flicker noise from those of 
“white,” band limited noise. A formal expression for the relation 
of flicker noise to white noise is given. An approximate method, 
amenable to the use of digital computers, is also given for the 
generation of flicker noise modulated numbers from random, 
independent numbers. 


Key Words: Statistical, model, flicker, noise. 


An automatic ellipsometer. Automatic polarimetry by 
means of an ADP polarization modulator III, H. Takasaki, 
Appl. Opt. 5, No. 5, 759-764 (May 1966). 

An analytical and experimental study of an automatic ellipsometer 
is reported. The polarized beam is modulated simultaneously by 
two ADP cells and the corresponding two signals are separated and 
are used to adjust polarizer and analyzer settings in any of the four 
zones. The reproducibility of the setting has been demonstrated 
to be of the order of .01 degree. 


Key words: ADP modulator, automatic ellipsometer, optics, 
polarimetry, thin film, ellipsometer. 


An intercomparison of atomic standards, R. Beehler, D. 
Halford, R. Harrach, D. Allan, D. Glaze, C. Snider, J. Barnes, R. 
Vessot, H. Peters, J. Vanier, L. Cutler, and L. Bodily, Proc. IEEE 
54, No. 2, 301-302 (Feb. 1966). ‘ 

A group of two cesium beam frequency standards and three hydrogen 
masers were brought together for extensive intercomparisons at 
the National Bureau of Standards. After thorough evaluation of 
each of the two types of frequency standards a preliminary value of 
the free-space frequency for the transition F = 1<— F=0 of the 
ground electronic state in hydrogen has been obtained. This 
preliminary value is 


1420, 405, 751.7860 + 0.0043 Hz 


The uncertainty of 0.0043 Hz corresponds to an uncertainty of 3.2 
parts in 10!?. 


Key words: Hydrogen maser, cesium, frequency comparison. 


Bone-air cancellation, E. L. Smith, N. S. Bowman, P. G. Weissler, 
and R. K. Cook, Bull. Lab. Electroacoustique, No. 9, 45-50 (Apr. 
1966). 

An experimental study using bone-air cancellation to correlate the 
sensation levels of the air and bone signals and the occlusion effect 


of audiometric earphones has been made. Attempts were also made 
to use cancellation in lieu of masking in bone-conduction threshold 
determinations. 


Key words: Air-conduction, audiometry, bone-air cancellation, 
bone-conduction, hearing, masking, psychoacoustics, 
threshold. 


Convenient safety cutoff device for water cooled equipment, 
J. J. Ritter and T. D. Coyle, Rev. Sci. Instr. 37, No. 4, 523 (Apr. 1966). 
This note describes a device to detect interruptions in the flow of 
cooling water through laboratory apparatus. A control circuit is 
provided to deactivate equipment and to shut off water supply to 
prevent flooding. 


Key words: Safety, flow-sensitive switch, cutoff, electric. 


Study suggests value of shared computers, A. E. Rikli, S. I. 
Allen, and S. N. Alexander, Mod. Hosp. 106, No. 5, 100-108 
(May 1966). 

A survey of real-time hospital information systems in advanced 
stages of development has been made and is reported briefly here. 
By definition, such a system is one that is in direct communication 
with portions of the hospital environment by receiving data, process- 
ing transactions, and returning results with sufficient responsiveness 
to affect patient care and administrative decisions near the time of 
occurrence of medical or hospital events. Eight hospitals qualified 
under the definition, and the general conclusion to be drawn is that 
more design studies and experimentation with equipment and tech- 
niques are essential to solving the special communications and 
information problems of hospitals. 


Key words: Real-time hospital information system, hospital data 
communications, remote data input, nursing station, 
laboratory, pharmacy, administration. 


The NBS-A time scale—its generation and dissemination, 
J. A. Barnes, D. H. Andrews, and D. W. Allen, /EEE Trans. Instr. 
Meas. IM-14.,, 228-232 (Dec. 1965). 

The NBS-A time scale has been generated at the Boulder Labora- 
tories of the National Bureau of Standards since April 1963. This 
time scale is based on the United States Frequency Standard, 
located in Boulder, Colorado, and is now disseminated from radio 
stations WWVB, WWV, and WWVH. 

The precision and estimated accuracy of the time scale are dis- 
cussed along with its comparison to the A.1 time scale of the Naval 
Observatory and TA, time scale of the Laboratoire Suisse de 
Recherches, Horlogeres, Neuchatel, Switzerland. The present 
technique for disseminating the NBS-UA time scale is discussed. 


Absolute determination of refractive indices of gases at 47.7 
gigahertz, A. C. Newell and R. C. Baird, J. Appl. Phys. 36, No. 12, 
3751-3759 (Dec. 1965). 

The refractive index of a gas is determined to high accuracy by 
measuring the change in the resonant frequency when the gas is 
admitted to a previously evacuated microwave cavity. This meas- 
urement was made at a frequency of 47.736 GHz, and is the first 
such measurement in this frequency range. 

The resonant frequency of the cavity is detected by means of a hybrid 
tee bridge circuit. When the cavity is matched to the wave guide 
the null at resonance is very sharp, and the frequency can be set 
to +5 parts in 10°. 

Two types of resonant structures were used in this measurement, 
a conventional cylindrical cavity and a Fabry-Perot resonator using 
one plane and one spherical mirror. Measurements were made with 
both structures on each of the gases, and the agreement between the 
two was within the limits of experimental error. Because the 
Fabry-Perot resonator has no side walls, air could be passed through 
it much more easily, and it is therefore preferred for use with moist 
air. 

The values of (n— 1) < 10® for the gases at 0 °C and 760 mm Hg. are 
as follows: Dry CO, free air . . . 288.1+0.1: Argon . . . 277.5+0.1: 
Carbon Dioxide . . . 493.9+0.4: Helium . . . 34.5+0.1: Hydro- 
gen... 135.8+0.1: Nitrogen .. . 293.8+0.1: Oxygen... 267.0+0.1. 
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International comparison of measurements at high frequen- 
cies, M. C. Selby. JEEE Spectrum 3, No. 1. 89-98 (Jan. 1966). 

The status of international agreement in measuring electrical 
quantities at frequencies above approximately 30 kHz, as of 1965, 
is analyzed. An attempt is made to tell (1) why it is necessary to 
agree internationally on measurements of high frequency electrical 
quantities, (2) for what specific quantities and related ranges (““meas- 
urands”) agreement is desirable at present, (3) what has been done 
in the past towards that end, and (4) what is being planned for the 
future. Aims and relative responsibilities and authorities of the 
International Electrotechnical Commission (IEC), of the Inter- 
national Radio Scientific Union (URSI), and of the International 
Bureau of Weights & Measures (BIPM) are briefly indicated. 


LF-VLF frequency and time services of NBS, D. H. Andrews, 
IEEE Trans. Instr. Meas. IM-14, No. 4, 233-237 (Dec. 1965). 

The United States Frequency Standard and National Bureau of 
Standards Time Scale are described and the techniques by which 
they are used to control the broadcasts from WWVB and WWVL, 
Fort Collins, Colorado are presented. A practical method for con- 
trol of frequency and time broadcasts from WWY,. Greenbelt, 
Maryland is described and actual results shown. 


The system of electromagnetic quantities at 30 kHz to GHz, 
M. C. Selby, Metrologia 2, No. 1, 37-45 (Jan. 1966). 

The system of electromagnetic quantities, units, and standards at 
frequencies of 30 kHz to approximately 1 GHz is briefly described. 
Its application at the National Bureau of Standards is used as an 
example. The system with some likely modifications seems in 
operation in several technically developed countries and is tenta- 
tively suggested for newly developing countries. Basic and derived 
physical quantities, and conceptual bases of standards and methods 
of measurement are described. Ranges of accuracies are given 
and suggestions for extensions of the system and for establishing 
validity via international comparison are pointed out. 


The system of electromagnetic quantities at frequencies 
above 1 GHz, R. W. Seatty. Metrologia 2, No. 1, 46-54 (Jan. 1966). 
The system of electromagnetic quantities at frequencies above 10° 
cycles per second is briefly described, giving emphasis to those 
techniques and standards which have been most thoroughly evalu- 
ated and tested. The quantities included are power, noise tempera- 
ture, reflection coefficient and impedance, phase shift, attenuation, 
and field strength. The measurement of frequency and wavelength, 
and of the microwave properties of materials is not within the scope 
of this paper. 

In addition to a brief review, some of the problems involved in 
obtaining a higher order of accuracy are discussed. Comparisons 
of standards of different nations are mentioned and the benefits to 
be derived from such comparisons are described. Areas needing 
further research. as well as those which would benefit from inter- 
national cooperation are designated. 

Although a highly detailed and technical presentation is avoided, 
numerous references are given. In addition, a selected list of 
definitions is appended. 


Key words: Electromagnetic quantities, microwave standards, micro- 
wave measurements, international comparisons of standards, coaxial 
connectors. 


A nomogram for determining azimuth and horizontal trace 
velocity from tripartite measurements, H. Matheson, Earth- 
quake Notes XXXVI, No. 1, 33-37 (Mar. 1966). 

The direction normal to a traveling plane wave front and the wave’s 
horizontal trace velocity at the earth’s surface, are easily determined 
by observing the wave at three noncollinear points. Reduction of 
the data involves the solution of several trigonometric equations. 
This note describes a nomogram which reduces this task to thirty 
seconds with ruler and pencil. 


Key words: Nomogram, tri-partite nomogram. 


Extended dislocation nodes in a silver-tin alloy, L. K. Ives and 
A.W. Ruff. Jr.. J. Appl. Phys. 37, No. 4, 1831-1837 (Mar. 15, 1966). 
A study of dislocation nodes formed in silver-tin solid solution alloys 
has been conducted using the technique of transmission electron 


microscopy. In alloys of 6 and 8 at. % Sn, both extrinsic and 
intrinsic stacking faults were observed at three-fold dislocation 
nodes. Burgers vector data obtained from diffraction contrast 
studies were in agreement with predictions based on geometric 
considerations. From a qualitative comparison of the node exten- 
sions, the extrinsic fault energy was found to be at least twice the 
intrinsic fault energy. Both fault energies were observed to de- 
crease with increasing solute concentration. At many of the nodes 
containing extrinsic faults, the dislocation curvature was small or 
negative; at the intrinsically faulted nodes, positive curvatures were 
always observed. It is believed that the dislocation constrictions or 
cross-linking dislocations required at the extrinsically faulted nodes 
were responsible for this difference. Motion under stress of the 
extrinsically faulted nodes was observed, consistent with the co- 
planar Burgers vector geometry. The general occurrence of both 
types of faults in low stacking fault energy FCC materials is 
predicted. 





Key words: Dislocation, alloys, transmission electron microscopy, 
stacking faults, dislocation nodes. 


On the definition of the stacking-fault: energy in binary 
alloys, R. DeWit and R. E. Howard. Acta. Met. 14, No. 3, 431-433 
(1966). 

We wish to discuss the relation between two recent thermodynamic 
treatments of the segregation of solute atoms at stacking faults in 
binary alloys. In one of these treatments, the equation relating the 
concentration of segregated solute to the stacking-fault energy is 
derived by means of an application of the Gibbs adsorption equation, 
while in the other treatment, the segregation equation is obtained 
by minimizing an expression for the free energy of the crystal, in- 
cluding the faulted region. The latter method of obtaining the 
segregation equation was used earlier by Suzuki. The segregation 
equations obtained by the two methods are different in form. We 
wish to point out, however, that this difference is only apparent and 
is due in part to a disguised difference in the way the stacking-fault 
energy is defined in the two treatments. We shall show that the 
possibility of such an ambiguity in the definition of stacking-tault 
energy arises in the case of an alloy but not in a pure metal. 


Key words: Stacking-fault energy in binary alloys, thermodynamic 
treatments of the segregation of solute atoms in stacking faults. 


Relaxations in polyethylene: Orientation of the tamellar 
crystals, R. K. Eby and J. P. Colson, J. Ac ist. Soc. Am. 39, No. 3, 
506-510 (Mar. 19, 1966). 

A torsien pendulum has been used to m nechanical measure- 
ments on linear polyethylene between —’ » and +110° C at about 
one Hz. Samples which had been orie ed by crystallization in a 
temperature gradient were measured t — ther with unoriented but 
otherwise equivalent samples. The d are analyzed to show that 
the @ relaxation in polyethylene is se’ ve to the orientation of the 
deformation with respect to the stru ‘symmetry. For shearing 
in the a, c face of the unit cell an sociated lamella boundaries 
the average over all direction is dif' nt from the average for shear- 
ing in all possible planes in the ce d boundaries. The y relaxa- 
tion does not exhibit this differ e. The data show a lack of 
correlation between lamella th «ness and various viscoelastic 
parameters such as the real par: oi the modulus and the magnitude 
of the relaxations. The results show the need for further mechanical 
measurements on samples that are more nearly perfectly oriented 
and more fully characterized with respect to structural details 
and imperfections. 


Key words: a, c face of cell, compliance, crystallization-oriented, 
lamella thickness, logarithmic-decrement, modulus, polyethylene, 
relaxations, shear, temperature, torsion pendulum. 


The optical constants of some oxide glasses in the strong 
absorption region, G. W. Cleek, Appl. Opt. 5, 771-775 (May 
1966). 

The optical constants, n and k, have been calculated for fused silica, 
three silicate glasses and a phosphate glass from reflectance meas- 
ments made in the strong absorption region, from about 7 to 13 um. 
Plots of the 2nkv product, which corresponds to the conductivity, 
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as a function of frequency are used to determine the location of the 
resonance frequencies in the infrared. 


Key words: Fused silica, infrared absorption, infrared reflection, 
optical constants, oxide glasses, phosphate glass, silicate glasses. 


Use of ellipsometry in the study of corrosion, J. Kruger, 
Corrosion 22, No. 4, 88-97 (Apr. 1966). 

Ellipsometry by determining the change in the state of polarization 
(ellipticity) of polarized light reflected from a metal surface bearing 
a thin film enables one to measure the thickness of this film and its 
optical properties. Because corrosion processes are strongly 
affected by the presence of thin films, ellipsometry is an especially 
valuable tool in corrosion research because it allows one to observe 
film formation, dissolution, or property changes while these phenom- 
ena are actually occurring on a metal surface immersed in any 
transparent environment. Experimentally the use of ellipsometry 
mainly requires an experimental cell that allows one to reflect light 
from the corroding or passive metal surface under study. Thus 
it does not interfere with other techniques, such as, for example, 
ultra-high vacuum surface preparation or potentiostatic studies, 
allowing them to be carried out simultaneously along with ellipso- 
metric measurements. 

Some of the experimental complications that must be considered in 
applying ellipsometry to corrosion studies are surface roughening, 
changes in the optical properties of the solution, non-uniformity 
of films, effect of light on the corrosion process, and phase changes 
in the corrosion product films. 

The following specific examples serve as a small sampling of the 
many ways in which ellipsometry may be applied to corrosion studies: 
(1) The determination of the nature and rate of growth of passive 
films on iron in inorganic inhibitor solutions. (2) The separation of 
the different stages in the early growth process of films formed by 
potentiostatic anodic polarization on iron in neutral solutions. 
(3) The determination of the point in the corrosion of copper when 
a new phase begins to form in the corrosion product film. 


Wave propagation in a three-element linear spring and dash- 
pot model filament, J. C. Smith, J. Appl. Phys. 37, No. 4, 1697- 
1704 (Mar. 15, 1966) 

Mathematical expressions are derived for the particle velocity, stress, 
and strain distributions in the wave that results when a semi-infinite 
viscoelastic filament is subjected at one end to constant velocity 
tensile impact. The equations governing the stress-strain-time 
behavior of the filament are assumed to be those for a linear model 
consisting of a spring coupled in parallel with a spring and dashpot 
in series. These equations are so formulated that by changing a 
single parameter either the spring branch or spring and dashpot 
branch can be made to dominate. The general solution, and special 
solutions describing the behavior at the wave front and at the point 
of impact are derived. Approximate solutions for rough calculations 
are also given. The application of these solutions in the interpreta- 
tion of experimental data is discussed. 


Perturbation theorems for waveguide junction, with applica- 
tions, D. M. Kerns and W. T. Grandy, IEEE Trans. Microwave 
Theor. Tech. MTT-14, No. 12, 85-92 (Feb. 1966). 

Perturbation or “compensation” theorems are derived for wave- 
guide junctions. These express changes in immittance-matrix 
elements, due to changes in the junction, in terms of integrals over 
products of changed and unchanged basis-fields associated with the 
junction and its adjoint. Media involved are required only to be 
linear. 

Applications are made to half-round obstacles in rectangular wave- 
guide, to finite conductivity in the obstacle problem, and to the junc- 
tion of rectangular with filleted waveguide. It is observed that the 
common expedient of approximating perturbed fields by unper- 
turbed ones does not always lead to correct lowest-order results. 
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